Dear all,
Some time ago I asked a question concerning the algorythm for numerical dsolve and related issues here.
I got a nice answer from pagan and implemented the algorythm.
However, since my system is big, the iterations take long time.
Here is the code itself:
sys1:=[TRIC[1],TRIC[2],TRIC[3],TRIC[7],TRIC[8],TRIC[9]]:
> sys2:=[TRIC[4],TRIC[5],TRIC[6],TRIC[10],TRIC[11],TRIC[12]]:
> GRAD:=subs(u=FF[1](t),a=FF[2](t),k=k(t),tau=tau(t),psi[1...