Fix the (header of) hessian routine

Description of the algorithm seems wrong (check the reference to the formulas given in Abramovitz and Stegun).