104

I,
l Prohlt•111~ on l'lnlt<• J>Hf r 1Nw1· Motho«l Im ifolvh•M 'fu,o 1•11t111 J,11u•1u JWJ'
1,

I. Solve !1 11 i l/ I l () with u(n) l 1mJ 11( I) 0. 'l'11kc h :, •

I '.i
, ulvcn tlllll
2. Solve 11"(.r) .r.11(,r) • 0 for 11(r1), .r1
()
1
,
:J :1 0
Y( l) = I nnd y(O) 1 11' (O)
1
3. Solve the BYP 11 • , _ l , y(O) O, 1/( 1) () by f1t1itc <:Jllfo1c11cc 111 ;tll11d 1/ 111) , r 11fl nl
11 1 011 1 411
1
difference method with h
'1

~i" l- P(irJgf'+ c,.(xU~r ~ 9l(7fr)

~rt t ':!tt I - ~f - r Ii~

lJr ., :; ~t +r - 1yr f ~r . . , (h'J- ·

JOlu.H~ :

I • u''
-.i 4 IJ + t = D, ~ [o) ::-1 ~ Ct) = o. . •. Wl- {tw;Je t-D rViJ o i 9t,< oJ((/1,,(),
1-;: ( Y1t) 1Dtl r "' I er I 2- J £' .::, ~ .
8t)t :
--- ~i '' :: ~Y+1 -#,yr+ yr-1/ht. fut f &- 1 '.

fl,, (12.. - 'Jl/1 ·I ijo) 1' ~ I -I- I s;, Q ,

'dt ,, +y,· t-1 r;,f)
Ib (y 2- - .;1y, -I I) +~ 1
4 f c- O,,

~ft r - J'cJT tYf-1 lb Vi. -- 8J2gt + /(, 1y1 -t I ~ rt)

1 ~? f I 'C,'D •
Jtl- /b~..1. - 81Y 1 :: -1::t- _.:r{!) ,
""~ (¼) k <t-_ 1//1,.
Pitt r ~ [;):
l&('Jf+r - J~rt~H) +~t~ I i:-0.
/b ( 4 5 - t¼- Y,z.. ..f ~I ) f ':f ..2. 'f f .r O .
Jltl.h :
(btJ 3- B~<JL +- !6 tJt -t qg :: - t
o
. '
Ytt tb 9t - 8lV1. + 11, f/ 3 ~ - f - >@)

Y(C>) t I
~llJ ~u•
 Boundary Value Problems ,n Ordinnry and Pnrtw l D1/Jc rc nt1nl /:.,qu.Jtion ) 195

pu.t: p C 3 Pu.i- t 1 ===

~':. .1\(4.1-JY,-+ ~o)

11:,(!J~ r~g3 +y.2) ty ~ t/ ~o
~\ ~~~~-3~\
- 3.;)~3 -tfb';f 1-tl:33.-1 t -=O
-9.~t+k'.b--~O ~(D
(b \:f !2.- ~ BC \f ~-.:.. - l ~ © PU!t- f ~ ..2
cri<. s clvfltg ~ 3 - e9 uttH crw Y.2.. = A (4 ~ -~y.i-+ y,)
~ tKe.. co.J,c.~ "'.2-- ~ ,4~~ -S.tj !l.- ~ l-\~t
t.02~ A~~--t_\2._-\µ.~t:O ~
tj,~ o.. C/)50 put T == ~. _
1j2- :::- 0, =TDIJ., ~3= .A(YA-.J~3+ld2.)
\J; =- f).3Cf8'} ~; ~A~~ •~~3-+ l-\~1...

, y 11 .c <J ., ~ (v J !!-o y (2 ) r 3 ·b 2~ ~~ ~ A.(3-t>~ ➔) ~-a~3 <\- k'!2..

1

h = (~) J.\'!i ~q-\A 3 .: - f~ ·l5D~ ~©

Sd\)1W\ (D ~ @) ~ ®
we.. caet .
0 Ya 3/i :J 9' :'.: 0. ~ ~6 3

1J,.'v~ c· ljo
"'X.1

~I
")(. 2..

ffj!_
?C 3

y~
~~
g, 1f~"'1
J ~ ,__ ~ t'" f%-4 3
~~ ~ J. t 138; 8-
J ~~)
H{DJ=- D. = 3" t,.;i 1 J)if 1
'{x) -xg[x);;D
~ {, J = / <er J [o) + ~t ( e J =:: J
k= --k
~)I~ yl

l,/c,
v., - v•+'
- LJ., -).llj-1
'iJ
tlp_, /hi..
J • .. 3/J.
~

)\
L=- j_ n.:z -_ _lJ.t -Y3
0
1/3
rt 2.

