Hi Everybody,

Suppose I have this procedure assign to p:

>p:=MyProc(args1,args2,args3,h) statement end proc:

where h is the step size.

args1 = {ode,ic} , args2 = x(t) et args3 = t=a..b.  The output of MyProc is a list of list

I would like to do what dsolve do when you write:

p:=dsolve( {ode,ic},x(t),x(t),numeric)

and then evaluate p(2) for example.

Hi all,

Suppose that you want to solve x'(t)=f(x,t)   x(t0)=x0

The theory show that:

x(t+h)=x(t)+1/2*h*(F1+F2) where F1=f(x,t) and F2=f(x+F1,t+h)

SO I made procedure:

Hi all,

while studying this ODE

Hi all,

> deq := diff(theta(t), t, t) = -g*sin(theta(t))/L;

                           2                              

Hi all,

While solving a différential equation, I encounter an integral: