horner

来自「收集的C语言算法程序」· 代码 · 共 28 行

Horner算法2007-07-16 18:32#include<iostream>
using namespace std;
#define N 5                 

void main()
{
int i;
double x;
double a[N+1],b[N+1],c[N+1],d[N+1];
cout<<"请按降阶输入多项式的系数:\n";
for(i=0;i<=N;i++) cin>>a[i];
cout<<"请输入x0:\n";
cin>>x;
b[N]=a[N];
c[N]=N*a[N];
d[N]=a[N]/(N+1);
for(i=N-1;i>=0;i--)
{
b[i]=a[i]+x*b[i+1];
     c[i]=i*a[i]+x*c[i+1];
d[i]=x*a[i]*1/(i+1)+x*d[i+1];
}
cout<<"f(x)的在x0处的值为："<<b[0]<<"\n";
cout<<"f(x)的在x0处的一阶导数值为："<<c[0]<<"\n";
cout<<"f(x)的在x0处的不定积分值为（积分常数为0）："<<d[0]<<"\n";
}