newtb.c

求解差值真的范围场馆哟个自行车不是你沙发是但书

#include<stdio.h>
#include<conio.h>
#include<alloc.h>
double newtb(double x[],double y[],int n,double t)
{
int i,j;
double *s,p;
s=(double *)calloc(n+1,sizeof(double));
for(i=1;i<=n;i++)
s[i]=y[i];
for(j=1;j<=n-1;j++)
for(i=n;i>=j+1;i--)
s[i]=(s[i]-s[i-1])/(x[i]-x[i-j]);
p=y[n];
for(i<=n-1;i>=1;i--)
p=p*(t-x[i])+s[i];
free(s);
return(p);
}