atiken插值c++_c++插值代码

//数值计算实验 插值

#include &ltiostream&gt

#include &ltcmath&gt

#include &ltcstdlib&gt

#include &ltgsl/gsl_errno.h&gt

#include &ltgsl/gsl_spline.h&gt

#include &ltmylib/MyClock.h&gt

#include &ltmylib/InitData.h&gt

using namespace std;

//要插值的函数,为sin(x)，可改。

double f(double x)

{

return sin(x);

}

//产生函数区间为[a,b]上的n个点的插值函数数值，x,y分别存入横纵座标

bool InitFun(const double & a, const double & b, const int & n,

vector&ltdouble&gt & x, vector&ltdouble&gt & y)

{

//如果a大于b返回错误。

if (a &gt b)

{

return false;

}

//计算x值间隔大小，将[a,b]均分为n-1份，这样连端点，一共n个点

double gap = (b-a)/(n-1.0);

x.resize(n);

y.resize(n