Hi everybody,
I need to solve a System of Ordinary Differential Equations (ODE) by Runge-Kutta-Fehlberg method. I know that algorithm only uses with an ODE which has a function with form y'=f(x,y). Now I want to enhance this algorithm with a system of ODE, [y1',y2',...,yn']=[f1(x,y1,y2,...yn),f2(x,y1,y2,...yn),...,fn(x,y1,y2,...yn)]
How to do this? Hope you support it to me. Thanks.
What I have tried:
I tried with this algorithm for an ODE solving. But not know how to apply it for a System of ODE.