#include <stdio.h>

#include <gsl/gsl_sf_exp.h>


void test_x_y_e10(double x, double y)
{
    gsl_sf_result_e10 result;
    int status;

    printf("\nx=%f, y=%f\n", x, y);

    status = gsl_sf_exp_e10_e(x, &result);
    printf("gsl_exp_e10:      (%.15g +/- %.15g) * 10^%d\n", result.val, result.err, result.e10);

    status = gsl_sf_exp_mult_e10_e(x, y, &result);
    printf("gsl_exp_mult_e10: (%.15g +/- %.15g) * 10^%d\n", result.val, result.err, result.e10);
}

int
main(void)
{
    test_x_y_e10(10000, 1);
    test_x_y_e10(100, 1);
}