motion of projectile in matlab

main.m
% Define initial conditions of the projectile
v0 = 20;           % initial velocity (m/s)
theta = 45;        % launch angle (degrees)
h0 = 10;           % initial height (m)

% Set up time vector and time step
t = 0:0.01:10;     % simulation time vector (s)
dt = t(2) - t(1);  % time step (s)

% Define gravitational acceleration and air resistance (if applicable)
g = 9.81;          % gravitational acceleration (m/s^2)
rho = 1.2;         % air density (kg/m^3)
Cd = 0.47;         % drag coefficient
A = pi*(0.1)^2;    % cross-sectional area (m^2)

% Define function to calculate the equations of motion
f = @(t, Y) [Y(2); -g - (0.5*rho*Cd*A/m)*norm(Y(2))*Y(2)];

% Solve the equations of motion numerically
Y0 = [h0; v0*cosd(theta); 0; v0*sind(theta)];
[t,Y] = ode45(f, t, Y0);

% Plot trajectory of projectile over time
plot(Y(:,3), Y(:,1));
xlabel('Distance (m)');
ylabel('Height (m)');
title('Projectile Trajectory');

% Plot velocity of projectile over time
plot(t, sqrt(Y(:,2).^2 + Y(:,4).^2));
xlabel('Time (s)');
ylabel('Velocity (m/s)');
title('Projectile Velocity');
1078 chars
34 lines