,where a and b are the coefficients to find.

I'm using the following code to find $\displaystyle \theta = (a, b)$, and an initial value of: $\displaystyle \theta_0 = (0.000426, 1.53917)$

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 $\displaystyle x^p$, 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.