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
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..