- #1

- 7

- 0

## Main Question or Discussion Point

Hi there. I want to evolve a system of non linear coupled ODEs

[tex] \frac{dx}{dt} = \frac{-k}{x^5(56-y^8)^{9/2}}(85+y^{5} + y^{6}) [/tex]

[tex] \frac{dy}{dt} = \frac{-k}{x^4(56-y^5)^{7/2}}(44+y^2) [/tex]

Let's say I have the initial conditions. What numerical method someone could use to solve this? adaptive step size like RK dormand prince method? I am a bit confused cause i tried some with no good results (maybe cause I am not good programmer, not of the method's fault)

[tex] \frac{dx}{dt} = \frac{-k}{x^5(56-y^8)^{9/2}}(85+y^{5} + y^{6}) [/tex]

[tex] \frac{dy}{dt} = \frac{-k}{x^4(56-y^5)^{7/2}}(44+y^2) [/tex]

Let's say I have the initial conditions. What numerical method someone could use to solve this? adaptive step size like RK dormand prince method? I am a bit confused cause i tried some with no good results (maybe cause I am not good programmer, not of the method's fault)