牛顿插值法

#include<stdio.h>
#include<conio.h>
#include<iostream.h>
double cs(double f[],double x[],int n)
{
double s=0.0,t=0.0;
int i,j;
for(i=0;i<=n;i++)
 {
  t=1.0;
  for(j=0;j<=i-1;j++)t=t*(x[i]-x[j]);
  for(j=i+1;j<=n;j++)t=t*(x[i]-x[j]);
  t=f[i]/t;
  s=s+t; if(i>n)break;
 }
 return s;
}

double N(double f[],double x[],int n,double y)
{
 double s=f[0],t=0.0;
 int i,j;
 for(i=1;i<=n;i++)
 {
  t=1.0;
  for(j=0;j<i;j++)t=t*(y-x[j]);
  s=s+cs(f,x,i)*t;
 }
 return s;
}

void main()
{
double x[3]={11.0,12.0,13.0},
       f[3]={0.190809,0.207912,0.224951},
       s;
printf("result=%f/n/n",N(f,x,2,11.5));
while(1)
{
 printf("Input a number(press 0 exit!):");
 scanf("%lf",&s);
 printf("Result=%f/n/n",N(f,x,2,s)); 
 if(s==0.0)break;
}
getch();
}

文章引用自：