xsq = SQRT( ABS(exf_one - huol*16. _d 0) )
x = SQRT(xsq)
& *( LOG( (exf_one + exf_two*x + xsq)
& *(exf_one+xsq)*0.125 _d 0 )
xsq = SQRT( ABS(exf_one - htol*16. _d 0) )
& *exf_two*LOG( exf_half*(exf_one+xsq) )
_RL xsq ! = x^2 [-]
xsq = max(sqrt(abs(1. _d 0 - 16. _d 0*huol)),1. _d 0)
x = sqrt(xsq)
& 2. _d 0*log(5. _d -1*(1. _d 0+xsq)) -
& (2. _d 0*log(5. _d -1*(1. _d 0+xsq)))
xsq = SQRT( ABS(one - huol*16. _d 0) )
& ( LOG( (one + two*x + xsq)*(one+xsq)*.125 )
xsq = SQRT( ABS(one - htol*16. _d 0) )
& + (one-stable)*( two*LOG(half*(one+xsq)) )
_RL huol,stable,xsq,x,psimh,psixh
xsq = max(sqrt(abs(exf_one - 16.*huol)),exf_one)
& ( LOG( (exf_one + exf_two*x + xsq)*
& (exf_one+xsq)*.125 _d 0 )
xsq = max(sqrt(abs(exf_one - 16.*htol)),exf_one)
& ( exf_two*LOG(exf_half*(exf_one+xsq)) )