6. Example-3
Find y(0.1) for `y''=-4z-4y`, `x_0=0, y_0=0, z_0=1`, with step length 0.1 using Runge-Kutta 3 method (2nd order derivative)
Solution:
Given `y^('')=-4z-4y, y(0)=0, y'(0)=1, h=0.1, y(0.1)=?`
put `(dy)/(dx)=z` and differentiate w.r.t. x, we obtain `(d^2y)/(dx^2)=(dz)/(dx)`
We have system of equations
`(dy)/(dx)=z=f(x,y,z)`
`(dz)/(dx)=-4z-4y=g(x,y,z)`
Third order R-K method
`k_1=f(x_0,y_0,z_0)=f(0,0,1)=1`
`l_1=g(x_0,y_0,z_0)=g(0,0,1)=-4`
`k_2=f(x_0+h/2,y_0+(hk_1)/2,z_0+(hl_1)/2)=f(0.05,0.05,0.8)=0.8`
`l_2=g(x_0+h/2,y_0+(hk_1)/2,z_0+(hl_1)/2)=g(0.05,0.05,0.8)=-3.4`
`k_3=f(x_0+h,y_0+2hk_2-hk_1,z_0+2hl_2-hl_1)=f(0.1,0.06,0.72)=0.72`
`l_3=g(x_0+h,y_0+2hk_2-hk_1,z_0+2hl_2-hl_1)=g(0.1,0.06,0.72)=-3.12`
Now,
`y_1=y_0+h/6(k_1+4k_2+k_3)`
`y_1=0+0.1/6[1+4(0.8)+(0.72)]`
`y_1=0.082`
`:.y(0.1)=0.082`
`:.y(0.1)=0.082`
This material is intended as a summary. Use your textbook for detail explanation.
Any bug, improvement, feedback then
|
