Question
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...