covariance ellipse

plot covariance ellipse according to x'*inv(P)*x=c