In attached file is the plot of data from gsl_fit_linear
after modified output printf() function:
printf ("%12.6f %12.6f %12.6f %12.6f \n", xf, yf, +1.0*yf_err, -1.0*yf_err );
//printf ("%12.6f %12.6f %12.6f %12.6f \n", xf, yf, yf + yf_err, yf - yf_err );
Regards,
Sergey.
