* which is to use the bounds in Chick,Branke,Schmidt together with
* interval arithmetic.
*/
double z = -StrictMath.abs(s);
StudentDist dist = new StudentDist((int) dof);
double u = dist.cdf(z);
if (s <= 0) // s=z
return StrictMath.log(u);
else
// s=-z
return StrictMath.log1p(-u);