0 ratings0% found this document useful (0 votes) 46 views2 pagesNSM (Module 1) 12
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content,
claim it here.
Available Formats
Download as PDF or read online on Scribd
10.11. RUNGE-KUTTA METHOD FOR SECOND ORDER DIFFERENTIAL EQUATION
Consider the second order differential equation
a’. dy ,
oa = ds % 3 5 Wig) = Vo 39" Ho) =o
dy dy _de
Let aethen S3= 5
Substituting in (1), we get
dz
Fe HI 2) 5 Ht) = H9 izle) = Zp
The problem reduces to solving the simultaneous equations
ae =? = f(x,y, 2)
lz
and "ise f(x, y, z) subject to y(x) = Yo 3 2(%p) = Zp
and this can be solved as in the previous article.
Scanned with CamScanner em ttt teeetetee
senA_ SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS ind
Bxample. Solve y" = ay’ ~y:¥(0) =3, 9) = 0 tw approximate 90.1)
Sol. Given YP may’ ys 90) = 8, 710) = 0 a
Lety’'= thon yt=
‘The given equation (1) reduces to
Verh (uy.2h sexy ehiny2)
eatjet 10 (0) = 9 and 2(0) = 0 ie..39 = O.¥9 = Ie z9 = 0.
Now hy = Whe, You f9) = zy = (0.110) = 0
1 = HE, Yo, 29) = Mzgey —9—) = - 08
0.15)
2
= .1K-0.15) =- 0.015
Vi (20+ v0 ttsey ft) = (0.05, 3,
115)
) uy (20+ Bay Azo +44) = (005, 3,
= (0.1){(0.05\-- 0.16) — 3] = 0.30075
150375)
Wi (sean of ne) =m .005,20005,
= (0.1(- 0.150875) =~ 0.0150875
(cas Baae
vase) «170,208,030
= (0.11(0.05,-- 0.150875) - 2.9925) = —0.3000018
by = hea + h Yq + hy 29+ 1) = hf101, 2,9849625, — 0.000018)
(0.1, 0.800018) = ~0.0300001
Wflig + hy Yo + hy 24 +1) = M{0.1, 2.9849625, ~ 0.300018)
= (0.1(0.1%- 0.000018) ~ 2.9849625) = - 0.014962
1
= thy + Bk, + ky +h)
3 o+2- 001) 2¢ ont) —o0eto01 «0.150125
1
reBujeayeayery
«09-2 0s +21-0810016.-030148n = -0sto09
= (0.1) = yp +h = 8 ~ 0.0150125 = 2,9849875,
2, =2(0,1) = 5 +1 = 0 ~0,3004999 = ~ 0,9004999,
fe
Scanned with CamScanner
