0% found this document useful (0 votes)
23 views198 pages

Unit 4

Uploaded by

prempsm2005
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
23 views198 pages

Unit 4

Uploaded by

prempsm2005
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 198

Unit 4 - Numerical Solution of initial and boundary value problems

When you have completed this unit you will be able to

1. Describe and apply various Euler methods to solve first order


ordinary differential equations.
2. Describe and apply Runge-Kutta method of fourth order to
solve first order ordinary differential equations, simultaneous
first order differential equations and second order differential
equations.
3. Classify the partial differential equation of second order as
elliptic, parabolic and hyperbolic.
4. Solve the two dimensional Laplace equation in Cartesian
coordinates.
5. Solve the one dimensioanl heat equation and one dimensional
wave equation numerically.
INTRODUCTION TO ORDINARY DIFFERENTIAL EQUATIONS
INITIAL-VALUE AND BOUNDARY-VALUE PROBLEMS
SINGLE-STEP AND MULTI-STEP METHODS
COMPARISON OF SINGLE-STEP AND MULTI-STEP METHODS

NUMERICAL METHODS OF SOLUTION OF ODE


EULER'S METHOD
IMPROVED EULER METHOD
MODIFIED EULER METHOD
RUNGE-KUTTA METHODS
FIRST ORDER RUNGE-KUTTA METHOD
SECOND ORDER RUNGE-KUTTA METHOD(Modified Euler Method)
THIRD ORDER RUNGE-KUTTA METHOD
FOURTH ORDER RUNGE-KUTTA METHOD
RUNGE-KUTTA METHOD FOR SIMULTANEOUS FIRST ORDER DE
RUNGE-KUTTA METHOD FOR SECOND ORDER DE
6/3/2021 Swinging Pendulum Problem - ODE - Appln.png

https://classroom.google.com/c/MjgwMDEzMjk0MDkw/m/MzUzMDI3MDcxODE2/details 1/2
6/3/2021 Swinging Pendulum Problem - ODE - Appln.png

https://classroom.google.com/c/MjgwMDEzMjk0MDkw/m/MzUzMDI3MDcxODE2/details 2/2
¦§¨©ª©«¬­¦§®¯°±²³ÿÿµ¶®·¶©¸ÿÿª®¹·«ÿº»ÿÿº¨¼½ÿ¯®¬

012 3456ÿ589 65ÿ  4ÿ 6684ÿ6948 5

bvd–—2+5˜5™,*ÿ b;d–—*+4˜5™62ÿ b?d–—*+1˜2™*šÿ


