收集的C语言算法程序

#include<iostream>
using namespace std;
#include<math.h>
#define n 3
#define N 50
#define err 0.00001

double f(double *a,double x)                            //指定f的系数，求出其值
{
int i;
double r=0;
for(i=0;i<=n;i++)
r+=a[i]*pow(x,n-i);
return r;
}

double ff(double *a,double x)                      //求f的导函数的值
{
int i;
double r=0;
for(i=0;i<=n;i++)
r+=(n-i)*a[i]*pow(x,n-i-1);
return r;
}

void main()
{
int i;
double x[N],a[n+1];
cout<<"Please input the initialized x0:\n";
cin>>x[0];
cout<<"Please input the matrix of coefficients:\n";
for(i=0;i<=n;i++)
cin>>a[i];

for(i=0;i<N;i++)
{
x[i+1]=x[i]-f(a,x[i])/ff(a,x[i]);                      
if(fabs(x[i+1]-x[i])<err) {cout<<"The root is:\n"<<x[i+1]<<"\n";break;}
}
}