Question: Recursive procedures in t

Dear all,

I am stuck with my gradient projection algorythm. I asked several questions on this but nothing helped.

So, the basic problem I have is the following:

I need to make a procedure, which would take as an input some function of time, x(t),

calculate numerical solution for the ODE system

diff(y(t),t)=f(t,x(t),y(t))->y_1(t)

via dsolve/numeric

then update the value of the function x(t)=f(y_1(t)),

giving as a result x(t) NOT as a function of procedure being returned by dsolve/numeric, but as a function of t.

Then this new x(t) is iterated with some error bound, etc.

Does anyone have any ideas as to how I may code such an algorythm?

Thanks in advance..

Please Wait...