{GILSCÿ10@|›
~ÿœˆŽ’‰“ÿ‚‘‚Ž‚Œ‰Œˆÿ€ÿŒ‚ÿ‰Œ’‰“ÿ‹ †‚Žÿ€ÿŽ”ˆ“‰Œˆÿ€ÿŒ‚ÿž‚Œ Ÿÿ€ÿˆŠ•ÿ¡¢•£¤•ÿ¥ Œ‚
Œ‰ŒÿŒ‚ÿ†ˆ‰‹‚Œ‚Žÿ€ÿŒ‚ÿ”ˆ”’“‰ÿ‰Žÿ‰‚ÿ‘‘ Œˆ‰“ÿŒÿŒ‚ÿ‹‰ŠˆŒ’†‚Žÿ€ÿŒ‚ÿ”’‚ŒŽÿ€
‚‰”ÿ“ ‘•

ÿÿÿÿ!"ÿ"#$ÿ"%"&'"ÿÿ()*+,,,ÿ#ÿ-"ÿ'-!!'!"$ÿ$ÿ!"
"'!ÿ"&'!"$ÿ!ÿ."$
/*012*13*)34)4045637
8ÿ$"9#!"$ÿÿ:%+ÿ41+*,;ÿ!"ÿ!ÿ9ÿ#!"ÿÿ!"ÿ99!"ÿ$"#!ÿÿ!"ÿ"#<
$ÿ$ÿ!$+ÿ:.ÿ!"ÿ!$ÿ$"ÿ#ÿ-"ÿ$"!""$ÿ
/2055,3*)34)/2016*137
=">'"!.ÿÿÿ:%+ÿ41+*,?ÿ!"ÿ!ÿ$ÿ!$ÿ9ÿ#!"ÿÿ!"ÿ99!"ÿ$"#!
ÿ!"ÿ"#$+
10@2 ABCÿEFGHIGHIÿJCHKLMLNÿONCPBQHGPQMRQCSTEJQPC
CHIGHCCSGHIU
VWXYZ[\]^_` a"##ÿ"%""ÿbÿc"ÿÿÿ!"ÿ"%""dÿ"ÿ">'"!.ÿ#"$
c!ÿ9-"ÿ##"%ÿ!"ÿ9"$#ÿ!ÿÿ""ÿ-$"+ÿe"ÿ"%""%ÿ99#ÿ!
'#ÿ9-"ÿ'!!".ÿ">'"ÿ!!ÿ!"ÿ9!ÿ$ÿ&"#!.ÿÿ!"ÿ-$.ÿ-"ÿfcÿÿ
'#!ÿÿ!"+ÿe""ÿ'#!ÿÿ!"ÿ&-.ÿ"ÿ!"ÿ'!ÿÿ$.ÿ$""<
!ÿ">'!+ÿe"ÿ$""!ÿ">'!ÿ"ÿ''.ÿ-"$ÿÿg"c!hÿcÿÿ!+
{GILSCÿ10@|} 8ÿÿ"i9"ÿ#$"ÿ!"ÿ9"ÿ9"$''ÿcÿ9"&'.ÿÿ:%+ÿje6+*+ÿe"
~ÿ€‚‚ƒ„ †‡ÿ†ˆ‰Š‰‹ÿ€ÿŒ‚ 9 !#"ÿÿc"%!ÿkÿ'9"$"$ÿÿÿc"%!"ÿ$ÿÿ"%!ÿlme"ÿ.ÿ#"ÿ#!%ÿ
!"ÿ9!#"ÿ"ÿ!ÿc"%!ÿ$ÿ!"ÿ!"ÿnÿ!"ÿ$+ÿe"ÿ9!ÿÿ!"ÿ9!#"ÿ!ÿ.
ŽˆŠˆŠÿ‘‚†’“’‹ÿŽˆŠ !
Œ‚€”‚ŽÿÿŒ‚ÿ‘‰Œˆ”“‚ÿ‰† "ÿÿ#9"!".ÿ9"#"$ÿÿ!"ÿÿ!"ÿ%"ÿo$ÿlm
Œ‚ÿ‰””‚“‚‰Œˆ• e"ÿ""<-$.ÿ$%ÿÿ:%+ÿ41+*pÿcÿ!"ÿ#"ÿÿ!"ÿ9!#"ÿ$ÿ!"ÿ##""<
n
 !
 +ÿ!ÿÿ#&""!ÿ!ÿ99.ÿg"c!hÿcÿÿ!ÿÿ!"ÿq$"#!ÿ!%"!ÿ!ÿ!"
9!ÿÿ!"ÿ9!#"r
q st)/ko) kuv
v c""ÿu)!"ÿ%&!!ÿ#!!ÿb74+4ÿ!w4dÿ$ÿv)!"ÿ##""!ÿÿ!"ÿq$"#!+
y e"ÿ%'ÿ##""!ÿÿ!"ÿ9!#"ÿbxdÿ-"#"
k z x) vl
bcdefeghibcjklmnoÿÿqrjsretÿÿfjusgÿvwÿÿvdxyÿkjs

0123 456ÿ89 ÿ6  


ÿÿ!ÿ"#!$%ÿ&'()*+,)-*.
/0%+( 01 '(
2 2 )-*
01)*+

)*+3 2%+(4 &*5678.
)-* 1
%ÿ!!$9ÿ%:ÿ;<!$ÿ%ÿ!ÿ%"#=#ÿ!ÿ#$!ÿ;<!$6ÿ>ÿ?=
!ÿ%<"ÿ;<!$%ÿ!ÿ#"<$ÿÿ:%%@ÿ$ÿ%Aÿ!!9$"!96ÿB<ÿ!Aÿ$Cÿ""%
?!#?ÿ<$ÿ?%%6ÿD%$ÿ$ÿ#$!ÿ;<!$ÿ:?$ÿ@ÿ#<"#ÿ$ÿ!ÿ:ÿ$!$
"!ÿ@ÿ%A#ÿ!!9$"!9ÿ&"!ÿE"6ÿFG6767.ÿÿ%"#ÿ!ÿ<:"!ÿ!H:!$
$";<ÿ"!ÿ@ÿ<%#ÿ$ÿ%Aÿ$ÿ#$!ÿ;<!$ÿ#"$96ÿIÿCÿH!:ÿ@$ÿ
$%ÿ!$!$A%ÿÿ$%ÿH!:6
JKLMNOKPQ F"#?ÿC$ÿ$ÿ%$ÿ!!"ÿCÿ$ÿ$!$ÿ$ÿ%%ÿH!%ÿÿ%ÿ+
%ÿ?Aÿ@9
%+(+/ +RS3 +8S/ +GS3TTT
R 8 G
&*567U.
Dÿ%:!ÿ!?<!ÿ#%!":$%ÿ%ÿ+%ÿ!H:!$9ÿ;<!ÿ$ÿ+CÿH%%#ÿÿ!#=
!%6ÿÿÿ%:!ÿ#%!":$%ÿV;6ÿ&*5678.ÿ@":%
)*+3 2+(4 &*567G.
)-* 1
C"ÿ%ÿ!ÿ%"#=#ÿ!ÿ#$!ÿ;<!$6ÿ%ÿ!H:!$ÿ%ÿA9ÿ:$!$
@"!<%ÿV;6ÿ&*567G.ÿ%ÿ!%9ÿ$ÿ%Aÿ!!9$"!96ÿÿ%<$ÿ@!%#ÿÿ$ÿ$9ÿÿ#=
$!ÿ;<!$%ÿ%ÿ?Aÿ@9
Y
+W-X(+4"% 21- &*5675.
Cÿ+4($ÿ#%!":$ÿ!$ÿ-(4ÿ!#ÿCÿ$ÿ%ÿ!%%<:#ÿ$!$ÿ$ÿA"$9ÿ&Z()+,)-.
%ÿ[ÿ!$ÿ-(46ÿÿ$:ÿ;<#ÿÿ$ÿ#<<:ÿ$ÿ":$ÿÿ"9"ÿÿ%"!$ÿ%
"!#ÿ$ÿ#ÿ!#ÿ%ÿ?Aÿ@9
^
\(*] 21 &*567_.
D?<ÿ*567Gÿ%C%ÿ!ÿ$ÿÿ$ÿ#%!":$ÿ+!#ÿA"$9ÿ)+,)-!%ÿ!ÿ<"$ÿ
$:ÿ!%ÿ"!"<!$#ÿ:ÿV;6ÿ&*5675.ÿC$ÿ+4(],`ÿ!#ÿ1(*ÿ$6ÿÿ#ÿ!%ÿ"!"<!$#
:ÿV;6ÿ&*567_.ÿ%ÿ768U8_%6
!@A"!"<!$%%%$!9!!":$%<$$:$$#<=
<:6ÿaCA9<:<%$!%"%#$!""<!"9$%<$%@"!<%$!%%<:$%
„ †‡ˆ‡‰Š‹„ ŒŽ‘ÿÿ“”Œ•”‡–ÿÿˆŒ—•‰ÿ˜™ÿÿ˜†š›ÿŒ‰

012 3456ÿ589 65ÿ  4ÿ 6684ÿ6948 5

~ >~
>@
€‚
ƒ

€ @ € @
PQRSTUÿ10WXY
Z[\]ÿ\_ÿ`abc[defgfh]ÿ?dh` ƒ
if[\ea]jÿk?lkmdbÿdÿ_nhe]a\hÿ\_ €‚
]agfÿmoÿdbÿed[en[d]f`ÿ_p\g
qrstuvswvxsÿ?yabÿzl{ÿdh`ÿ]|f
[fh}]|ÿabÿuÿ_]s

ÿÿÿ!"#ÿ$%ÿ&'()'ÿÿ'**)'*+,ÿÿ-ÿ*--'+ÿ%ÿ%.'ÿ'ÿ)/*'(
-%()%ÿ0%ÿÿ!1#,ÿ2*ÿ-ÿ'ÿ/%ÿ*%/3(ÿ3+-*'(ÿ3-'%ÿ%0ÿÿ/%%
4+ÿ%0ÿÿ/%5-ÿ5-*)--5ÿÿ6'3-ÿ1ÿ'5ÿ7ÿ*%)(5ÿ.ÿ)-5ÿ0%ÿ-ÿ3)3%-80%
9'/3(,ÿÿ)(ÿ'5ÿ0%):%5ÿ;<ÿ/%5-ÿ)'%ÿ!1#ÿ/)-ÿ.ÿ'-0%/5
%ÿ2%ÿ0-:%5ÿ)'%-ÿ%ÿ.ÿ*%/3'.(ÿ2ÿÿ'.%&ÿ/%5-ÿ$-ÿ-ÿ'**%/:
3(-5ÿ'-ÿ0%((%2-ÿ$ÿ&(%*+ÿ=-ÿ505ÿ.+
>?A= B#
>@
'5,ÿ0%,ÿÿ!1#ÿ*'ÿ.ÿ93--5ÿ'-
>=ACD-? !#
>@ E
)'%-ÿB#ÿ'5ÿ!#ÿ'ÿ'ÿ*%)3(5ÿ-+-/ÿ%0ÿ2%ÿ%5'+ÿ500'(ÿ)'%-
$ÿ)/*'(ÿ-%()%-ÿ.+ÿÿ)(ÿ/%5ÿ'5ÿÿ0%):%5ÿ;<ÿ/%5ÿF&ÿ
-)(-ÿ-%2ÿÿ$'.(ÿ!,ÿ2*ÿ*%/3'-ÿÿ''(+*ÿ-%()%ÿ0%ÿÿ('ÿ)'%
%0/%%ÿGÿ!#Hÿÿ*%()/ÿI#ÿ2ÿÿ)/*'(ÿ-%()%-ÿÿ*%()/-ÿJ#,ÿK#,
'5>#
$ÿ)(ÿ'5ÿ0%):%5ÿ;<ÿ/%5-ÿ+(5ÿ500ÿ-)(-ÿ'5ÿ.%ÿ5-'Fÿ2
ÿ''(+*ÿ-%()%,ÿ'(%)Fÿÿ0%):%5ÿ;<ÿ/%5ÿ0%ÿÿ%('ÿ*'-ÿ-ÿ*(%-
%ÿÿ''(+*ÿ-%()%ÿ'ÿ-ÿÿ)(ÿ/%5ÿ$%ÿ3%3(+ÿ&'()'ÿÿ500*ÿ.:
2ÿÿ('ÿ'5ÿ%('ÿ/%5(-,ÿÿ-ÿ/3%'ÿ%ÿ5/ÿÿ'**)'*+ÿ%0ÿÿ):
/*'(ÿ-)(-ÿ$-ÿ-ÿ'**%/3(-5ÿÿÿ2'+-ÿL-,ÿÿ)(ÿ)/*'(ÿ-%()%ÿ-
'-(+ÿ*%FM5ÿ'-ÿ'5)'ÿ.*')-ÿÿ%&-%%-ÿÿ'(ÿ*%5%ÿ'ÿ@ABÿ-ÿ$-
*('(+ÿ&%('-ÿ*%-&'%ÿ%0ÿF+ÿN*%5,ÿ*%()/ÿK#ÿ'5ÿ>#ÿÿ$'.(ÿ!ÿ-%2ÿ
-%()%ÿ%0ÿÿ0%):%5ÿ;<ÿ/%5ÿ0%ÿ-3ÿ-M-ÿ%0ÿBB1ÿ'5ÿBB!ÿO*')-ÿ-ÿ&'+
¢£¤¥¦¥§¨©¢£ª«¬­®¯ÿÿ±²ª³²¥´ÿÿ¦ªµ³§ÿ¶·ÿÿ¶¤¸¹ÿ«ªº

0123 456ÿ89 ÿ6  


DEFGHÿJ KLMNOPQRLSÿLUÿOÿ[QSXOPÿOSO[^VQZO[ÿRL[_VQLSÿLUÿVWXÿR“QS\QS\ÿNXSY_[_MÿNPL][XM
“QVWÿVWPXXÿSLS[QSXOPÿS_MXPQZO[ÿRL[_VQLSR`
oejmdjblcÿotbcdslmÿ‚emkdejq
Gdjblc
Ejlm”kdslm Hmbc •k–—˜cfbcÿ™š •k–—˜cfbcÿ™š
Ddtbg ‚emkdej v›xuJu† v›xuJu† v›xuJu 
q vœ v vž vŸ
‹ ‹Ž’Œ’ ‹Ž’Œ’ ‹Ž’Œ’ ‹Ž’Œ’
‹‘ ‹ŒŒŽ’ ‹Š‰ŒŒ ‹ŒŠŠŒ’‘ ‹ŒŠŠŒŽ
‹ ¡ ‹‘Š’Œ‘ ‹Œ ‘‘’ ‹‘‰’Œ ‹‘‰’’‘
‹Š ¡ ‹Œ’‰  ¡ ‹ŠŠŒ‘ ¡ ‹ŒŒ’ ‘ ¡ ‹ŒŒ’‘
‹’ ¡ ‹Ž’ŒŠ‘ ¡‰‹Œ ŠŽ ¡ ‹Ž’‘Š ¡ ‹Ž’‘‘
‰‹ ¡ ‹Œ Œ‰‘ ¡ ‹ Š‘‘ ¡ ‹ŒŒŒ’ ¡ ‹ŒŒŒ’
‰‹‘ ‹’ ‰ ¡ ‹‘’‰ ¡ ‹ŠŒŠ‰‰ ¡ ‹ŠŒŒŽŒ
‰‹ ‹Š‰ŽŠ’ ‹Š‘‰Ž ‹Œ Œ‘ ‹Œ ‘
‰‹Š ‹ŽŽ’ Š‘ ‰‹‰ŠŽŒ ‹Ž’ ŽŠ‘ ‹Ž’ ŽŽŽ

DEFGHÿJKLMNOPQRLSÿLUÿVWXÿNXPQLYÿLUÿOSÿLRZQ[[OVQS\ÿ]LY^ÿZO[Z_[OVXYÿUPLMÿ[QSXOPÿOSY
SLS[QSXOPÿMLYX[R`
abcdefgÿh
ijdkdlm Gdjblcÿnefbm oejmdjblcÿnefbm
pdqrmlsbtbjkgÿ1u vwxyyz~
{|} €otbcdslmÿ‚emkdejÿeƒÿH„JÿvJ†‡
yˆ‰Š ‰‹ŒŠŒ ‰‹ŒŽ
yˆ ‰‹ŒŠŒ ‰‹Š
yˆ‘ ‰‹ŒŠŒ ‰‹’Œ

ÿÿÿ!"#ÿ$#"%ÿÿ&ÿ"&"'#ÿÿ"&&!ÿ"ÿÿ&#ÿ(ÿ"ÿ&$ÿ&)
*+*,ÿ&ÿ"#&ÿ" "ÿ(ÿ&ÿ-ÿÿ".+ÿ/%ÿÿÿ*+*,0&ÿ&$0&)ÿ"&%
1'"&ÿ"ÿ#"#ÿ!"2!!ÿ3"#ÿÿ*+456756ÿ"ÿ89,+:7ÿ&ÿ;ÿ&(ÿÿ/"'#ÿ<5+,=+ÿ/&
"&ÿ"ÿÿ$#!ÿ&ÿÿ&ÿ-"#ÿ$&ÿ(ÿ0$#"ÿ" ".ÿ(ÿ"
$ÿÿ,+:7ÿ&+ÿ/&ÿ&"&ÿ"##(ÿ.ÿÿ&"#.ÿ"&&!ÿ"ÿÿ
'(ÿ#!&ÿ;>=ÿ"ÿ;?=ÿÿ/"'#ÿ<5+,ÿ#.ÿ$&&ÿÿÿ'(ÿ
#"ÿ"ÿ#"ÿ!#+
@ÿ(".ÿÿ"")ÿÿÿ'(ÿÿ#"ÿ"ÿÿ#"ÿ!#
&ÿÿÿ'"&&ÿÿ$+ÿ/"'#ÿ<5+<ÿ&(&ÿÿ$ÿÿ&##"ÿ"&ÿ"##"ÿ'.ÿÿ#0
"ÿ!#ÿ"ÿ#"ÿ!#ÿÿÿÿ"#ÿ&$#"!&+ÿAÿ&ÿ&ÿ"ÿ
"##"ÿ$&ÿ"-ÿ#&#.ÿ(ÿ1&ÿ&!"##ÿ'"&ÿ1&ÿ"ÿ-ÿ"$$2!"ÿÿ&1
ÿBC+ÿ;<5+,:=+ÿ/&ÿ"$$2!"ÿ"&ÿ(ÿ1'!&ÿ#"-+
/&ÿ""#.&&ÿ"ÿ.$"#ÿÿ"&&ÿ.ÿ(##ÿ#.ÿÿ"&ÿ"ÿ-+ÿ/
#.ÿÿÿ!"#ÿC&ÿ'!&ÿ$"#"#.ÿ&-"ÿÿ#"ÿ$'#!&%
"ÿÿ!".ÿ"&&ÿ"#0(#ÿ$'#!&ÿ"ÿ#"+

You might also like