0% found this document useful (0 votes)
57 views13 pages

Solutions To Selected Problems of Assignment 5

1) The given document describes solving a second-order boundary value problem (BVP) as a system of first-order differential equations. 2) It presents the steps to solve the BVP using explicit Euler's method, including defining the variables, calculating values at increasing x values, and comparing the final solution to the given boundary conditions. 3) The solution depends on an initial assumption for the derivative value, which is adjusted to improve the accuracy of the final solution compared to the true boundary values.

Uploaded by

Shivani Kumari
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)
57 views13 pages

Solutions To Selected Problems of Assignment 5

1) The given document describes solving a second-order boundary value problem (BVP) as a system of first-order differential equations. 2) It presents the steps to solve the BVP using explicit Euler's method, including defining the variables, calculating values at increasing x values, and comparing the final solution to the given boundary conditions. 3) The solution depends on an initial assumption for the derivative value, which is adjusted to improve the accuracy of the final solution compared to the true boundary values.

Uploaded by

Shivani Kumari
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/ 13

2.

Assiqument 5
d -+0*5 y +Z =0 (o) 4, y'(o) =0
Ste size heo. , Loe need to ehate y at
OD E -IVP as both he
gven frsklam rotlem
Comditiony ane specifed at aa single value of xfat ODE 0)
But, aete metht that oe have "stuoied fe sinynve roblms
ae meut fnfo sebing firnstorden oDEs ¢ the fon m def(x)
os, sinee the given "oDEs is becond -oPder ODE , we need to
a

dnwentitinto a system of tuo firt-order ODES We do that


by toting dofiing tuo vaniales y y2,sueh thot y
dz
With thaet assnfto,he gin
. WiK
2 nd mdar becomes
t0-5y, +7y =00
' h e system f DEs that we need to sdve may be given
d
at t=9y 4 -2a
at <=0,y, =0
Say thet fi= & f Sowe =-y-0-5y
Cmfaningwi't shadadthefon my-f(»),
above ODES any method that can
te ill bove
sbacified , we can chese
DE-|VP. inee no shecifie metid is
exbit Extor metod x=0.]
Caleulation of 4, &.y, at
Yo= 4 u the initial vale of
o0

Then, y, =yo thyse Atox0 =4


Jo-Bo0-5y..) = 0to(-7X4 - 0-5 xo) =- 2-8
at O2.
Caeulation of ky
=4
,=-2-8
3-72
, t hy,, =4t o.|X (-2"8) =4-0-28
=
Yn
t-y-0-5y.,)= -4-8 +o-1-7x4 -0o5x (-29)
J2 , =-5.46

nefens too y,y at x=02, so


we
Mot that y,, tt
hawe got
We ottain y =3-72 at x=02 (A)
mettol instead
ofRen method of Bplict
Note that if we consider any value f
l:ke RK4 n Heun' metkod, the
Eler methot, I;ke
Howeva, tRis differne
Jifeent.
tRe final ayer ill be indeel smal stet size
wiil neluce shen we teke very
(ia-1e)
B - 0|5T-0 j T(0) =240,T(0) -\50
T=ce
m
’ dT Cmem ’dT nemi
(a) Let
-05o m =0
’mo-15 =0 =t03873
’ m=±o 5
’ m= 045

he e*+,e"a -03873 x
T= 0-3873 x
+e
(2)
T=ce
we
Can oite
the brunday conditoms -03873x 0
Gte=240
03873X0
Ce3 +ee = 240, ie,
to-2c4
=\50,iey 48-09, =|50
-0:3873 X|O
e03873xIO
te,e
and,
We necl to solee tfe
C+2=240
48-09c,+o02 c,= |S0 -31)
Q"02.
240 X0.02. |50 x =3.0
0 02. - 4809
4809 0.02
240
48 09 150 |X 150 -24 0 x 480g = 236. 98
0 02 48 Og
|48-09 o-o2

The anaytieal sution


0 3873x +236.98%-0:3873x
T =3:02 e
ODE BVP

) The given frobtom [refen to equating la - le] ia than veluc

problem, as the Conditions areNZ=0rorided at mere Since itu an


