Consider the general 2D Gaussian function, centered at (0.5,0.5),
A*exp(-a*(-0.5 + x)**2-b*(-0.5 + x)*(-0.5 + y)-c*(-0.5 + y)**2)
where the covari