,where a and b are the coefficients to find.

I'm using the following code to find

, and an initial value of:

So here's the problem. In the first iteration, the Jacobian grows to be a very badly scaled matrix:

note that the left column is

, which makes for a huge number compared to the right column, making the matrix very badly scaled, so i never get to the solution.

I've tried everything from QR decomposition to logarithmic scaling, but nothing seems to work for me.

Any help will be much appreciated.