Function polynomial_fit(vector<double>&, vector<double>&, const double, double&, double&)

Function Documentation

void polynomial_fit(vector<double> &xa, vector<double> &ya, const double x, double &y, double &dy)

Performs a polynomial fit of a series of points CF Numerical Recipes p.

119 , Neville algorithm

Parameters:

  • xa – array of abcissas

  • ya – array of values

  • x – projected abcissa

  • y – requested extrapolated value

  • dy – requested error on extrapolated value