关闭

牛顿插值法

标签: input
2115人阅读 评论(0) 收藏 举报
分类:
#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();
}
 
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:329581次
    • 积分:3941
    • 等级:
    • 排名:第8009名
    • 原创:64篇
    • 转载:67篇
    • 译文:2篇
    • 评论:114条
    最新评论
    Opengl网站
    其它
    优秀编程网址