Question

For this homework you will be creating functions in SciLab. You will turn in electronically, each of the functions requested below. Part Problem Statement: Create a function called Projectile that will give you the X and Y location of an object while in flight. The cquations for the X and Y locations are given below function [X, Y]= Projectile (Velocity, Angle, Gravity) Input variables: Velocity velocity in units of m/sec Angle angle of the launcher from 0 to 90 Degrees Gravity gravity of the planet where the object is launched in units of m/sec X-the X location output as a column matrix in units of meters Y the Y location output as a column matrix in units of meters C C C Calculations: Use a WHILE loop to do the calculations while the Y value is greater than -0.01 Start your time at zero and increase in increments of 0.1 seconds You may want to use a variable and increment it to help make the output matrix. C C C Remember to initialize time, X, & Y, to zero o Part 2 Graphing Write a script in SCILAB that will graph the Projectile path using the points you calculated in part 1. Graph four different projectiles using Subplots Note this script should include the function from the first part and be used to calculate the points for the graphs Nar e the file HW7_Prob_01.sce For the 1st projectile use Velocity 90 m/sec, Angle 20*, Gravity = 9.81 m/sec For the 2nd projectile use Velocity = 90 m/sec, Angle = 30", Gravity = 9.81 m/sec* For the 3rd projectile use Velocity 90 m/sec, Angle 40', Gravity 9.81 m/sec For the 4th projectile use Velocity 90 m/sec, Angle 50'. Gravity 9.81 m/sec* For the legerd use the angle irforr.ation to distinguish the graphs . Remeber to include the graph title, label x and y valucs appropriately, include units, and grid lines

Answer 1

Scilab code for projectile motion clear all close all %1st projectile v90;angle-20;g=9.81; [X,Y]=Projectile (v,angle , g); subplot(4,1,1) plot (X,Y) legendFor angle 20) xlabelX in m.') ylabelY in m. ') title'Projectile motion') ylim(0 250]) xlim0 900]) 82nd projectile v90 angle 30;g-9.81 [X,Y Projectile(v,angle,g); subplot(4,1,2) plot (X,Y) legend For angle 30') xlabelX in m. ) ylabelY in m. ) titleProjectile motion ylim(0 250]) xlim(0 900]. %3rd projectile v-90; angle - 4 0; g-9.81; X,Y Projectile (v,angle,g) subplot(4,1l,3) plot (X,Y) legendFor angle 40') xlabelX in m. ') ylabelY in m.') titleProjectile motion') ylim(0 250]) xlim (0 900]) 4th projectile v-90; angle 50;g-9.81; X,Y-Projectile (v,angle,g); subplot(4,1,4) plot (X,Y) legendFor angle 50' m.) XIаbe. ylabelY in m.') title'Projectile motion ylim0 250 ] ) xlim( 0 900]

Function for projectile motion function X, Y]=Projectile (Velocity,angle,gravity) vo-velocity th-(pi/180)*angle; t 0; k-1; yy-1; X(1)-0;Y(1)0; whi . ; k-k+1; X (k)-v0*t*cos ( th ) ; Y(k) v0*tsin ( th)-0.5*gravity" t^2; yy=Y (k) end end 응용용용 용 응용용 End of Code 용용응용음용용용용웅응용용 Projectile motion E 200 For angle 20 100 100 200 300 400 500 600 700 800 900 X in m Projectile motion 200 or angle 30 600 0 100 400 500 200 300 700 800 D00 X in m. Projectile motion E 200 For andle 40 드 100 300 500 600 700 200 0 100 400 800 900 Projectile motion 200 For angle 50 100 100 200 500 900 300 400 600 700 800 X in m. Published with MATLAB R2018a 2

%%Scilab code for projectile motion
clear all
close all

%1st projectile
v=90;angle=20;g=9.81;
[X,Y]=Projectile(v,angle,g);
figure(1)
subplot(4,1,1)
plot(X,Y)
legend('For angle 20')
xlabel('X in m.')
ylabel('Y in m.')
title('Projectile motion')
ylim([0 250])
xlim([0 900])

%2nd projectile
v=90;angle=30;g=9.81;
[X,Y]=Projectile(v,angle,g);
subplot(4,1,2)
plot(X,Y)
legend('For angle 30')
xlabel('X in m.')
ylabel('Y in m.')
title('Projectile motion')
ylim([0 250])
xlim([0 900])

%3rd projectile
v=90;angle=40;g=9.81;
[X,Y]=Projectile(v,angle,g);
subplot(4,1,3)
plot(X,Y)
legend('For angle 40')
xlabel('X in m.')
ylabel('Y in m.')
title('Projectile motion')
ylim([0 250])
xlim([0 900])

%4th projectile
v=90;angle=50;g=9.81;
[X,Y]=Projectile(v,angle,g);
subplot(4,1,4)
plot(X,Y)
legend('For angle 50')
xlabel('X in m.')
ylabel('Y in m.')
title('Projectile motion')
ylim([0 250])
xlim([0 900])

%%Function for projectile motion
function [X,Y]=Projectile(Velocity,angle,gravity)

    v0=Velocity;
    th=(pi/180)*angle;
    t=0; k=1; yy=1;
    X(1)=0;Y(1)=0;
    while yy>-0.01
        t=t+0.1;
        k=k+1;
        X(k)=v0*t*cos(th);
        Y(k)=v0*t*sin(th)-0.5*gravity*t^2;
        yy=Y(k);
    end
end

%%%%%%%%%%%%%%% End of Code %%%%%%%%%%%%%