newton.cpp

数值计算方法中的,牛顿向前和向后算法. 以前花了很长时间编写的,希望对后来的师妹师弟有用.

float f(float x[],float y[],int s,int t)//牛顿插值法，用以返回差商
{
	if(t==s+1)   
		return (y[t]-y[s])/(x[t]-x[s]);
	else
		return (f(x,y,s+1,t)-f(x,y,s,t-1))/(x[t]-x[s]); 
}