TMAддз Differential Equations and Scientific Computing, part A

 !"$#%&')( *,+-.)+
/ #0 123"45 6718'9:
;<>=@?BADCE<GFIHJAK<MLND<>OP=EQSRTHVU QWQYX[Z
\)]_^a`cb>dfehgfikj>lmiiWl
!T "#%$'&)(+*,(.-0/21435768"&!9:&<;>=@?BADCFE!GIHJ?KADC!E!GIL"M
N 0O
* NQPSR C PT O= P T NQPSRVU>W0XY O[ZF\^]`_#La$Fb
cedf+d0gihDD$'&jg^lk 14mFnpoqmsrp5utF3wv * NxPSR C PT O v4y gikz_#gj_#{+&|$'8D&
}"~:d0#z$'&%#D&d€7&!8f€%+$'T d9‚d€ƒ*„~:&!T h9…%$'&+†\
N }O‚* NQPSR C PT Oƒ=‡EE XY0ˆ‰UŠP RS‹ PSRUŒP T‹Ž PT.XY \,]`_#La$Fb2dgih"z$'&}d$'8gik14mFn‘omsr‘5utF3‰* NxPSR C PT O
L"M’gj_#L 14m n oqm r 5utF3 * NQPSR C PT O.\
“ ”•—–l˜™,š
N 0OV›z_L"!&j(+*2( - / 1435 =œgik 14mFno msr‘5xt!3wv * NxPSR C PT O v 6&`6%Lž$$'d|ŸLDM^$'8"& g^lkD_g¡"g¢d€£$'8"&‰}"~…d#D$…&
%#D& v * NQPSR C PT O v d€* NxPSR C PT O.\7¤9…dg¥$'8"&V8"_#La$I6&~:$%9:$}ž†
gikz_#gj_#~:_#L"¦$…8"& PT…§ M"&!h&!LDM"&!La$7€%+$'d9
d¨0&+9V$'8D&‰_#La$…&!9:%?BADCFE!G©b7c‚8D&‰9…&+~…"q$_~2$…9…_qz_##q†^E N €Qd9 PT =ªElO.\«c‚8"&‰gikz_#g|"g¬d€S$…8"& }"~…d#D$…&
%#D&d€$…8"& PSRp§ M"&+h”&+L"M"&+Lž$`€s$'d09‚_~  X0W0­ €Qd09 PSR =®AD\ƒc‚8"_~gi&!L"~$'8"%$w(+*2( -0/21435 =  XW0­ \
N }O¯°&±!d0gjh"&+$'& $'8"&±~…²ž%9…&!~`$'dw¦0&s$Fb
* NQPSR C PT O2=ªE0E XY0ˆU°P TR ‹ PSR«UJP TT ‹³ PTsX%Y = ­0XYU´NxPSR«U E XW O T U´NxPTVU¶µžX%Y O T
¯°&‰F%L^L"d¨6´M"&s$'&!9:gi_L"&$'8"&`g^lkD_g¡"g¬}a†jgj_#L"_gi_~…_#LD¦‰$'8"&`$p6d|L"&!¦a$'_&`$…&!9:gi~2d¨0&!9;b,·¸kz_ §
g¡"g¢d%€e*’$'8ž"~Vdz+!"9:~V€Qd09 PSR =¹E X0W L"M PT =‡E 68D_#8)¦0_q0&+~"~2$'8%$Vgik 14mFn‘o msr5xt!3 * NQPSR C PT Oƒ=
­0XY,U E Xº = W A Xº \»pLj$…8"&`~'gj&‚6F†wgj_#L"_g¡"gªdz+!"9:~268"&+Li$…8"&u~¼$ƒ$p6d
$'&!9:gi~29:&g^kz_g^ y
_½K&½ y €Qd9 PSR =@A)d9 PSR =¾E¡L"M PT =@Az\±] &!LD!&¿gj_#L1ÀmFnpo msr‘5utF3‰* NxPSR C PT O= ­0XY
U E X%µ|U E ˆaXº =
E0E X%Y0ˆ \‚›_#L"+&
$…8"&<gj_#LD_#g¡"gÁ_#~hd0~…_q$'_q0& y * NxPSR C PT O= v * NQPSR C PT O v _#Lf; y %L"Mf6V&|+d0L"!"M"&±$'8$
(+*2(.-/21Â35I=œgik1ÀmFnpoqmsrp5xt!3* NQPSR C PT Oƒ= W A Xº \
ÃS!"#T %$'&)(+*,( - r 1435 68"&!9:&|;Ä=¬?BADCFE!GSHŽ?KADC!E!GIL"M
N 0O‚* NQPSR C PT Oƒ= PSRP T \
N }O`* NQPSR C PT O,=Å~:_#L RpNQʞƔËQÇSÌpÍ PS1 T‘R Î OD5 ~:_#L NÈ)ÇSPT OV6_$…8 Æ %L"M È 9:}"_$…9'9¼†)_#La$…&!¦0&+9…~!\
] _Lž$!bƒ~…_L TÉ =
T
n
T
“ ”•—–l˜™,šc‚8"&<Ï N ;O § L"d09:gÐd€«*„_#~M"&sŸL"&!Mf}a†”b(+*2( - r 143z5 = NÑaÑ 3 * NQPSR C PT O TÒ PSR Ò PT O r \
N 0O
R R T Ò Ò
R TÒ
R Ò
T
(+*2( - r 1435 =¥Ó Ó P R PÕT PSR PT =œÓ P R PSR Ó PÕT PT =¬? P ZR XY G ÔVR Ö ? P×T X0Ø G ÔR = E E Ø
Ô Ô
Ô
Ô
R
~…dÙ(+*2( - r 1Â35 = Ú R \
N }O£»‘€ Æ L"MÛ%d09 È × _~S&!²ž%ž$…d‚{!&+9…d$'8"&+L¡*¡_#~—_M"&!La$…_#F%#†`&+²"0$'d{!&+9…d`_#gjh"†_L"¦V$…8%$(+*2( - r 1Â35 =
AD\2܉$'8D&!9:6_~…&6V&<¦0&s$Fb
R R T
T
r
(+*2( - 1Â35 =®Ó Ó ~…_#L NƐÇSPSR OD~…_L T NÈ)Ç£PT O Ò PSR Ò PT
ÔR Ô
R E U !d0~ N©W%È)ÇSPT O Ò
U

N
W
”
Æ
S
Ç
S
P
R
E
!
d
~
O
Ò
PSR Ö Ó
PT
= Ó
W
W
Ô
Ô
R
R

N
W
”
Æ
S
Ç
S
P
R

N
%
W
)
È
S
Ç
P
T
…
~
_
L
O
…
~
#
_
L
O
=ÞÝ PSR…X0WU
µ0ÆÇ ß Ô Ö Ý PTsXWU
µÈ)Ç ß Ô
E
= —E X0W‰U …~ _L µ NÆ”WÇÆ”Ç O Ö —E X0W‰U ~:_#L µ0N©È)W%È)Ç Ç O =¹E X%µ C
L"M$'8ž"~¿(+*2( - r 143z5 =¹E X0W _€ Æ =œAwL"M È =œAD\
S —&s$ N O2= NQP O2= Ô ‹ R[PSR ‹ T…PT C
„C «=@E0C W C Y P = R NQPSR C T PT O
}&^$'8"& ~:h+& d€#_L"&F9<h”d†L"d0gj_u~±M"&+ŸLD&!MŽd0Lޏ$'9:_uLD¦0#& 6_$'R 8¶+d09…LDR &!9…~ R y T y LDT M žT Z¨\
&!9:_&&skzh"_#!_q$&.kzh"9…&+~…~…_d0L"~ N _#L^$…&!9…gj~d€—$…8"&+d09…L"&+9!dd09:M"_#L$'&!~ = N R C T O y = N R C T O y
N O‚MD&+ŸL"&+M¸}a†
L"M aZ= N aZR C aZT O…OV€Qd09$…8"&±}~:_#~‚€Q"L"s$'_#dL"~ R y T y Z
N !O,= E = 0C
N E¨O
A = 0C
6_$…8 :C ^=¬E0C W C Y \]`_#La$FbV~…&+$‰"h¸$…R 8"&|_#L"T &!9~:†~:$…&!g d€ƒ&!²ž$'_#dL"~`68D_#8„9:&!u$'&!~ Ô y R y L"M T
N Os\`›d0&
€Qd09$'8D&
$'dj$'8"&±%#D&!~ %$$'8"&<!d09:L"&!9:~ N O y N O y L"M N žZOd€ƒw€Q"LD+$'_d0L
!d& j+_#&!La$…~ +d09…9:&!~…hd0L"MD_#L"¦¿$'dw+d09…L"&+9‚%#"&+~`d%€$…8"&±}~…_~‚€Q"L"+$…_#d0LD~!\
“ ”•—–l˜™,š Idd w%$,$…8"& }"~…_~2€Q"LD+$'_d0L R Ÿ9…~¼$F\,›z_L"!& R _#~
d0L 6V&‰gi &$'8D& R L"~'%$…{
R+NQPSR C PT T O‚= Ô ‹ R[PSR ‹ T:PT \ !!d9…M"_L"¦^$'di$'8"&<M"&+ŸLD_$'_d0L R 8%~ $'8D&|%#"&
d0L"&<_#L %L"M
{!&+9…dj_#L %L"M aZF\ N ›z&!&
¤7_¦0"9…&¿E\ÂO¿] &+L"!& y 6&±8F0&
_#L $…8"&!~:&<+d09…L"&+9…~9…&+~…h&!s$'_&!q†”b
E = Ô ‹ R TRR—‹ T TTR
A
= Ô ‹ R R ‹ T T
A
= Ô ‹ R žZR—‹ T aZT
Ü9‚_#Lgi%$'9:_qki€Qd09…g’b
E
E TRR TTR
Ô
R
A = E R T
T
E žZR aZT
A
"! ! #!
$ &! ! !
' ' ' & ( %! ! ! ! ! !
') !* ,+ .. -/ 0 !
!
/
!
/
1
2
43 5 76
8'
' 9;:=<?>A@CB C6 D
' F! E/ ! D
'
/!
HJGI KK !! !!
K ! !
L
L L
! ! M CNO M !! !! NO M NO
P QSU R T P QSV R T P QSW R T
¯Œ&|8F0&|$…8"9…&+&¡&!²ž%$…_#d0LD~L"Mf$'8"9:&!&|DL)6žL"d¨6L"~ N Ô y R L"MX T O.\`¯°&¡!L„~…d0&
$'8"&|#_#LD&F9
~:†~¼$'&!g d€&!²ž$'_#dL"~}dl&<}a†F±
Y "~:~…_#L&!#_gi_L%$…_#d0LS\ec‚8"&±9…&+~…"q$`_#~
TR a! ZT U a! ZR ! TT
!
Ô = M"&+$Z
T
R = ! M"T &+U $Z a! ZT
TR
U
!
ž
Z
!
R
T = M"&+$Z
68"&+9…&
M"&+$œ
Z
=[a! ZR ! TR ‹ ! TR a! ZT U ! TR ! TR U a! ZR ! TT U ! RR a! ZT ‹ ! RR ! TT \
¤d09 $…8"&¿}~:_#~ €QDL"+$…_#d0L\' T 6T &¿¦&+$‰$…8"&¿~'%gi&<g^%$…9…_ kF¢
Z
~}d¨0& y }"D$8"&+9…&^±
]
= N ADC—E0CA0O0_
N ~:_#L"+&`' T _~±d0L"&j_#LÙ$'8"&jL"dM"&8! L"MŠ{!&!9:d¸_LŠ$…8"&id%$'8"&+9
$p6Vd¸L"dM"&+~Os\’›d0_#LD¦$…8"&j~:†~:$…&!g d€
&!²ž$'_#dL"~¦0_&!~
W
Ô
R
T
D‰L"M’~…_gi_u9:†¿€Qd09 ' Z 6_$'8F]‚=
=
=
=
_
!žZR ! TR U ! RR !aZT
! ZT U M"&+! $TR Z
M"R &+$Z
! R U a! ZR
M"&+$Z
3
N ADCADC—E¨O ¦0_&!~‚$'8"&±+dz& i+_#&+Lž$…~
! RR ! TT U ! TR ! TR
R! T U M"&+! $TT Z
R = M"&+$Z
TR U ! RR
!
T = M"&+$Z
> @CB `d$'&$…8%$ƒMD&+$Z &!²ž%#~ W ,N O—68"&!9:& ,N O_~7$'8D&9…&!±d%€ \«›zT &!& B S9 > />>
\ d%$'&€Q"9¼$'8"&+9$'8$_$gj_#¦8ž$L"d$‚}&L"&!+&!~:~'9¼† $'dws$'†j!dgih"z$'& ' %L"M&' Z \5Y_&!L $'8D&
&skzh"9:&!~:~…_#dLf€Qd09 ' R _q$‚_#~‚hd0~:~…_#}D#&$'djgiC6 &±wh”&+9…g¡z$%$…_#d0L d€e$'8D&
LDdzM"&±_L"M"_#+&!~+\
I $&!9:_&wL &skzh"9:&!~…~:_#d0LـQR d09‰$…8"&¿R 9:&FR d€,T $'8D&¡$'9…T _#L"T ¦0& _#L %B S9 >
/
>
>
_#L $…&!9…gj~‚d€e$'8"&±+d09…L"&+9!dd09…M"_L%$…&!~ ! = N ! R C! T O y ! = N ! R C! T OL"M&! Z = N ! ZR C ! TZ Os\
“ ”•—–l˜™,š
Ô
=
!%$
O
&
!"
46
!#
µ N°
¤7_¦0"9:&¿E0b'ƒ9…d0}"&!g Y LDM('ƒ9:d0}"#&+g f
¯ &!& ˆ Os\
(
¤9…d0gФ7_#¦"9…&¿E6V&
!#+"u$'&$'8D&<%9…&F ,N ‚O ~€Qd0#d¨6~!\
(
,N Oƒ= TR v )7v+* = TR v )7vv v ~…_L-,<= TR v ) H¸ v
Y
NW O
d¨6 y +#&F%9…†i$'8"&&!+$…d09…~
)
L"M’¶9…&±¦0_q0&+Lf}a†
=
)
! T U ! RR = N ! TR U ! RRR C! TT U ! TRR OsC
! Z U ! = N ! ZR U ! R C! ZT U ! T Os½
N©Y O
Nµ O
Š=
7kzh"#_!_$…†j$'8"&
9:&F¿_#~$…8ž"~ ¦_0&+Lf}a†
!! !! !! !!
! ! ! !
T R T R
,N O,= TR v ) H¸ v = v ZRR UU RRR TZT UU TTR v
NØ O
N©ˆ O
= TR v N TR U RR O N ZT U TR O U´N TT U TR O N ZR U RR O v ½
d%$'&$'8$‚$'8"&±!9:d0~…~ § h"9:dzMD"+$`}&+$p6V&!&!Lf0&!s$'d09:~ _L $p6djM"_#gj&!LD~…_#dL"~‚_#~¿Lž"g¡}&!9+\
I¯´_$…8 ) L"M^Ùd09:_#&+Lž$…&!M^~_L^¤7_¦0"9:&
E`$'8"&`!9:d0~…~ § h"9:dzMD"+$ ) HjŠ_#~,hd0~…_q$'_q0& L"Mi$…8ž"~
R
,N Oƒ= T N ) H¸7Os\
(
! ! ! !
(> @CB
— Jd0L"~:_#M"&+9|$'8"&j$'9:_uL"¦"u%$…_#d0LÙd€`;Á= ?KAzC W G,H´?BADC—E!GV_#La$'d Y $'9…_#L"¦0&!~
M"9'F6L¶_L
¤7_¦0"9…& W \
46
Ø N°
¤7_#¦0D9…& W b'ƒ9:d0}"#&+g ¯ &!& ˆ Os\ƒc‚8"&±$'9:_uLD¦0"u$'_#dL^d€7;\
*
0-
N 0O<d0gjh"D$…&j$…8"&^gi~:~<gi%$'9:_qk 6_$…8Œ&!&!gj&!La$'~ È = ÑaÑ 3 NxP CDO NQP CDO Ò P Ò y …C Ù=
E0CF½!½F½!C Ø \
] _Lž$!bŒc‚8"&Ù&!~…_&!~:$ 6F† _#~i$'d¶"~:&Ù$…8"& ²žM"9…%$'D9…&„€Qd09:g¡"uŠ}~:&!MÅd0L>$'8"&„%"&„d€
$'8D&
_#La$…&!¦09…L"M y NQP C DO NQP C DO y %$^$…8"&„gj_#M § h”d_#La$'~jd0LÄ$…8"&„$'9:_uL"¦#&~…_M"&!~ y ~:_#L"+&¸$'8"_~^€Qd09:g¡"u
_#~±&.k"%+$±€Qd09±h”d†L"d0gj_u~‰d€M"&!¦9…&!& W \w»‘$
_#~
~…dh”d~…~…_}"#&¡$…d’69:_$…&jMDdl6L°&skzh"_#!_q$±L%†ž$'_F
&skzh"9:&!~:~…_#dL"~7€Qd09I$'8"
& ¼$'&+La$ § €QDL"+$…_#d0L"!~ d0L<&F8|$…9…_#L"¦0& N +€[½ OIL"M|_La$'&!¦9'%$…&
$'8D& h"9…dM""s$'~L%†ž$'_F†0\ c‚8D_#~ y 8"d¨6V&+0&+9 y _~ Jg¡"8´89:M"&!9)6F†0\ Ü}"~:&!9¼0& $'8$ y "~…_L"¦
²žM"9…%$'D9…& y 6&fM"d0#L "B$^LD&!&!M³$'d LDdl6 $…8"&fLq†$…_#!&skzh"9:&!~:~…_#dL"~ y d0L"q%
† $'& )( +* -, .$/, 0 '( 21 $',±68"_#8f9:&
g|"8f&!~…_&!9‚$'dw+d0gihDD$'&%\
*
: C> < :=<
*
\6
B S9 >
µ
X>>
> @C9 > /@ 8>
N }O«d0gjh"D$'&2$'8"& :#"gjh&!M g^~:~«gi%$'9:_qk y 68"_8¿_#~e$'8"&MD_u¦0dLžg^$'9…_ k±6_$…8|$'8"&M"_u¦d0L
&!&!gj&!La$_#L&!8¸9:dl6®}”&+_#L"¦¿$…8"&±~…"g d%€$…8"&±&!#&+gi&+La$'~_#L $'8D&
+d09…9:&!~:h”d0LDM"_#L"¦w9:dl6®d€‡\
N .O 'ƒ9…d¨0&±$'8$ y "~:_#L"¦wL"dM£²"M"9'$'"9:& y $'8D&<%h"h"9…dFkz_#gi%$…&g^%~…~‚gi%$'9:_qki†0d0’¦0&+$_~+$…†
$'8D& ¼#"gjh”&+M ¿gi~…~gi%$'9:_qk”\
] _Lž$!b × R NQP C DO @E
“ ”•—–l˜™,š
N 0OV¯Œ&±~:$'9:$$…dj+d0gjh"D$'&$…8"&
9…&! ƒN xOd€e$'8"&$'9:_uLD¦0#&+~ y e=‡EC W C Y b
E ÖE = EC
,N R O,=
W
W
E ÖE = EC
,N iT O,=
W
W
W E
,N
Z O,= W Ö =‡E0½
c‚8"&+L y 6V&‰+d0gihDD$'&‰
€Q&+6>&!#&+gi&+La$'~d%€e·°b ȸRpR y ȸRT y ȸR Z y L"M È TpT \ `d$'&$'8"%$V$…8"&‰_La$'&!¦9'L"MD~
R %L"M iT \2ÜL $'8D&
R R LDM T R 9…&LDd0L § {+&!9:d¿d0L"†id¨0&+9 R y %L"M T T _#~LDd0L § {+&!9:d¿dl&!9
d$…8"&!98LDM Z R _~‚L"d¨68"&!9:&|L"d0L § {!&+9…djL"M$'8D&!9…&s€Qd09…& ȸR Z =œAD\
R T
R T
R R T
ȸRpR = Ó^Ó R R Ò P Ò w= N R!N T CA0O…O ‹ N R+N ADY C T O:O ‹ N R+N T C T O…O ,N R O
3
R TR ‹ TR TR ‹ A A
T
Ö
Ö ,N R O,= E ,N R O,= E C
= Ö
Y
ˆ
EW
R R ‹ R A ‹ A R
T
Ö T ,N R Oƒ= E ,N R O,= E C
ȸRT = N ~¼†gigj&+$…9…_#
ÂOƒ= È T…R = Ö T T Y Ö
W%µ
EW
R R ‹ R R ‹ A
R R ‹ R R ‹ A
E N ƒN R O ‹ ,N iT O:O2= E ½
T ÖT T ÖT
È TpT = T Ö T Y T Ö T
,N R O ‹
,N jT Oƒ=
Y
ˆ
ˆ
d0La$'_Lž"_#L"¦jL"#d0¦d0"~…q†^¦_0&+~!b
R T R A A T R RT
R
R Õ T R A RR ÕT T
Õ
A T R Õ ÕRR Õ RRR T RR
=
A R RR T R RR T AR
T Õ RT RT Z
N }O¤9:d0g È £= × R È C =¹E0CF½F½F½+C Ø y 6V&
!d0gjh"D$…&b
ȸ RpR = E W ‹ W%E µ ‹ A ‹ A ‹ W%E µ = ˆE ½
E
‰L"#d0¦d0"~…q†b
È TpT = YE È ZpZ = WE È ÕpÕ = YE È ×p× = YW ½
* *
(
(
(
(
(
D
* * Ø
(
(
(
(
c‚8ž"~!b
R A A A A R A A A A
A A Z TR A R A
=
A A A Z AT A A A A Z
N O $ ‰MDM"_#L"¦¿$…8"&±&!#&+gi&+La$'~‚_#L9:dl6®Lž"g¡}&!9 ,¦0_&!~!b
×
È £=®ÓiÓ NxP C DO NQP C DO Ò P Ò ¡=œÓ^Ó NxP C"O Ò P Ò ½
3
3
R
d¨6´"~:&‰$'8D&‰€Qd09…g|"u€Qd09,$'8"&`0d0"gi&`d€e
hž†9…gi_M y L"M)+d0gih"9…&`$'8"& 9…&!~:"$V$'d
68%$†d0 ¦0&s$
68"&+L¸"~:_#L"¦wL"dMS²ž%M"9'%$…"9…&%\
:=< D
* *
ˆ