allow tighter tolerance in htol if requested by the user.

matlab/newrat.m

@@ -46,7 +46,7 @@ icount=0;
 nx=length(x);
 xparam1=x;
 %ftol0=1.e-6;
-htol_base = max(1.e-5, ftol0);
+htol_base = max(1.e-7, ftol0);
 flagit=0;  % mode of computation of hessian in each iteration
 ftol=ftol0;
 gtol=1.e-3;