The differential equation numerical analysis

On 05/06/2012 02:35 AM, Part Marty wrote:
> Please tell me the differential equation numerical analysis like this!
>
> atvalue(y(x),x=0,0);
> atvalue(z(x),x=0,1);
> s1:'diff(y(x),x)=y(x)^2+x*z(x);
> s2:'diff(z(x),x)=x^2+y(x)*z(x);
> desolve([s1,s2],[y(x),z(x)]);
sol: rk ([y^2+x*z,x^2+y*z], [y,z], [0,1], [x,0,1.5,0.01])$

plot2d ([[discrete, makelist ([sol[i][1], sol[i][2]], i, 1, length(sol))],
              [discrete, makelist ([sol[i][1], sol[i][3]], i, 1, 
length(sol))]]);

See "? rk" for details.

Regards,
Jaime