How would you set this problem up in MATLAB using RK4?

3. (20 points) For projectile motion in the absence of air resistance, the velocity of the projectile obeys the differential equation do dt where g is the acceleration due to gravity. If one is to take air resistance into account, the velocity of the projectile obeys the differential equation dt where u = , and c is a coefficient of air resistance. In terms of the velocity compo- nents, this vector equation can be written as u. dt dur Write or modify an Octave program that uses RK4 for systems to approximate the solution to this system with g 9.8m/s2 and c- 0.1 m-1