function PlotEllipse(m,a,b)
% DRAWELLIPSE plot an ellipse
%   DrawEllipse(m,a,b) plot an ellipse with center in m and
%   excentricities a and b.
    
theta=[0:0.02:2*pi];
u=[m(1)+a*cos(theta); m(2)+b*sin(theta)];
plot(u(1,:),u(2,:));

    