Scribd
Upload
214 views

UTP CM Matlab Assignment and Test Questions - May 2015 Sem

This document contains 12 questions regarding computational methods and MATLAB assignments. The questions involve implementing numerical methods like bisection, false position, Newton-Raphson, and secant methods to solve nonlinear equations. Other questions involve solving systems of linear equations using Gaussian elimination, LU decomposition, Gauss-Seidel iteration and Thomas algorithm. Additional questions involve interpolation, numerical integration techniques like trapezoidal rule, Simpson's 1/3 rule, and solving ordinary differential equations using Euler's method and 4th order Runge-Kutta method. M-files are to be developed to implement these numerical techniques.

Uploaded by

Nurul Aishah
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOC, PDF, TXT or read online on Scribd
214 views

UTP CM Matlab Assignment and Test Questions - May 2015 Sem

This document contains 12 questions regarding computational methods and MATLAB assignments. The questions involve implementing numerical methods like bisection, false position, Newton-Raphson, and secant methods to solve nonlinear equations. Other questions involve solving systems of linear equations using Gaussian elimination, LU decomposition, Gauss-Seidel iteration and Thomas algorithm. Additional questions involve interpolation, numerical integration techniques like trapezoidal rule, Simpson's 1/3 rule, and solving ordinary differential equations using Euler's method and 4th order Runge-Kutta method. M-files are to be developed to implement these numerical techniques.

Uploaded by

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

FCM2043/FDM2043 - COMPUTATIONAL METHODS

JANUARY 2015 SEMESTER


MATLAB ASSIGNMENT/TEST QUESTIONS

Q1.

Develop an M-file to implement the bisection method. Using this program solve
the following problem.
The velocity of falling parachutist is given as
v (t )

gm
(1 e ( c / m )t ) .
c

Where v(t ) = velocity of parachutist = 40 m / s ,


g = gravitational constant = 9.8 m / s 2 ,
m = the mass of the parachutist = 68.1 kg .
Find the drag coefficient, c at the time t 10 seconds using the initial bracket of
the root as [13, 16] and iterate until a 0.001 %.
Q2.

Develop an M-file to implement the falseposition method. Using this program


solve the following problem.
The velocity of falling parachutist is given as
v (t )

gm
(1 e ( c / m )t ) .
c

Where v(t ) = velocity of parachutist = 40 m / s ,


g = gravitational constant = 9.8 m / s 2 ,
m = the mass of the parachutist = 68.1 kg .
Find the drag coefficient, c at the time t 10 seconds using the initial bracket of
the root as [13, 16] and iterate until a 0.001 %.
Q3.

Locate the root of f ( x) 2 sin( x ) x


(a)
Using the MATLAB function fzero with an initial guess of
x0 2 .
(b)
Using Newton-Raphson method by writing a function M-file. Use an
initial guess of x 0 0.5 and iterate until a 0.001 %.

Q4.

Develop an M-file to implement the modified secant method. Using this program
determine the loest positive root of f ( x) 8 sin( x ) e x 1 with an initial guess
of x0 0.3 and 0.01 . Iterate until a 0.000001% .

Q5.

Find the solution of the following set of linear algebraic equations

x 2 y 3z 1
3x 3 y 4 z 1
2 x 3 y 3z 2

(a) Using the left-division \.


(b) Using Gaussian elimination.
(c) Using the LU decomposition.
[Hint: Since [LU]x = b, let [U]x = y, so that [L]y = b. Now, first solve for
y and then for x.]
Q6.

Develop a function M-file


with the Thomas algorithm.
f1 g 1

e2 f 2 g 2

e3 f 3 g 3

.
.

Tridiag.m to solve the following tridiagonal system

.
.

e n 1

f n 1
en

g n 1

f n

x1 r1

x 2 r2
x3 r3

. .

. .

. .
x r
n 1 n 1
x n rn

Thomas Algorithm:
(i)
Decomposition:
e
ek k and
f k 1

f k f k ek . g k 1 , where k 2,3,4, , n .

(ii)
(iii)

Q7.

Forward substitution:

rk rk e k . rk 1 , where k 2,3,4, , n .

Back substitution:
r
xn n
fn
(rk g k . x k 1 )
and x k
, where k n 1, n 2, ,2,1 .
fk
Using your program, solve the following tridiagonal system.
x1 4.175
0.020875
2.01475

x
0.020875

2.01475
0.020875
0

0.020875
2.01475
0.020875
0
x3 =

0.020875
2.01475
2.0875

x4
Develop a MATLAB script file to determine the solution of the following system
of linear equations using the Gauss-Seidel iteration method by performing first
seven iterations.

9 x1 2 x2 3x3 2 x4 54.5
2 x1 8 x2 2 x3 3x4 14
3x1 2 x 2 11x3 4 x4 12.5
2 x1 3x 2 2 x3 10 x 4 21
Q8.

Develop a script M-file to estimate f ( 2.75) using Lagrange interpolating


polynomials of order 1, 2 and 3 for the following data.
x
f(x)

0
0

1
0.5

2
0.8

3
0.9

4
0.941176

5
0.961538

For each estimate find the true percent relative error if the try function is given by
f ( x)

Q9.

x2
.
(1 x 2 )

The force on a sailboat mast can be represented by the following function:


F

z 2.5 z / H

200
dz
e
7 z

where z the elevation above the deck and H the height of the mast.
Compute F for the case where H 30 using
(i)

the M-file for Trapezoidal rule with the step size h 0.1 .

the MATLABtrapz function.


Q10. Develop an M-file to implement Simpsons 1/3 rule. Using your program solve
the following problem.
The velocity of falling parachutist is given as
v(t )

gm
(1 e ( c / m ) t ) .
c

Where v (t ) = velocity of parachutist,


g = gravitational constant = 9.8 m / s 2 ,
m = mass of the parachutist = 45 kg ,
c = the drag coefficient = 32.5 kg / s .
If the distance, d, traveled by the parachutist is given by
6

v(t )dt ,
0

find the distance using Simpsons 1/3 rule for the segments 10, 20, 50, and 100.

Q11.

Develop an M-file for Eulers method to solve a first order ordinary differential
equation (ODE).

i(t)
E

The current around the circuit at time t is governed by the following differential
equation
3

di
2i 3e 2 t ,
dt

i (0) 2 .

Using your program, solve the above initial value problem over the interval from
t 0 to 2 with the step size h 0.1 .
Q12. Develop an M-file for Fourth-Order Runge-Kutta method to solve a first order
ordinary differential equation (ODE).
Using your program solve the following initial value problem over the interval
from x 0 to 2 with the step size h 0.2 .
dy
x2 y ,
dx

y (0) 0.8 .

You might also like