Numerical Methods

Check the file: Questions.pdf

• Questions.pdf

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

This is only a preview of the solution. Please use the purchase button to see the entire solution