int gsl_sf_pow_int_e(double x,int n,double* result);
double gsl_sf_pow_int(double x,int n);