Messing with the names a bit
An implementation of the Gauss-Newton algorithm

This is preparatory to trying out the Laplacian approximation in nonlinear mixed-effects models.