Solving a system of ODEs using Runge-Kutta

  • Thread starter P_Sampras
  • Start date
  • #1
1
0
Using the 4th-order Runge-Kutta method, I have been able to successfully compute the solutions to a coupled pair of two first order differential equations using the following formula:

T5JedKmm.png

MBFNciXm.png


When solving systems of ODEs with more than two equations I am unsure if I am properly expanding on the above ideas.

For example I am trying to find a system of three first order equations:

da/dx = 2bc + 3a

db/dx = 4ac +11b

dc/dx = 5ab + c

Using and expanding on the equations above would I be correct in the following:

K0 = 2bc + 3a
L0 = 4ac +11b
M0 = 5ab + c

k1 = 2*(b+(L0)/2)*(c+(M0)/2) + 3*(a+(K0)/2)
L1 = 4*(a+(K0)/2)*(c+(M0)/2) + 11*(b+(L0)/2)
M1 = 5*(a+(K0)/2)*(b+(L0)/2) + (c+(M0)/2)

And so on for K2,L2,M2 etc?

Hopefully that's fairly readable.

Is this the correct approach to take or am I making a glaring error somewhere?
 

Answers and Replies

  • #2
AlephZero
Science Advisor
Homework Helper
7,002
293
That looks rignt.

The basic idea is you first calculate all the "_0" quantities, then use them to calculate all the "_1"'s, etc.

Trying to do the calculations in a different order won't work, except in special cases and if you get lucky.
 

Related Threads on Solving a system of ODEs using Runge-Kutta

  • Last Post
Replies
0
Views
8K
Replies
4
Views
7K
Replies
0
Views
2K
  • Last Post
Replies
4
Views
2K
  • Last Post
Replies
3
Views
3K
  • Last Post
Replies
5
Views
6K
  • Last Post
Replies
16
Views
2K
  • Last Post
Replies
6
Views
19K
  • Last Post
Replies
7
Views
2K
Top