牛顿插值法

#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();
}
 
阅读更多
文章标签: input
个人分类: 数值计算程序专栏
上一篇追赶法解三对角阵方程组
下一篇最小二乘法程序
想对作者说点什么? 我来说一句

牛顿插值法的c++程序

2009年01月04日 47KB 下载

拉格朗日插值法 牛顿插值法

2009年09月20日 1KB 下载

matlab实现牛顿插值法

2011年12月13日 612B 下载

数值分析 牛顿插值法c++代码

2010年10月10日 746B 下载

没有更多推荐了,返回首页

关闭
关闭