newton_honer.cpp

数值计算中的一些算法用C语言实现

#include <iostream.h>
#include <math.h>

float Honer(float x0)
{
int i;
float a[5],b[5],d[4];
a[0]=-4.0;
a[1]=3.0;
a[2]=-3.0;
a[3]=0.0;
a[4]=2.0;
b[4]=a[4];
d[3]=a[4];
   for(i=3;i>=1;i--)
   {
   b[i]=x0*b[i+1]+a[i];
   d[i-1]=x0*d[i]+b[i];
   }
   b[0]=x0*b[1]+a[0];
   return b[0]/d[0];
}
void main()
{
	float x0,x1=3.0f;
	do
	{
		x0=x1;
		x1=x0-Honer(x0);
	}while(fabs(x0-x1)>1e-6);
	cout<<"2*x^4-3*x^2+3*x-4=0的解是:"<<x1<<endl;    
}