def mathcontour(x, y, errx, erry, cov):
    """
Suppose I want to draw the error ellipses for two parameters 'x' and 'y' 
with 1 s.d. uncertainties and the covariance between these uncertainties. 
I have a mathematica notebook that does that: 'chisq contour plot.nb'.

This method takes the values 'x,y,errx,erry,cov' and outputs Mathematica
code that I can just copy and paste in the appropriate notebook in order 
to draw the ellipses.
	"""
    print('x0=', x, ';')
    print('y0=', y, ';')
    print('\\[Sigma]x=', errx, ';')
    print('\\[Sigma]y=', erry, ';')
    print('\\[Sigma]xy=', cov, ';')
