牛顿插值法
来源:互联网 发布:unity3d ui 高科技 编辑:程序博客网 时间:2024/05/01 14:20
#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;
}
#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;
}
{
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();
}
{
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();
}
文章引用自:
- 牛顿插值法
- 牛顿插值法
- 牛顿插值法
- 牛顿插值法
- 牛顿插值法
- MATLAB 牛顿插值法程序
- matlab牛顿插值法
- 牛顿
- 牛顿插值法 matlab版算法
- 牛顿插值法 C语言版算法
- 插值法:拉格朗日插值、牛顿插值
- 《拉格朗日和牛顿插值法》 …
- 【数值分析】插值法:拉格朗日插值、牛顿插值
- 牛顿插值法的C语言实现
- uva 1546 - Complete the sequence!(牛顿插值法)
- 【数值分析】插值法:拉格朗日插值、牛顿插值
- 利用均差的牛顿插值法(Newton)
- 利用差分的牛顿插值法(Newton)
- 快速傅立叶变换程序与反变换程序
- 单反相机
- ERP面试题之解答[作者:Yuanhy]
- Guass消元法和Guass列主元素消元法解线性方程组
- 追赶法解三对角阵方程组
- 牛顿插值法
- 解决关于页面空间disable之后服务器端取值问题
- 最小二乘法程序
- 处理JQuery和My97DatePicker的冲突
- mount/umount 详解
- 在ubuntu 下升级内核就是这么简单。。。。
- Windows的远程注册表服务(RemoteRegistry)的两个隐藏的依赖关系
- 推荐曹济的FPA培训课程
- unicode文本文件存储类(C++)