fleg.cpp

提供了几种数据拟合的c++代码：1 直线拟合(FIT)；2 线性最小二乘法((LFIT, COVSRT), (SVDFIT, SVDVAR)(oddity), (FPOLY, FLEG)(exampl

void fleg(double x, double pl[], int nl)
{
	pl[1] = 1.0;
	pl[2] = x;
	if (nl > 2)
	{
		double twox = 2.0 * x;
		double f2 = x;
		double d = 1.0;
		for (int j = 3; j<=nl; j++)
		{
			double f1 = d;
			f2 = f2 + twox;
			d = d + 1.0;
			pl[j] = (f2 * pl[j - 1] - f1 * pl[j - 2]) / d;
		}
	}
}