Yr~ tJ-, . ~ (x?}; ~ r-1 r - JP;;!

Hr == A{~ r+• - a-~r-+ ~ r-,) :J '{o) =
('
teh
/
~'1.
196 Work Book on Numerical Methods

9
1akl~ / )·

: , )j{oJ -t ~
1
(0) ~,
'/ ('Ji -Jy3t\fr J- iY'- "D.
~c LJ ~J ., I 2 ~lJ 1 B~ / - 3 ~ ,., I ~ tV qr:f{,/ -Jl,U -o
..f 1- ::: I
1 J-:f y3--!:H-1.~2,-fc,c o 1 JJ.. -
;;../ 3
;1.-:f [1 )- '361 .. -r 2. 1'~1 = 0
~- 1 = f (-2 -t:i~o t-B~,)
- .1.'f~ I + 5 by z " ,¥.1 ----'/"@
~ (,) ::: I V3 -: I

<Jr '1 = Yr+t -;.2.~ r -t ~f~r /h

2
50111 I Wj (f)., ~~~
-~y we(fif
~ YP+ I - d ~ t ~ ~ r-1 I (t) 2 0

r
~v = -0· qg,g-O
ta.,k:.~t,1.J T ~0.
~I~ .-().3iq 3
q (Yc·t'-~otY-1)-'XoYD =0, y2-;; O-JJr:53
C/~r - r&~o +°I (-& t ~~v ·t 3 ~, )-0 ~1:
~ i, ~ ,1 .,_ 5~, + A':J ==-,

Cf~ I .,- f 3~ 0 - b -t b 40 t Cf y, s:- 0 y,~, ~ t ~ (,) ::. 0 .

-t~yv 1 r8 Lj r -; b ; - t2!j O +iy, !;: f h =- ~ .

~co
1a)<1 veg 1 :;; t
'7/ {ft -~g I t{Jv} ~ K,, y, :.: 0
1 [~ 2- - q~ f - ~ 0 ) - ¾LY I ) ~ l)

Of ~o f 1Y2 - 18 YI -- JfJ_ = V
3

qJo +t~ .i - f-%!1-1 ~6

3
J;r yo ·t-2- :::r~ ~ -- 0 5 v, : :- o .
/::f ljv - 5 6 if 1 +21Y;, ::[) ~(!J
198
Work Book on Numerical M etlwds

II - Problems on Finite Difference Method for Solving Two - Dimensional Laplace•

Equation (Elliptic) s

l . Solve 'v' 2 u = o numerically, using S-point formula and Liebmann 's iteration, for the following Illes
with uniform spacing and with boundary conditions as shown below: h

0 It> ;o 30
U1 U,1- - M
,'¼ UJt.,

bO
60 bO 6ti

2. Solve 'v' 2 u = 0 numerically, using 5-point formula and Liebmann 's iteration, for the following Illesh
with uniform spacing and with boundary conditions as shown below:
f) 11'.)0 zoo to'D 0

U, \lL. u,~

40-0 Cl11 ti.I; I¼,

- )too

