Solving second-order ODE with Runge-Kutta 4

  • Thread starter Thread starter LANS
  • Start date Start date
  • Tags Tags
    Ode Runge-kutta
LANS
Messages
24
Reaction score
0

Homework Statement


0l0AByW.jpg

Note: I think there is a typo here but I'm not sure. Is there supposed to be a comma between the delta t/2 and y_n on K2 and K3, and delta t and y_n on K4?

Homework Equations


See above.

The Attempt at a Solution


Substituting dy/t = z gives
\frac{dz}{dt} = 3z - 2ty - cos(t)<br /> <br /> \frac{dy}{dt} = z

I'm not sure where to go from here. I can find K_1, but I'm not sure how to find K_2z as it depends on t, y, and z. What do I choose for z in K2? Do I need to redefine K_2 as

f \left ( t_n + \frac{\Delta t}{2}, y_n + \frac{K_1_y}{2}, z_n + \frac{k_1_z}{2} \right ) Is there some other way I should approach the problem?

Any help is appreciated, thanks.
 
Physics news on Phys.org
In K1, there should be a comma between tn and yn.

In K2, K3, and K4, there should be commas between the delta t terms and the yn terms.

Remember z = dy/dt, and you are given dy/dt = 0 at t = 0.
 
SteamKing said:
In K1, there should be a comma between tn and yn.

In K2, K3, and K4, there should be commas between the delta t terms and the yn terms.

Remember z = dy/dt, and you are given dy/dt = 0 at t = 0.

This is aimed towards SteamKing, K1 is found with dy/dt=0, but how is K2 solved when f(0.05,1) is undefined by the question and there is no obvious way to determine the f(0.05,1).

Thanks
 
You are incrementing t between calculating K1 and K2, which is not what is called for in the definition of K1 ... K4. For K2, you use the same values for yn and t as for K1. Once you have calculated K2, you use this value to calculate K3, and you use K3 to calculate K4. Once you have calculated K1 ... K4, then you increment t by delta t and determine y(n+1).
 
Prove $$\int\limits_0^{\sqrt2/4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx = \frac{\pi^2}{8}.$$ Let $$I = \int\limits_0^{\sqrt 2 / 4}\frac{1}{\sqrt{x-x^2}}\arcsin\sqrt{\frac{(x-1)\left(x-1+x\sqrt{9-16x}\right)}{1-2x}} \, \mathrm dx. \tag{1}$$ The representation integral of ##\arcsin## is $$\arcsin u = \int\limits_{0}^{1} \frac{\mathrm dt}{\sqrt{1-t^2}}, \qquad 0 \leqslant u \leqslant 1.$$ Plugging identity above into ##(1)## with ##u...
Back
Top