newton.cpp

大量数值计算的算法

#include<stdio.h>
#include<math.h>

double f(double x)
{
	return (676039*pow(x,12)-1939938*pow(x,10)+
		2078505*pow(x,8)-1021020*pow(x,6)+225225*
		pow(x,4)-18018*x*x+231);
}

double f_dao(double x)
{
	return (8112468*pow(x,11)-19399380*pow(x,9)+
		16628040*pow(x,7)-6126120*pow(x,5)+900900*
		pow(x,3)-36036*x); 
}

main()
{
	double a=-1,b=1,h=0.01,a1,b1,x0,x1,op,i;
	a1=a;
	b1=a1+h;
	while(b1<b)
	{
		while(f(a1)*f(b1)>0)
		{
			a1=b1;
			b1=a1+h;
		}
		x0=0.5*(a1+b1);
		do
		{
			x1=x0-f(x0)/f_dao(x0);
			i=x1-x0;
			x0=x1;
		}while(fabs(i)>0.5e-10);
		printf("x=%1.9f\n",x1);
		a1=b1;b1=a1+h;
	}
}