diff --git a/cephes/cprob/stdtr.c b/cephes/cprob/stdtr.c index 9e5c833..9f78a3e 100644 --- a/cephes/cprob/stdtr.c +++ b/cephes/cprob/stdtr.c @@ -1,3 +1,4 @@ +#include /* stdtr.c * * Student's t distribution @@ -171,7 +172,9 @@ else f += tz; j += 2; } + printf("f: %lf, x:%lf, z: %lf, rk:%lf\n", f, x, z, rk); p = f * x/sqrt(z*rk); + printf("p: %lf\n", p) } /* common exit */