function h = plot_ellipse(xc,yc,Q,color)

theta = 0:2*pi/60:2*pi;

Ellipse = sqrtm(Q)*[ cos(theta) ; sin(theta)];

h = plot(Ellipse(1,:)+xc,Ellipse(2,:)+yc,color);


