newtbyuxiang.cpp

来自「求解差值真的范围场馆哟个自行车不是你沙发是但书」· C++ 代码 · 共 32 行

#include<stdio.h>
#include<conio.h>
#include<math.h>

double newtbYX(double x[],double y[],int n,double t,double yp)
{
	int i,j;
	double w,r;
	x[n+1]=t;
	y[n+1]=yp;
		for(j=1;j<=n+1;j++)
	{
		for(i=n+1;i>=j;i--)
		{
			y[i]=(y[i]-y[i-j])/(x[i]-x[i-j]);
		}
	}
	t=y[n];
	w=1;
	for(i=0;i<=n;i++)
	{
		w=w*(t-x[i+1]);
	}
	r=t*w;
	return(r);
}