- xbt_assert0(x,"Please report this bug.");
- xbt_assert0( (x<0.0), "Please report this bug.");
- res_fpip = sqrt(1/x - 2/(3*var->df*var->df));
- return -(1/(2*x*x*res_fpip));
+ double critical_test;
+
+ xbt_assert0(var->df>0.0,"Don't call me with stupid values!");
+ xbt_assert0(x>0.0,"Don't call me with stupid values!");
+
+ res_fpip = 1/(var->df*var->df*x) - 2/(3*var->df*var->df);
+ xbt_assert0(res_fpip>0.0,"Don't call me with stupid values!");
+ critical_test = (2*var->df*var->df*x*x*sqrt(res_fpip));
+
+ return -(1.0/critical_test);