C++ Wrapper of gsl non linear least squares Easy curve fit using gsl non linear least squares: auto result = curve_fit(gaussian, {1.0, 0.0, 1.0}, xs, ys); The example in the example.cpp file is adapted from the gsl webpage.