(元記事 2012/02/05) #include <gsl/gsl_interp.h>#include <gsl/gsl_spline.h> gsl_interp_accel *acc = gsl_interp_accel_alloc();const gsl_interp_type *T = gsl_interp_cspline;gsl_spline *spline = gsl_spline_alloc(T, data_size); gsl_spline_init(spline, f0, fn, data_size); i=0;for(xi = inif[0]; xi < inif[data…