//牛顿迭代法!
/*
============================================================
题目:用牛顿迭代法求解3*x*x*x-2*x*x-16=0的近似解。
============================================================
*/
#include
#include
#define E 1e-8
double hs(double x)
{
return(3*x*x*x-2*x*x-16);
//原函数
}
double dhs(double x)
{
return(9*x*x-4*x);
//导函数
}
void main()
{
double x1=1.0,x2;
x2=x1+1.0;
while(fabs(x2-x1)>E)
//不管怎么变,仅仅要x1和x2能进入循环,就能够完毕求解。即x2!=x1.
{
x1=x2;
x2=x1-hs(x1)/dhs(x1);
//此处是用泰勒展式保留前两项推倒得出的。
读者可查找相关文献具体理解之。
}
printf("方程:3*x*x*x-2*x*x-16=0\n");
printf("解得:x=%.4lf\n",x2);
}
/*
============================================================
评: 牛顿迭代法非常好掌握。有用性较强!编程时仅仅需求得导函数。 将原函数和
导函数在程序中替换,就可以求得方程的解!
============================================================
*/