24 double dif = fabs(x-xPoints[0].first),dift;
25 const unsigned int _K=xPoints.size(),_KP=_K+1;
26 std::vector<double>c(_KP),d(_KP);
28 for (
unsigned int i=0;i<_K;i++) {
29 dift=fabs(x-xPoints[i].first);
34 c[i]=d[i]=xPoints[i].second;
36 y = xPoints[ns--].second;
37 for (
unsigned int m=0;m<_K-1;m++) {
38 for (
unsigned int i=0;i<_K-m-1;i++) {
39 double ho = xPoints[i].first-x;
40 double hp= xPoints[i+m+1].first-x;
45 <<
"Error in polynomial extrapolation"
51 deltay = 2*(ns+1) < (
int)(_K-m-1) ? c[ns+1]: d[ns--];