syms x
y=dsolve('Dy=x*y','y(0)=1','x')
Y=subs(y,x,0:0.1:4);
plot(0:0.1:4,Y)