#include <gsl/gsl_cdf.h>
#include <gsl/gsl_sf.h>
#include <stdio.h>

int main() {
  unsigned int ks[6] = {909922, 909954, 991071, 991073, 991074, 991076};
  double ts[6] = {909007.8, 909004.6, 990089.29, 990089.27, 990089.26,
                  990089.24};
  for (unsigned int i=0; i < 6; i++) {
    printf("%f\n", gsl_cdf_gamma_Q(ts[i], ks[i], 1.));
    printf("%f\n", gsl_cdf_poisson_P(ks[i]-1, ts[i]));
    printf("%f\n", gsl_sf_gamma_inc_P(ks[i], ts[i]));
  }
}