Code illustrating that the gradient of a function f defining
an implicit surface S = {p : f(p) = C} gives the direction of
the normal. The normals are estimated numerically in the
above example, although they could have been computed
analytically by using the chain rule.