Welch's t-test: Difference between revisions
Content added Content deleted
(→{{header|C}}: forgot to include main function in update) |
|||
Line 95: | Line 95: | ||
sum2 += (pow(h * i,a-1))/(sqrt(1-h * i)); |
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))); |
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))); |
||
if ((isfinite(return_value) == 0) || (return_value > 1.0)) { |
if ((isfinite(return_value) == 0) || (return_value > 1.0)) { |
||
return 1.0; |
return 1.0; |
||
Line 102: | Line 102: | ||
} |
} |
||
} |
} |
||
//------------------- |
|||
int main(void) { |
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}; |
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}; |