2.0 0 ll:J b5l ltt4

D
0 ro o zvv (0 l)

3. Solve 'v' 2 u = 0, for the following square mesh with boundary conditions as shown below:

4. Using 5-point formula and Liebmann's iteration method solve 'v' 2 u = 0 over a square region of side 4
satisfying the boundary conditions

i) u(O,y) = 0,0 :Sy :S 4

ii) u(4,y)=12+y,0:Sy:S4
iii) u(x, 0) = 3x , 0 :S x :S 4
iv) u(x, 4) = x 2 ,.0 '.S X '.S 4
 ~----------=====-------- lfo1111dtll) V11 /11r 1'1111>11•11111; I , I l

UJ ft! ,
1il/11111y I

I \
1111 / / '11111111 11111
f lrt/11111 I '/lliili1111 11 , 9(1

), ' \I '\ ~O \ t1 11u\ i ll ~ ''-' "\

,1i\\~:

I ·H'
~<'Ill
t·it(I ll<.o.
11
litt \tC\\ ~(tlU Q1 °' ~ 1 \, t> '\ \\\) , \ ~, t o l
I' ti } ' (I ~ I u~ . t" ti •1~
111 1 I)
101
LI 10,. lt\ (0 1 u~ v , we COM U?, tt) :: \ lAtl'i bO -l lJ}
0
'...\ ~klll)l
I1, • J '
,~rff1t'I 1.lle i /a,uclm1 cl " polt1t~
• ttrUL~ l\ WO'i 5· pot//. t fonWJ-litJ ,
:: 1 LA tH bt t \5 -1 D 1
fr
: :rn.1~ -
ftf l~ h,t,l<l!
U>i ,~, c 0 ~~~
LlJ\ Cl] ;:. -t, [
60 i l O i ll '-tb) .\ U~lo']
sy DfPF

u)o' -= .lJ.., flJ..\ t> .,.~ oD,/J

o1
~ 1 ('rSc-t b l':I +li -~s t ..28 ,'1 '3J

. . . '?i (be)
l
u~ :: D
\ I>
-= s~. 3~r6 .
Jteno.H tnll a, :
u/01
B9 SFP f
,c 1'5
u}'l "'* f:2.0 t I~ t U1ult ll~ ~l J
u/0 ' ~ 1;:, LJo t J-t l) 1 u: 0
\ Llk\c1
t ~ [~ ~II> I l&.'113 t ,i-<1,;,]
u.t'l.);... \9, ~ ~ 6.
~+i [10t)rO tlfJ-tDJ
lht,) e \ \!D t RM ll1 {t l llt/''} ~
=- :f!>

Ll;i
lo)
~l8,-:J6
1-l ·
-: ~ r~ tH J.{ t> <\ 19 • S 1f5 f 3 • 9 1 ~ ~1
~ J't, 68';r5.
8\j SFPF
lt3 lo) 1F~ H.e/0\
';:. D -~ (,t ij~lo~
~ ~ t:,.) :: 1,, [.4 a t iv 4 U I t\\ u. ~t\) J
~ 1 [1-1c+-bo ➔ rct-WJT :nwsJ
~ ~' fJ-lt)-t bD ➔ 16 ·tO J -:;: sq, b9f6.
~· lff>/A 4M M -.. AL
Im vu +Go+ ,s-~~ +2Q.'l'=>}
~l) ~ JB, -=1s ~ 81 -3 .f6 .
1teJtaflet< f) ,' Tte.n ti~~tJV.. ~ :
11 (1)~ '
' A t JOtlO + u.,_to)+ ll3(0JJ
U t~) -- l r.
' -- Ji t20flo ~ .tJ ,G&:/5 1 S~•h&=r,]

-;: 1 [;vt 10 ·t /8 , ,;f ~ -t ttPlq

~~J.l-81-\3":f- ,

~I 11) ~ )1 •3 16
 a ,,, e. , l 7
