Matlab: Write a "right-hand-side" m-file for projectile dynamics (with drag). Required function: 'AtmDensityEarth.m' included here.

Purpose of the file is to test modified eular and runge kutta functions on.

4. Projectile with drag % Test both methods for the projectile with drag model % Parameters and initial conditions: mu = 398600.4; % Earth gravitational constant (km ^3/s^2) Re-6378.14; % Earth equatorial radius (km) B = 1e6 ; rhoEun-@lh) AtmDensityEarth(h); % Density function (Be sure to get this m-file from the website.) x0 = 0; y0 = 0; vx0 = 30; vy0 = 40; T = 2*vy0/9.81; % simulation duration (sec) % Ballistic coefficient (kg/m^2) % initial x (horizantal) position (m) % initial y (vertical) position (m) % initial x (horizontal) velocity (m/s) % initial y (vertical) velocity (m/s) 4. Projectile with drag % Test both methods for the projectile with drag model % Parameters and initial conditions: mu = 398600.4; % Earth gravitational constant (km ^3/s^2) Re-6378.14; % Earth equatorial radius (km) B = 1e6 ; rhoEun-@lh) AtmDensityEarth(h); % Density function (Be sure to get this m-file from the website.) x0 = 0; y0 = 0; vx0 = 30; vy0 = 40; T = 2*vy0/9.81; % simulation duration (sec) % Ballistic coefficient (kg/m^2) % initial x (horizantal) position (m) % initial y (vertical) position (m) % initial x (horizontal) velocity (m/s) % initial y (vertical) velocity (m/s)