## Question

## Transcribed Text

## Solution Preview

This material may consist of step-by-step explanations on how to solve a problem or examples of proper writing, including the use of citations, references, bibliographies, and formatting. This material is made available for the sole purpose of studying and learning - misuse is strictly forbidden.

function main()clear all

close all %empty variables in memory then close all windows, pictures etc start fresh

options=odeset('RelTol',1e-12,'AbsTol', 1e-12);

time1=8147; %time it takes for earth to do one revolution in seconds

mu=3.986004415*10^5;

r0=[6678.1363, 0, 0];

v0=[0, 7.725760634075587,0];%

[propagatedtime,propagatedposvel] =ode45(@(t,x)twobody(t,x,mu),[0, time1],[r0, v0]',options);

N=length(propagatedtime); %length of the vector propagatedtime

size(propagatedposvel)

plot3(propagatedposvel(:,1),propagatedposvel(:,2),propagatedposvel(:,3))

% we plot the values x,y,z in 3d

hold on

xlabel('x')

ylabel('y')

zlabel('z')

disp('-------Part 3----------')

options=odeset('RelTol',1e-12,'AbsTol', 1e-12);

a=propagatedposvel(N,:) % propagatedposvel is of size N times 6...