g=0.2;
m=0.6;
delt=0.00001;
xo=0.00001;
yo=0.00001;
zo=0.00001;

for n=1:1:1000000
   xtpt=xo+(m*xo+yo-xo*zo)*delt;
   ytpt=yo-xo*delt;
   ztpt=zo-(g*zo+g*0.5*xo*(xo+abs(xo)))*delt;
   xo=xtpt;
   yo=ytpt;
   zo=ztpt;
   disp(n);
end
a=1;
sech=0;
for n=1:1:1000000
   xtpt=xo+(m*xo+yo-xo*zo)*delt;
   ytpt=yo-xo*delt;
   ztpt=zo-(g*zo+g*0.5*xo*(xo+abs(xo)))*delt;
   if(abs(itog(n,2))<delt&ytpt>0)
   sech(a,1)=ytpt;
   a=a+1;
   end
   itog(n,1)=xtpt;
   itog(n,2)=ytpt;
   itog(n,3)=ztpt;
   xo=xtpt;
   yo=ytpt;
   zo=ztpt;
   disp(n);

end
plot(itog(:,1),itog(:,2),'r.');