nr_iterate function

Compute Newton-Raphson Parameter Update with Numerical Stabilization