Index: ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp =================================================================== --- ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp (revision 17739) +++ ../trunk-jpl/src/c/shared/Numerics/NewtonSolveDnorm.cpp (revision 17740) @@ -31,7 +31,7 @@ while(true){ /*Newton step*/ - y2 = y1 - (c1*pow(pow(10.,y1),s-1.) + c2*pow(10.,y1) - c3)/((s-1)*c1*log(10)*pow(pow(10.,y1),s-1.) + c2*log(10.)*pow(10.,y1)); + y2 = y1 - (c1*pow(pow(10.,y1),s-1.) + c2*pow(10.,y1) - c3)/((s-1)*c1*log(10.)*pow(pow(10.,y1),s-1.) + c2*log(10.)*pow(10.,y1)); if( fabs(y2-y1)/(fabs(y2))