f (Conditions pgvickd at and' x= lo) .
s e this
Can be wwd to
froblem, shooting method
ODE-BVP þroblem, shroting
freblem .
Naw, we wmte Ens (la-le) here once grin
T(o) =24o, T0o) =(50
JT -o-|5T=0;
We define T, =T &T2=dI
d
Uhe above 2nd- order oDE
be exfnemed an

ay
tuo finst- onda oDEs:
a ystem of
at <=o,T=240
dT at =l0, T= (50 -4)
=0|57
methal inswes tran foming the bondary vkie
Sineaing we need to asme
ooblemninto intil vae froblem (IVP)
vale of , at =0, ho that ean (4) Can he selvel an an

IVP
S t i e i A pertinend quetion is we

the value of T at =0.inee T decreanes frrom 240 to 50 as


incneaes foom O to O, wwe ean say tht dTdr is hegative
refenhed an
as
be taton
m, katen on
negative. 4)
Let us take 0 at R e |50- 240
0 -9,50
Can tade Jometting elose to that, to sot witk
S deroted asTyo
subseyoty
=T
(at i0, T, =240)
=6: 15T (a a=0, T; =-|0)
denated as T2o
subsugyanty
Let ws ove teis AIVP Bpo blm by exbicit Eulen's method.
Netake Ste size = 2·5
as given in tee question
;
Galculation of &T, at i=2:5
Tio = 240, Tzo = -|0
. By erplicit Euler's meto&
T, =To +h Tao= 240 +2-5x-10) = 215
2 =T2o tho15T,)=-l0+2·5 x o·15X (240)=80
Calelathon of T & T, at =5
72=lu thT = 215 + 2·5 X80 = 415
T, =T, th (o-15 Ti) -=80+ 2·5 X0.\I5x215 = l60.625
Clalati on of T & T at x=75
T,=
13 t T = A15 +2·5 X160.625 =86.5625
(23 =Tb, th(o"45 T2)=|GO-625 +2-5 X0-|5 X4I5 =3ls-25
Caleulatien of T, & T, at =10
2 refenre das
75) P later
4=3 h7,, = 8\656 2.55 4+2-5 X3I6.25 =1Go7 18

T4 =T23 th(o-15T)=316-25 + 2-5 Xo- 15 x8i.5625


=622- 44
Calculatel
asumption
Zhl 4 is theavalue o T, at x=|0,baned on the ammtin
ef =-|0 at x=0 , and stef site (h)= 2:5. On the
vfhen han,the actual value
(50. Cleny , teone a
Auge diffone .So, foo scomsidn
ts
$2o mating tke second ansumption (i.e, m) oe
m= =-40 hi it again a gues detila om
this to be diucunad later.
soved
he TVe Brotlem to be
at na0, T =24o

= D.(5T, at x=0, T2 =-40

let
Caleulation of T & T at z=2-5
ho =240, To -40

Ti =To t ATzo = 24O +2-5 xxc40)


ho) =140
T =t2o +h-(0.15Ti) =-40 t 2-5 Xo.5 X240 = 50
Caldation ef & 1) at =5

Ti2 +h2 = \40 t2-5 X50 = 265


Te2 = T2 th (015Ti)= 50+2-5 Xo-(5 xl40= 102-5
Caleulation 1; & T, at 7-5
Tia=T2t h122= 245 + 2-5 X\02.5 = 52|·25
l23 T22th. (o (5)=f02-5 +-2.sx0-15 X265 =201 B/5
Culeulathon
refermed
= 52): 25 +2-5 X 20|·8 75
=(l025 9375
T24, = T23 th (0-)5Ta) = 20| · 875 +2-5 Xo-)E X52|- 25
397.34
ls, swe get anthen Tated a neu asfian
of Blzso (ie, m) -40
So, we Khawe the

Assumed vale of ,m,=-\o m,=-40 m, ?to ke


fund)
|607 1875 025.937s 150 (giver
| =) () browndany
Çonihiom
fuesin)
From the abrve tabe, m, may be calelate g
via linean inteshelation technique , as not intabeltim
(stricty speting, e we ae doing inen erte apolatinin ttin
deln
,asit
on chice

150 -(025-93 75 l025-9375 - l6 07. |8Z5 Tdealy,


be chen
’ m3 = m, t (15 0-1025-9375) (mm) Such that
(lo25-93 75-1607. (875) intetbadatan
’ m = -40 +((50 - |o 25.9375) -40 +|0) Hruevn,
these
(1025-9375 -(607. |875) do hot
’ y -85.2
Ly Caleation of my constittes the Anst Setion,You
provided
;since the
thewe given qustion aoks penfom a
wtoiteationThus
þenfonm 3iterctions, wiu bo calculte m, ng of iteiations)
Che T ve 6reflem to be

at =0, T= 240
=0-|5T, at r=0, I--85- 2(ie,
in the

Let soue this ootem


Calculation of T & T, at x=2-5
Tio =240, 20 =-85-2
Ih=o+hTzo =24042-5 x-85-2) = 27
T4 =T2o th (o-15Tio) =(-85-2)+2-5x0-} 5X240 =4-8
Caleulation ef T &T at =5
Th=T th = 27+2-5X4·8 =39
T = th. (o-5T) =48 + 2-5 XO·15X27 = |4:925

a=T2thTig =39+2-5 X|4.925 =76-3


T23 =T22 th(0-15T,) = J4-925+2-5 X0-|5 x39=29.55

y = 3 thT23 = 76-3 +2" 5X29-55 =150.I75


4 =l23 th-(0|s T)=2955 +2·5x o. 15 ×76-3=58-16

So, we get |S0.|75, banel 4he aswnpim


thus, we Rawe the folruig tatle
Assumed le of Eenm-1o m=-4omy85-2|my=?
|1co7.1875 (025-9325 I50-|75 150
8)
Note that athouh gueses of m, m e almey
Nou, We wle be B my in onderto andmygi
caleDto
my vá the te chnique of linan itenpotatn Ths,
may wite:

150-)50. |75
|50|75- l025-9375
’my= M, t (i50 -|50 .|75) (ng- m)
(15o-(75 - l025.9375)
my=-85.2 + ((50- 150-175) (-85-2 440)
(150- 175-l025.9375)
)ny =-85. 209

Calculation of m, contitut the second iteation


he IVP prroten to be soved
at =0 T=240

=015T, at =0, T =-85-209

sove this Brtlem


Caloulatim of T T at =25
Tio =240, Tao =-85-209
h =Tio thTao = 240 t2-5 x85- 209) =26977
I =2o th. (0-45Tie) =85:203) +2:5 Xo-(5 X240
= 479)

