there is any criterion for "stpmax=STPMX*MAX(sqrt(sum),DP(n))"?
Is there any criterion to set the stpmax using following expression:

stpmax=STPMX*MAX(sqrt(sum),DP(n));

in routine newt()?