牛顿迭代法和二分法求方程的根

使用牛顿迭代法（在1.5附近）和二分法（在（-10，10）区间内）求2*x*x*x-4*x*x+3*x-6=0的根

牛顿迭代法：

算法分析：假设方程f（x）=0，在x0附近有一个根，那么用xn+1=xn-f(xn)/f’(xn)。依次计算x1，x2，x3，…………，那么这个序列无限逼近方程的根。原理是根据f（x）在x0附近的值和斜率，估计f（x）和x轴的交点。

#include<stdio.h>
int main()
{
double x=1.5;
double fx,f1x,eps=0.00001;
do
{
fx=2*x*x*x-4*x*x+3*x-6;
f1x=6*x*x-8*x+3;
x=x-fx/f1x;
if(fx<0)fx=-fx;
if(fx<eps)break;
}while(1>0);
printf("%lf",x);
return 0;
}

二分法：

算法分析：对于区间[a，b]上连续不断且f（a）·f（b）<0的函数y=f（x），通过不断地把函数f（x）的零点所在的区间一分为二，使区间的两个端点逐步逼近零点，进而得到零点近似值的方法。

#include<stdio.h>
int main()
{
double l=-10.0,h=10.0;     //划定区间，默认函数在-10出取负值，在10处取正值
double m,fm,eps=0.000001;  //规定逼近值eps
do
{
m=(l+h)/2;         //找到区间中点
fm=2*m*m*m-4*m*m+3*m-6;    //计算区间中点函数值
if(fm<0)l=m;       //如果区间中点函数值为负值，将区间最小端点赋为区间中点
else h=m;          //如果区间中点函数值为正值，将区间最大端点赋为区间中点
if(fm<0)fm=-fm;    //如果函数值为负值化为正值便于和eps比较，此处和上述任何一个步骤无关
if(fm<eps)break;   //如果函数值趋近与eps，即得到根m
}while(1);
printf("%lf",m)
return 0;
}

91b1