Numerical Methods (Lab.
)
GENG 300
Course Number: GENG 300
Course Name: Numerical Methods Lab
Summer 2023
Lab Assignment #2
Group Number in Blackboard: 4
Student Name 1 : Abdallah Al-Amiri
ID Number 201901683
Student Name 2 : Murshed Al-Muhanadi
ID Number 201706102
Student Name 3 : Khalifa Al-Malki
ID Number 201803375
Instructor: Eng. Nasser Alnohmi
Nasser Al-Nohmi Pa ge|1
Numerical Methods (Lab.)
GENG 300
A3-1. Find out the velocity of a free-falling bungee jumper using subfunction
structure, which is defined as:
𝑔𝑚 𝑔 𝐶𝑑
𝑣=√ 𝑡𝑎𝑛ℎ (√ 𝑡)
𝐶𝑑 𝑚
m=70 g, t=10 s, cd=0.30, g=9.81 m/s2
Hints: Find out the velocity using subfunction and call it from main function.
Code:
Output:
Nasser Al-Nohmi Pa ge|2
Numerical Methods (Lab.)
GENG 300
A3-2. Modify If.. Structure as used in exercise 3.3.1 by input and disp function. Test
the code by changing the variable ‘grade’
Code:
Output:
Nasser Al-Nohmi Pa ge|3
Numerical Methods (Lab.)
GENG 300
A3-3. For a scalar, the built-in MATLAB sign function returns the sign of its
argument (−1, 0, 1). Here’s a MATLAB session that illustrates how it works:
>> sign(25.6)
ans =
1
>> sign(-0.776)
ans =
-1
>> sign(0)
ans =
0
Develop an M-file to perform the same function by using if structure.
Code:
Output:
Nasser Al-Nohmi Pa ge|4
Numerical Methods (Lab.)
GENG 300
A3-4. Write a program using for loop to compute factorial. Suppose 5!=1x2x3x4x5.
Hints: x=1;
for i=1:n
x=x*i;
Code:
Output:
Nasser Al-Nohmi Pa ge|5
Numerical Methods (Lab.)
GENG 300
A3-5. Write a program to find out positive values from the matrix below:
10 −3 5
A=[−5 2 8 ]
4 0 −7
Hints:
[m,n]=size(A);
for i=1:m
for j=1:n
if A(i,j)……….
Code:
A = [10 -3 5; -5 2 8 ;4 0 -7]
positive_values = A(A > 0);
disp("Positive values in matrix A:");
disp(positive_values);
Output:
Nasser Al-Nohmi Pa ge|6
Numerical Methods (Lab.)
GENG 300
A3-6. Calculate value of k2-50 for all integers in [-10, 10] domain but only until k2-50
becomes negative. Use while --- break structure and fprintf for displaying results.
Code:
k = -10;
while true
result = k^2 - 50;
if result < 0
break;end
fprintf('k = %d, k^2 - 50 = %d\n', k, result);
k = k + 1;
end
Output:
Nasser Al-Nohmi Pa ge|7
Numerical Methods (Lab.)
GENG 300
A4-1. Using if…..elseif else structure, write a program for the following grading
policy of a University.
A (Excellent) = 90 to 100
B+ (Very Good) = 85 to <90
B (Good) = 80 to <85
Fail = <80
Note: If there is a negative value as an input, it will show an error message “negative
value is not allowed “.
Execute the program and display the results, whenever the grade is 90, -89, 85, and 79
Code:
function Grade
grade=input('Enter Your Gade: ');
if grade > 100
error('Invalid');
elseif (grade >=90 && grade <= 100)
disp('A (Excellent)');
elseif (grade >= 85 && grade < 90)
disp('B+ (Very Good)');
elseif (grade >= 80 && grade < 85)
disp('B (Good)');
elseif grade >= 0
disp('Fail');
else
error('negative value is not allowed');
end
end
Nasser Al-Nohmi Pa ge|8
Numerical Methods (Lab.)
GENG 300
Output:
A4-2. Create an index of a book by using Switch Structure with the condition below:
Chapter 1: Introduction
Chapter 2: Modeling
Chapter 3: Programming
If the chapter is not between 1-3, it will show an error message “The Index
number is not vailed”.
Execute the program and display the results, whenever the chapter number is 1, 2, 3,
and 4
Code:
function Book
chapterNo=input('Enter chapter number: ');
switch(chapterNo)
case 1
disp('Chapter 1: Introduction');
case 2
disp('Chapter 2: Modeling');
case 3
disp('Chapter 3: Programming');
otherwise
error('The Index number is not vailed');
end
end
Nasser Al-Nohmi Pa ge|9
Numerical Methods (Lab.)
GENG 300
Output:
A4-3.
(a) Plot the bungee jumper velocity by using fplot over the range from t=0 to 12 sec
by using @.
9.81 ∗ 80 9.81 ∗ 0.25
𝑣=√ 𝑡𝑎𝑛ℎ (√ 𝑡)
0.25 80
(b). Pass the function to the m-file and find out average velocity from t=1 to 12.
Code1: fplot
Code2: Passing the func. To m-file and finding the avg of velocity
function Vavg(f,xl,xu,step)
t=xl:step:xu;
v=f(t)
V_average= mean(v)
end
Nasser Al-Nohmi P a g e | 10
Numerical Methods (Lab.)
GENG 300
Output1:
Output 2:
Nasser Al-Nohmi P a g e | 11
