fleg.cpp

来自「工程算法 这是一个很有用的工程数值算法集锦」· C++ 代码 · 共 19 行

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;
		}
	}
}