> #This procedure uses the numerical method to solve the differential equation in #demo2.mws. > with(DEtools); # We firist define the initial valued problem: > diffeq:=diff(y(x),x)=x+sin(y(x)); > inits:=y(0)=1; # Then we use dsolve, but include the option "numeric": > sol:=dsolve({diffeq, inits}, y(x), numeric); # The approximation value of the solution at x=1: > sol(1); #If you want to have a picture of the solution, then do: > with(plots):odeplot(sol,[x,y(x)],-1..1); >