.nci0c\'-~ c, :
j, \'l'~Jtn"' f-1 ~f~H ~•1•\'Pf'-il
1 c, 0 ~ \ ~ 1· ~2 ·l-1t ISJ t- tr
l

'- "1 !,.~ \f \~.

u.,l(?) C:
2 · /1.? le]
(t~ ( ,, ~ t f,to~ l:,r 1 ,, . i:;;...,~ . . :\'J, n0 :::= ,_ b . .2 m'-' .
I. t.h ) ~ t [20 ..,. I{ 0 -t 16 . .2.I 09 ., ,4 ~ .1.1 t),]
~ y), £,g ~i:, u,2. -~

~ S3 · IOfUt
Ct.t:/Tt.J r. J [ D0-16'0 -t .J't, t8 i'=> t .3"1 ,b8 qF;]

c. A 4, 8A 3 ;r
]}01,o._~ eM. A :
J. 8htct tcw_ ho W,{~mltf VtLl~ ~
u/A l ~ ! po .,. , ➔ o 29. h~-l15 + 3 q. bs "f~ J vw1tl°l) ot-Lto))JJ t ~ wet,uc ,
U ' ;: :. U b :: Lt ':f t::. l.AC, .
=~A.8'13-;f.
gfVlh. fCa. bo W1ciaJVJ v(J)»lh a.JUI
ll:2. CA} == °i [Zt> -t JH>-t ~. ~1-,3 :f + AR• t$ k3:{1 vm.Hccul)J J ~ ru_ e.txl c .
.. l{_b ::. {). g.
= 3J, J\ 2..\8 · t,t,_!:.U.~ -

U"b (II).::. t [1\o -t tot ~t-l•%k3::t ~J\1-1-8n3==tJ

lJ.fj = J_
r.,
,t t!V 0
.,.~o01 IID07.ltoo]

=- A a, l-\2.18 . tt5 _-:;::. soo

uI\ tH) ~ 1 [so + bo + i <1, b8 ; f> -i 3-, , b8=njJ u, :: *r )ft>O ttOO t D t lu,]
-::. A~,8~)f . '""i [1'00 tZOD ' t- 3 ooJ
~affr-M JS :
u., :: 2 16 . ~ J)FP f
llt(fi)!:-¾f.2.0110-t 3.2,~2181~2. ~2.l&J t, U-( "" U°l = J'25,
" lL1 !:: IA.~ e '

~ e} 6. d IV °I .

