Modify this code in MATLAB :https://www.mathworks.com/matlabcentral/fileexchange/2270-numerical-methods-for-physics--2e?focused=5047636&tab=function

Take the acceleration due to the external force to be in the x-

direction. Make the necessary changes only for the 4th order RK. In order to

check your program, use r =1 AU, v=2? AU/yr, ? = 0.005yr , and 2000 steps.

Plot both the orbit and a graph of the angular momentum as a function of

time.

12. Suppose that our comet is subjected to a constant force in one direction (eg gravitational attraction of a large but distant object).(85] Modify the orbit program to simulate this system. Set the strength of the perturbing force to be 1% of the initial gravitational force. Using fourth-order Runge-Kutta, show that an initially circular orbit is transformed into an elliptical orbit with the semimajor axis perpendicular to the perturbing force. Produce a graph of the angular momentum as a function time. Computer of 12. Suppose that our comet is subjected to a constant force in one direction (eg gravitational attraction of a large but distant object).(85] Modify the orbit program to simulate this system. Set the strength of the perturbing force to be 1% of the initial gravitational force. Using fourth-order Runge-Kutta, show that an initially circular orbit is transformed into an elliptical orbit with the semimajor axis perpendicular to the perturbing force. Produce a graph of the angular momentum as a function time. Computer of