Now using cannonical L-M formula.

The previous version is actually an approximation. I don't see why we should use
the approximate expression when the exact expression just requires two more
matrix multiplications. The multiplies take < O(n^3), while the least-squares
solution of the normal equations is O(n^3). So, it doesn't seem like there would
be much run-time cost. Will have to test in practice.