#include <stdio.h>
#include <gsl/gsl_sf.h>
#include "../gammai.c"
#include <math.h>
#include <float.h>
#include "gammai.h"
#include "../msg.c"