1, = T thTz =26977+ 25x479| = 38-954


22=T2 th(o-15)T, =479)+2-5xo-15 X26-577=l4.907
Caiculation efof T,& 1,T, at 1=75
y =7+hT = 38-954 +2·5 x J4:907 = 76· 221
23= 2 th (o)ST,)=14:907 +2 5x0-15 X38954
= 29.5)5

Ty= a+hg 76-222 +2-5 X29-515 =|50. o095


T2, =Tza h{o1ST:) =23-5I5+2-5x0-}5 X76-222.
= 58-D98 2.

So, we fad and get i . ie.


Te= |50.0095 with thut
So Since actual value of TL = |50,
may sy that r reults have cowergd
ith o3
cuneey (cack fm yonef) , and
may stef ke
n calealationyNefe fhat the
do antter
ioton,
leh
bu as
wtion asks w to do 3 iteratony tut we are
wndestandyate
there il}) to get highly accnate neydls with juf
hardy be itentiony.Tlis is becaue ow choice of m
any hang
tn the viuy
m ongh andom ,wm made inttively and
takenas
nConidenin the bruday comdtiony That elly
hehed w in eting coengane wit jwt 2ijenat
We derstond the poen af sheeting method,
the acuaey of nesult that it can offes
Solution
with
f the sien frootlen by ghoing metid
stef size 2.5 a TIeehe =26.97,T=38.954,
Tl=75
e7s = 76:222
(Ans).
(C) es, fimite can be wl foo tkis Brotlemy
difference
finite diffenca metho Can sowe ODE - BvP

We have been ahel to tahe a stef size ot 2,


there ill
the domain of x (os*<ID) Thur
be 4 intien ndes, as defietel belos

i=h i5

hteron ndey

Nos, the given


T 0. 15 T=0
viA Cudral differee
Vpm dacneising the eauatin
fonmulation, we get : -0-15T; =0 6)
A) T refesto the
tee steb size, and
shere ) s
ia gn
node y=2, abeady
femßerate at e ih writen fom all thefou
eqmation is
The disenetizel
intenion nodes (i=l, 2,3,4) -
T, - 27, flo -o-\5T, =0 6)
2

-2T, +T -0-15 T; =0 -)
1=3 , -21; t 2 -0. 15 Ta 0 (e)
i=4: Ts - 2T, t T3 - 0·15 T, =0 )
Thus, e Rase To, T, Ta, Ts, Ti, s to obtain from
f esuationsEs(c) -(). Hasewen T, &Ts
nefu to ga T120 Tlaao retectivey, and
we adready
kns hem To =240, T =\50 .

voe need to selve for fown wnknowns :


Ta), wing fo 13,
eqnations Eas ()-).
Ea: (o)may be
,-27, +240 - 0·\5TI =0
(2)
’-2T, +240 - o.6T, =0
’ 2-6T -1 24 o -Go)
t fonn ey ),
0-15T, =0
-T +26T, -3 =0
-6)
-Tat6T, -T -0 -(2)
Lastty
(2
Tt- 2T4tl o-15T, =0

’ \50 -2T4+ l5. -0:|5T, =0


22
-Tg +2-6T =|5o 3)
.".Ean (o) -0) the set
f equatims
hat neels to be stwed to ottain Ti, T, Ts,T,
Lat wmite them beles fom Covenience, so
tat we can fo m the system of eqwtims
in matain nottm, AX = B
2-6T-1, =240 -(0)
-T t26T; -; =0 - -6)
-T +2:6T;-TË =0
l3 t26T, =l50 --(2)
Egotios 0) -(13) can be nefruented u
240
-| 26 -! O
-| 2.6 -
T
Lo - 2"6
150
he ahve

ohere
qpation is of the fom
foom AX=B,
26 - |
A=
- | 2.G-| 0
-| 2-6 --
-| 2-6

X= T
13

240
B
LISD
The matix ean AX=B needsto be soed to qet xX.
It üw easy to otsowe that A à a a tridiagonal
matix . Ts,tehis froem can be sowed coweniety
by Thomas adynitn, alyo knn as
Matrix Aly orfthm

You might also like