This article deals with a problem of the contour lines of cartogrpahic distortions computation. This method is based on the Bezier cubic surface and ensures C2-continuity between surfaces.
Numerical derivatives are computed with the Stirling formula.