C = [0 0] ; % center
a = 2 ; % major axis
e = 0.4 ; % eccentricity
b = a*sqrt(1-e^2) ; % minor axis
%
R = 2; % Radius if cricle
th = linspace(0,2*pi) ;
% Ellipse
xe = C(1)+a*cos(th) ;
ye = C(2)+b*sin(th) ;
% Circle
xc = C(1)+R*cos(th) ;
yc = C(2)+R*sin(th) ;
% plot
plot(xe,ye,'b',xc,yc,'r')
axis equal