- 5

- 0

## Main Question or Discussion Point

Hello All,

I need to utilize a Runge Kutta second order approach to solve two coupled first order DE's simultaneously given some initial conditions and a conservation relationship.

The DE's are as follows:

[tex]\frac{dp}{dt}[/tex] = aq - bp

[tex]\frac{dq}{dt}[/tex] = -aq + bp

Where a and b are constants, a=3, b=4, and the conservation relationship is n

I'm having trouble getting the ball rolling with writing the R.K. algorithm for this, will be programming in matlab. I have solved the equations analytically by hand and have found the following:

p(t) = -0.055972e

q(t) = -0.055972e

Any help on how to get started with the Runge Kutta algorithm to solve numerically would be appreciated!

I need to utilize a Runge Kutta second order approach to solve two coupled first order DE's simultaneously given some initial conditions and a conservation relationship.

The DE's are as follows:

[tex]\frac{dp}{dt}[/tex] = aq - bp

[tex]\frac{dq}{dt}[/tex] = -aq + bp

Where a and b are constants, a=3, b=4, and the conservation relationship is n

_{t}= p(t) + q(t). Here, n_{t}= 3.I'm having trouble getting the ball rolling with writing the R.K. algorithm for this, will be programming in matlab. I have solved the equations analytically by hand and have found the following:

p(t) = -0.055972e

^{-3t}+ 3.055972e^{-4}q(t) = -0.055972e

^{-7t}+ 3.055972Any help on how to get started with the Runge Kutta algorithm to solve numerically would be appreciated!