int gsl_sf_multiply_e(double x,double y,double* result);
double gsl_sf_multiply(double x,double y);
int gsl_sf_multiply_err_e(double x,double dx,double y,double dy,double* result);