I'm hoping someone can help check whether my final contour plots look plausible based on the surface. I haven't done too much differential geometry but I've needed to work with Gaussian/Mean curvature for a simple 3D gaussian surface. Here's an example: (A = 7, a=b=1/(3.5)^2) It's parameterised using Monge patch representation; the z axis is f(x,y) = Ae^-(ax^2+by^2). The Gaussian curvature K and Mean curvature H can be approximated via the following: A little bit of calculation later, I end up with: and These don't look right to me, considering on my example surface there's clearly curvature beyond x=y~4 range. Also, should the Gaussian and Mean contours look more similar? If it seems like my method is sound yet these still look wrong then the error must be in my calculation code. I haven't been able to spot a mistake in it, though. Many thanks.