>

>

> with(DEtools);

> DE:=diff(y(t),t)=y(t)^2+t^2;

> DEplot(DE,y(t),t=-2..2,[[y(0)=1],[y(1)=-.5]],y=-2..2);

> dsolve(DE,y(t));

> DE1:=diff(y(t),t)=2*y(t)/t;

> DEplot(DE1,y(t),t=.1..5,[[y(.1)=0.1],[y(1)=-.5]],y=-5..5);

> dsolve(DE1,y(t));

> DE2:=diff(y(t),t)=t-3*y(t);

> DEplot(DE2,y(t),t=-2..2,[[y(0)=1],[y(1)=-.5]],y=-2..2);

> dsolve(DE2,y(t));

> DE3:=diff(y(t),t)=2*y(t)*(1-y(t)/3);

> DEplot(DE3,y(t),t=0..5,[[y(0)=.5],[y(0)=2]],y=0..4);

> dsolve(DE3,y(t));

>