function X = dglv (t,y,Z)

ylag1 = Z(:,1);

gLa=3;
kgamma=0.1;
No=4;
ggamma=0.01;
kgammar=0.01;
Resc=0.05;
alpha=3;
tec=50;
k=0.15;
BN=(Resc*rho)/(1-rho);


X= [ -1*kgamma*y(1)+kgamma*sqrt(k)*e^(gLa*(1-1i*alpha)*(2*rho*(t-tec)-1))*ylag1(1);...
     -1*kgammar*rho+BN*(1-rho)-Resc*rho-(e^(2*gLa*(2*y(2)-1))-1)*abs(y(1));...
     No-ggamma-2*BN(1-rho)+2*Resc*rho]
 


end