% PLot for # 3 in HW 2, Spring 2013

clear all

ee=0.093;  % eccentricity of Mars
D=1;       % arbitrary parameter
r=ee*D;  e2=ee^2;

t=[0: 0.01: 2*pi];
st=sin(t);  ct=cos(t);

figure(231);
plot(r/(1-e2)*ct, r/sqrt(1-e2)*st,'k--',  r/(1-e2)*ct, r/(1-e2)*st, 'r')
hold on
plot(-e2*D/(1-e2),0,'ko',  0,0,'ro','markersize',6)
hold off
axis('square')
legend('ellipse', 'circle', 'focus', '(0,0)')
set(gca, 'fontsize', 20)

