The solution of a constrained optimization problem is easily solved by
linearly approximating a point of interest to a first hypersurface given
as a constraint and finding the extreme value of the objective function
by moving a point r along a curved line having second-order osculation
with the geodetic line of a hypersurface S(r)=c and passing through the
point r.