Note: The subject is Matlab...

clc

clear

close all global g theta

g = 9.8; % m/s^2 - acceleration due to gravity theta = 60; Vo_arr = [60 80]; % m/s

for k =1:length(Vo_arr) Vo = Vo_arr(k);

% Compute intial velocity components Voy = Vo*sind(theta); Vox = Vo*cosd(theta);

% Computation of the time, tmax to attain maximum height tmax = Voy/g;

% Total flight time = ttol =2*tmax ttot = 2*tmax;

% Computation of maximium height hmax = (Voy)^2/(2*g);

% Let us compute the range, R of the projectile R = (Vox^2)*(sind(2*theta))/g;

% Generate arrays % The time array tarr = 0:0.002:ttot; % the instantaneous height and x arrays harr = Voy *tarr - 0.5 * g* tarr .^2; xarr = Vox *tarr;

% Plots plot(xarr,harr) hold on; grid on; end legend('60','80')

%%%%%%%

Test the script by running it (click on run). It will display a trajectory for a projectile launched at 60

o

with

speed array, Vo_arr of 60 and 80 m/s. This is for you to see the trend.

You will break up the script into two User-defined functions.

Using the script, create

(a). User-defined function that computes and outputs R, ttot, Vox and Voy using the speed, Vo as the

only input.

theta

and

g

are the only global variables. No other global variable should be declared. Note

that R stands for the range and ttot is the total time of flight.

(b). User-defined function that computes harr and xarr and plots(xarr, harr) using Vox,Voy and ttot as the

only input variables. This function should have no output

(d). Write a script that calls these two functions. Use for loops for 4 different Runs of the script for Vo =

70, 80, 100 and 125 m/s. These values should be supplied to the script via the keyboard using the input

command. Using legend label your plots. The plots should be on one figure which you should declare as

Figure 3

(e). Modify your program to use subplot to show each of the above plots on one page. Name the figure as

figure 4

.