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.