Example 3: x*y

> restart; with(DEtools): with(plots):

> de3:=diff(y(x),x) = x*y(x)/4;

> L:=dfieldplot(de3,y(x),x=-2..4,y=-2..8,arrows=medium):

> display({L});

> dsolve({de3,y(0)=1},y(x));

> y0 := x -> exp(1/8*x^2):

> dsolve({de3,y(-1)=4},y(x));

> y1 := x -> 4*exp(1/8*x^2)/exp(1/8):

> M0:=plot(y0(x),x=-2..4,colour=blue,thickness=3):

> M1:=plot(y1(x),x=-2..4,colour=yellow,thickness=3):

> display({L,M0,M1});

>