Hi there. I want to evolve a system of non linear coupled ODEs(adsbygoogle = window.adsbygoogle || []).push({});

[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)

**Physics Forums | Science Articles, Homework Help, Discussion**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# C/++/# Numerical methods for system of ODES

Have something to add?

Draft saved
Draft deleted

**Physics Forums | Science Articles, Homework Help, Discussion**