orbit formula graphing orbits

%Orbit One

thetadeg = 0:1:360; h = 51593.15284; mu = 398600; e = 0; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)

%Orbit Two

thetadeg = 0:1:360; h = 56935.55284; mu = 398600; e = 0.217819; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)

%Orbit Three

thetadeg = 0:1:360; h = 62277.95284; mu = 398600; e = 0.45708; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)

%Orbit Four

thetadeg = 0:1:360; h = 67620.35284; mu = 398600; e = 0.7177924; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)

%Escape thetadeg = 0:1:120; h = 72962.75284; mu = 398600; e = 1.00; r = (h^2/mu)./(1-e.*cosd(thetadeg)); polarplot(thetadeg,r)

[Steven Lord]

Rather than flagging your own question as Unclear, you should edit the question or add a comment with more information clarifying what you've asked. I've removed the flag.

Answer 1

[David Hill]

%Orbit One
 
thetadeg = 0:.1:360;
h = 51593.15284;
mu = 398600;
e = 0;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
hold on;                                    %need hold on
 
%Orbit Two
 

h = 56935.55284;
e = 0.217819;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
 
%Orbit Three
 

h = 62277.95284;
e = 0.45708;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
 
%Orbit Four
 

h = 67620.35284;
e = 0.7177924;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
polarplot(deg2rad(thetadeg),r)
 
%Escape                   need to reduce the range so it fits on the plot
%

h = 72962.75284;
e = 1.00;
r = (h^2/mu)./(1-e.*cosd(thetadeg));
Rad=deg2rad(thetadeg);
polarplot(Rad(484:3118),r(484:3118))