function plot_covar(mu,P,col)

A = sqrtm(P);
t = [0 : 2*pi/200 : 2*pi];
y = [cos(t); sin(t)];
x = 2*A*y + mu*ones(1,length(t));
plot(x(1,:), x(2,:),col);
