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