Chapter 5, Problem 66P

This is a tilt and rotation question. Here are notes attached for reference.

I need help with a MATLAB code. For question b.6 I have the MATLAB code shown below. How do I edit the code to answer question b.7. Please make sure the plots are reasonable. clc; clear all;   % Constants mu = 398600; % Earth gravitational parameter, km^3/s^2   % Initial chief and deputy positions and velocities in ECI frame % Assume circular orbits in equatorial plane for simplicity a_c = 10000; % km a_d = 11500; % km   r_c0 = [a_c; 0; 0]; v_c0 = [0; sqrt(mu/a_c); 0];   r_d0 = [a_d; 0; 0]; v_d0 = [0; sqrt(mu/a_d); 0];   % Initial relative state delta_r0 = r_d0 - r_c0; delta_v0 = v_d0 - v_c0;   x0 = [delta_r0; delta_v0]; % 6x1 initial relative state   % Time span tspan = [0 3600]; % 1 hour in seconds   % Damping cases cases = struct( ... 'name', {'Critically damped', 'Under-damped', 'Over-damped'}, ... 'Kr', {eye(3)*2.5e-3, eye(3)*0.001, eye(3)*0.01}, ... 'P', {eye(3)*0.01, eye(3)*0.0006, eye(3)*0.02} ... );   % Simulate each case for i = 1:length(cases) Kr = cases(i).Kr; P =…

Just do Questions 7, 9, 11. Here are notes attached for reference.