Welch's t-test: Difference between revisions

(→‎{{header|C}}: forgot to include main function in update)
Line 95:
sum2 += (pow(h * i,a-1))/(sqrt(1-h * i));
}
double return_value = ((h / 6.0) * ((pow(x,a-1))/(sqrt(1-x)) + 4.0 * sum1 + 2.0 * sum2))/(expl(lgammal(a)+0.57236494292470009-lgammal(a+0.5)));//0.5723... is lgammal(0.5)
if ((isfinite(return_value) == 0) || (return_value > 1.0)) {
return 1.0;
Line 102:
}
}
 
//-------------------
int main(void) {
const double d1[] = {27.5,21.0,19.0,23.6,17.0,17.9,16.9,20.1,21.9,22.6,23.1,19.6,19.0,21.7,21.4};