%y = dsolve('Dy=.5*y-y*y*.5/1.44','y(0)=.1');


for i=1:60
t(i)=(i-1)*.25;
z(i)=36/(25+335*exp(-.5*t(i)));
r(i)=sqrt(2*(1-exp(-.5*t(i))));
end
plot(t,z,'r',t,r,'c')
hold on;
for i=1:21
    tt(i)=(i-1)*.25;
    zz(i)=.1*exp(.5*t(i));
end
    plot(tt,zz,'b');
    for i=1:60
        ttt(i)=(i-1)*.25;
        zzz(i)=1;
    end
    plot(ttt ,zzz,'g');