积累系统性知识
积聚技术精华
  首页    个人中心    撰写积文    建立课题    订立目标    整理积文    管理课题    管理目标    技能Get    代码积累 
牛顿插值法
error997 (error997)    2014-11-20 20:14:34      目标    课题
切换到: 纯代码  
   
#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();
}

文章引用自:

转自 http://blog.csdn.net/y___y/article/details/1788619
(+0)技能Get

建议楼主:搜索关键字 |参考其他资源 |回复 |追问
  error997(error997):   个人中心    课题    目标    代码积累