,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.