#include <stdlib.h>
#include <stdio.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_sf_laguerre.h>

int main () {

    double x = 0.0, step = 0.2;
    int    k = 0, l = 0;

    gsl_sf_result res;

    for (k = 0; k < 4; k ++) {
        for (x = 4.0; x < 6.1; x += step)  {
            gsl_sf_laguerre_n_e(k,0.5,x,&res);
            printf("L_%d(%.1f) = %.8f        error = %.8f\n",k,x,res.val,res.err);
        }
        printf("\n");
    }


    gsl_sf_laguerre_n_e(2,0.5,4.97290e+00,&res);
    printf("L_%d(%.5f) = %.8f        error = %.8f\n",2,4.97290e+00,res.val,res.err);

    gsl_sf_laguerre_n_e(3,0.5,4.97290e+00,&res);
    printf("L_%d(%.5f) = %.8f        error = %.8f\n",3,4.97290e+00,res.val,res.err);

    return 0;
}