u.,_ (5 J ~ 1-
j.,
[ :iv .f Jto rl- .2>-1·~ H3=7 -t AH•~ Jf3-0

==-% 3~• )l2. l8 ·

u3 ('3l ~ l [J-tD t6o-l.2J..t,%'13q4lf1-t,S'1~}

J..f

-:: :. ll_2., lf219.

ll~ (-,J ~ 1II Urrso t 6-0 .f 3:i. I A2 ra t J\1' tt218 l

 Ctl\4 c
JrnifiVI ·
tf(U1~) ' O . (.) / l J ~ /1

Cl1,'t.l 15 u (l.~,l\l tt12"''-f, DLi?:J ✓ ~

u. r')(,o) :-- ..; x: , o ✓,.. /"

u(Y. ,k) ~ .,1. 10 ~ x r A

(O,O) (1,0) Jl Uk :: J.,.. [o+ e, t J· 5 + 3)

O 3
~ ~.%"cf~.
U,; .,, ~ ( ~ + I.It .i. b ,t- oJ
Jte,,ll.\t~ \:
U.5 : 6 I (1 +o t 5 , b ~ B 4 d•i i ~ ) ~
tt, ~ 1
U-a -:: 1 (A+ll> + lk+b)
::: J· 3:f'3.
(D~PF) Ll.,.•i(.l \+a,5+1 n+Gl: 5 . bd5 .
~ 10 .
U '1 = { ( 6 t I~ +11 H )
Uo ={ (q + I '5 t'6· b~5 ,i. Cf, 8:f'.\) =- q,~ :f~ .

~ q.~ ui\ ~ 1 ,o ~ at '6 t (, + 3) = 4'· 8 1 5 .

lti =.. * (~t 4 t 6 ,tO)

lJ. 1'> = *( O, b J 5 1- J . t,<; ,t <j , 8:f 5 t 1, , I A '3 ]

{~Fpr) = l, ~ ~ ~~ t 6'6 ·
1 (tl\l-
~
U.1 =- *( ~+6-t J,5 tlD}
Uc,=

u:~ ~i. (o+ 5-+ ~-&~5+b•1~~~

\0 +t +~-~) :: q .g-;f~ ·

e S,
::: ~. 615.
u8 ~ i (+ b b t- j +Cf '~ J ::: 6 It r; '
I

ll tii c ~ (Is .f- 'f + 'f · 8 -tB -t 6 .,i 6)

J.\

~ "I• S"f5 . = 9. 6,,

{;. ;r.
{Yf pf)
~ *{ o + b +-t to)
lfe)l cul rN a. :
u,: ~ (H-t) 4- ti• b!~4 a.&~~)= J. 3~5 -
~J U1 -:.1 (1\+ a•r~~ +9,~1~-+"•'d6) c5o513
lls -;; 1(l, ·t ~ 6 +3 t b)
0
_lti~i(9t l~+ fi ,bA5+'i ,8:f.S)= 9 .~:p5
,;l c. JJ5 UJ\: i (o + ~-i~5 +6,l~~t 3 \ c J.~:r~ .
 lt1' { (~ -l-A.~4 ~ ~ l~-\ '1 , i';f~-4 6, IR~)
~ 6,1.~"' ·
~6 ~ i (1)..\..\ 9 ,8!f6-4 G, r~,f 1,5) -=~-8tt'5 .
( l1 - *( 0--\ 3 -¼ Q , 8 l'o -\ e,. p~ F5 ) ~ 3 ·
la ~ ~ ( 6 + 6,la6 ➔ ~ <1- C/ • 9 J:::. b• f 6 b ·
I

+ ~ · & 15-f t'- f,;z 5 J ::

1

6Cf """ j; ( r 3 t °1 Cf· 5 ·

 - nnd JJ11Hfol Difforcntfo/ 1'
BoundMy ½1/uc Problems i11 Or' II no,y .
:l/UIJ.LIOnH 205

for Solvi ng p oisson's Equation (Elll ti

JJl - probleJllS on Finite Difference Method
•
~ ·
~
With h I, solve u,, +u,. -lO(x +v' + 10) o< x < 3 0
2
l a and L'1ebmann 's iteration ' -
· method.
3
and u - 0 th
-
P c)
' ~ y -< - on · e boundary
1. using 5-pt 1ormu
lVldtng the sq .
z. Solve v' 2u = 8x 2y2 for square mesh gjven u = 0 on the boundari es, d' ..
uare mto 16
sub-squares of length 1 unit.

lJ.;d)

U.1 ;; A "' [1) ~)

Lb =-: B =-- ( 21 2,)
U; ~ C~ (I; I)
u~ =p :;- {i) , ) ~ ~ L-o,ol
(d,u.Q, to 1~tttVU~Wu~ ) b Mr»i take liA ~o .
PD,i~• s equa.tt~ f~ give.u. Y:

i.l?,l = tr Llf-1 ,Y t U.1tl, 'j -l-\!~, f-1 +u~' rt ,1- r rt ( 11'.f, YI l]

At pof~~ A => ( L~J. /' yFf'F

u, ~ Ji fD+1HOt%l + !£.
t,
1"+1-""t10J r
-;. {to)+': (P5) c ~1~ ~,

A+- pol i~ .B ~>- [ :21 iJ

112. :::-l fo+o+u.~1U~ -t ~ (2-l-+~+\tl ~ J_ c~~-~t~r+,~\-t ~~
A ~ 4

U:2-~69 . 8 ,
At-- p1H vtl- e 9 ( t , 1)

u6 " *r uI t- U4 H -t-0l + l! (t ~ ltro ) " 1(~-pi +'3'P)) H D

u~:::. A~,g
u~ ~ u, l <Wt tfJ Jr;,<,u~eb<l tlfy ) .
u~ ~ 3:r-, 5
 206 Work B ook on Numerirnl M ctho<is

lt2 Ct I ~ ::ff ,°\

l\3 ° ~ i [o +o -+ u, ;-u 1+ ~ Ct+,~

1
4 10 '

u.~tt)t:: 6.2. "l ,

U.N tr I ~ 0 5. ':+ •

J. tJ~ r: <t Z:11'-- ll.'•= V

d1vtdlt/4f tf<JZ ,3 ff ttcVW !'t1.lfJ Ib J uh - iq_ua.x~
Df i.al9M t tJ.Lt.t f: ,
y
( ~ ~x..
0 f) -r;< '
D 1)
,,
,I
0
J tfM--VUlfflJj ~J t- LU(.
ll ! ✓

... ' 0 ti/-1, l),,~ U,b V

,,
4t '!::.UJ -c t1 ::r- ~ Ucr {
✓
t>
0
ftt Ll1- ,u
U) ,;:_ ll. 1-1. !:_ tl b t U.g- '
V
t) 0 V D

let- tM ta_~ [u~ l::- 0 l

u, ._ i [ otM o -1-- u-sl - \"- ({(11 , ~1 l l =
\Y-
I J - t)
A
lt.t t ~ - ~]
'"1
*fo 1 - :. [2C-1J""t-n" J
"

...
6 • u, ~ u~ !: tl.-::r c:. u't <;: - ~

U~ -: . U~ ~ ll ~ ~ Ua, ~ - /
 ~
1
lY \\n,hlt"ml'l(<Hl Ont- l>l111rni.l1111nl 11,·ttl 1110111 J 11untlou (1'1u nholl,) U!J l~s.plldt Pti<•LI

.
(l\r-u1lr1 ,-.., hm tit M1, 1l111tl)
. ~
h>d

'
I !'i1lht• /Ju 1//J", ,\Ith tht' l11ntl1t1on" u{ll /) 0 11( I, t) , 11(, , 0) J'( 11 J) tuk1no 1,
lit _Ii), • t)

\"mrh,, ln~ l\c:ndcr Schmidt rcc11rrcnc1• rq111\110n Continue 1hc aoluuon thrnugh I() rime 11tcp'I.

" 1 ,oil the 1-olut1on to 11 1 • "" subJt:Ct to IL(T, 0) si1urJ" , 0 ,., .r. ~ I, 11(0, t) u(l ,I) () u,,0 0
Schm1<.lt method "

1 Solve by Bender-Schmidt fonnu lo upto / =5 for the equation Ut = 11u subjecl lo i, (O, I) u(r,, t) • ,,
2
nnd 11(.r, 0) - x 2(25 - x ), t.aking h = l

~~ut-\m.u ,
t-J
~ ~(o, t-) =- D !: U.(~,
' •
d t- L ..L~
'l. r}X t. U ('X , D) t 'X.. (/.\ -xJ
~ :: 1 -=-) 1D ttwe. ste,ps -

~.\:: ,'\
a'~ ~ c. 2. l~
Jt ;r1..i.
cow.p wlnt '}lVW =-> we ~ c-z--= ; .

1 t = c21c
h1. -
';\_ :: .L
R)91 B,s &"
(.L) k
J__ =- ~
J. Cr) :i

k t
I
-')(. -
.2.
2.-
I
I§]
0
=
0 8
0
0

e,
J "' If