Write a Matlab program:

%% 2 . Write an m-file that implements the 4th order Runge-Kutta method. % do another quick check rhs @(t,x) -x ; % d/dt(x) --x the solution is: x(t) = x(0)*exp(-t) time -0:.01:1; x0 = 1; [~,xout1] = ode45(rhs,time , x0); xout2 - ODENumIntModEuler(rhs,time,x0) xout3 =0DENumIntRK4(rhs,time , x0 ) ; % just as before % the exact solution at the end time is: exp(-1) 0.367879441171442 % the last element of xout3 should be close to this value