Scribd
Upload
112 views

Assignment Msg389

This document contains information about an assignment for the course Engineering Computation II at Universiti Sains Malaysia. It lists the names and student details of three students assigned to group 4b, and provides two mathematical problems to solve involving differential equations. The solutions to the problems are presented using methods like the Jacobi iteration, Gauss-Seidel method, and the fourth-order Runge-Kutta method.

Uploaded by

Madi Bidin
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
112 views

Assignment Msg389

This document contains information about an assignment for the course Engineering Computation II at Universiti Sains Malaysia. It lists the names and student details of three students assigned to group 4b, and provides two mathematical problems to solve involving differential equations. The solutions to the problems are presented using methods like the Jacobi iteration, Gauss-Seidel method, and the fourth-order Runge-Kutta method.

Uploaded by

Madi Bidin
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 16

SCHOOL OF MATHEMATICAL SCIENCES

UNIVERSITI SAINS MALAYSIA



MSG 389: ENGINEERING COMPUTATION II
Semester 2, Academic Session 2013/2014

ASSIGNMENT
GROUP: 4b (Adam-Bashfort)

LECTURER:
DR.SARATHA A/P SATHASIVAM


NAME

IC NUMBER
MATRIC
NUMBER
1. MARINA AISHAH BINTI MOHAMAD 921027-11-5534 112556
2.MAZIATUL ATHIRAH LONG HUSSIN

921026-11-5938 112557
3.MARDHIYANA BINTI ABIDIN 921116-12-5216 112555

DATE OF SUBMISSION: 31
st
March 2014







2.
(b) Find an approximate value of for the initial value problem



using the multistep method

)
with . Calculate the starting values using third order Taylor series
method with the same step length,
[40 marks]



SOLUTION 2(b)
Initial value problem,



Using the multistep method,



For starting values, using third order Taylor series method :



We have :


For h=0.1,



Then, use multistep method,




For :



For :



For :














SOLUTION 1(a)

(



)(

) (

) (

) (



i) Jacobi Iteration Method

K=0



K=1



K=2

( )



K=3

( )




K=4

( )


K=5

( )



K=6

( )



OR



k


0 0 0 0
1 1.333333 0.25 0.2
2 1.316667 -0.466667 0.366667
3 1.611111 -0.50 0.276667
4 1.592222 -0.624722 0.722222
5 1.782315 -0.726667 0.806352
6 1.844340 -0.842746 0.847130






ii) Gauss-Seidel Method


K=0



K=1



K=2

( )

K=3

( )

K=4

( )

K=5

( )

K=6

( )


OR

k


0 0 0 0
1 1.333333 -0.416667 0.633333
2 1.683333 -0.75 0.836667
3 1.862222 -0.890278 0.928556
4 1.939611 -0.951945 0.968700
5 1.973548 -0.978949 0.986289
6 1.988413 -0.990779 0.993994


























SOLUTION 3(a)

(



)(

) (

) (

) (

)

Gauss-Seidel Method :




K=0



K=1



K=2



K=3



OR



k


0 1 3 5
1 -2.833333 -1.433333 -1.972727
2 1.495959 -0.904646 -0.849146
3 0.906663 -1.011503 -1.024290
















Additional Questions.
Compute the solution of


with initial values

up to the third zero of . Use the Runge-Kutta(4


th

order) formulas for two first order equations.
Solution :
Convert into two system of first order differential equation
Let



Therefore, we got

] [

] , where [

]



The formula for Runge Kutta 4
th
order method

z
i+1
= z
i
+ (h/6)(k
1
+ 2k
2
+ 2k
3
+ k
4
) where

k
1
= f( t
i
, z
i
)
k
2
= f( t
i
+ (h/2) , z
i
+ (hk
1
)/2 )
k
3
= f( t
i
+ (h/2) , z
i
+ (hk
2
)/2 )
k
4
= f( t
i
+ h , z
i
+ hk
3
)

For ,

] [


] [

] [

]

[


]

[

] [

]

[


]

[

] [

] [

]

[


]

[

]


So,

] +

] [

] [

] [

] ) = [

]


For ,


(__,[

]

[


]

[

(__,[

]+

[[

]]) = (__,[

])

[


]

[

(__,[

]+

]) = (__,[

])

[


]

[

(__,[

]+ [

]) = (__,[

])

[


]

[

]

So,

]+

] [

] [

] [

])

[

]


For ,


(__,[

]

[


]

[


(__,[

]+

]) = (__,[

])

[


]

[

(__,[

]+

[ [

]]) = (__,[

])

[


]

[

(__,[

]+ [[

]]) = (__,[

])

[


]

[

]

So,

]+

] [

] [

]

[

] )

[

You might also like