mas01cr@280: #include <gsl/gsl_sf.h>
mas01cr@280: #include <stdio.h>
mas01cr@280: #include <stdlib.h>
mas01cr@280: #include <math.h>
mas01cr@280: 
mas01cr@280: int main(int argc, char *argv[]) {
mas01cr@280:   if(argc != 4) {
mas01cr@280:     fprintf(stderr, "Wrong number of arguments: %d\n", argc);
mas01cr@280:     exit(1);
mas01cr@280:   }
mas01cr@280: 
mas01cr@280:   long int meanN = strtol(argv[1], NULL, 10);
mas01cr@280: 
mas01cr@280:   double d = strtod(argv[2], NULL);
mas01cr@280:   double sigma2 = strtod(argv[3], NULL);
mas01cr@280: 
mas01cr@280:   double logw = (2 / d) * gsl_sf_log(-gsl_sf_log(0.99));
mas01cr@280:   double logxthresh = gsl_sf_log(sigma2) + logw
mas01cr@280:     - (2 / d) * gsl_sf_log(meanN)
mas01cr@280:     - gsl_sf_log(d/2)
mas01cr@280:     - (2 / d) * gsl_sf_log(2 / d)
mas01cr@280:     + (2 / d) * gsl_sf_lngamma(d / 2);
mas01cr@280: 
mas01cr@280:   printf("w: %f\n", exp(logw));
mas01cr@280:   printf("x_thresh: %f\n", exp(logxthresh));
mas01cr@280:   exit(0);
mas01cr@280: }