C语言编程：牛顿迭代法求方程的根

程序流程分析：

① 赋值x0＝1.5,即迭代初值；

② 用初值x0代入方顷游中程中计算此时的f(x0)及f’(x0)，程序中用变量f描述方磨启程的值，用fd描述方程求导之后的值；

③ 计算增量d=f/fd；

④ 计算下一个x，x=x0-d

⑤ 把新产生的x替换x0，为下一次迭代做好准雀山备

⑥ 若d绝对值大于1e-3，则重复②③④⑤步。

源程序代码：

#include <math.h>

main()

{

float x,x0,d,f,fd

x0=0

do {

f=2*x0*x0*x0-4*x0*x0+3*x0-6

fd=6*x0*x0-8*x0+3

d=f/fd

x=x0-d

x0=x

}while(fabs(d)>1e-3)

printf("x=%f\n",x)

}

x=0,f(x)=-6,f1(x)=-3,x1=-2,

x=-2,f(x)=-46,f1(x)=39,f(x)/f1(x)=-1.179

因好迹此，在第一次循环友高并是就会跳出。把fabs(f(x)/f1(x))>1e-5改成fabs(f(x)/f1(x))<念乎1e-5试试（针对的是问题补充）

你的算法用的有问题

#include <stdio.h>

#include <math.h>

int main(void)

{

    double a=0.0

    double b=0.0

    double c=1.0

    double x//----

    const double e=0.00001

    printf("input num:\n")

    scanf("%lf",&a)

    getchar()

    x=a //负数时要进行处理

    if ( a< 0 ) a=-a

    c=1,b=a

    while(fabs(c-b)>e)

    {

        c=b

        b=0.5*(c+a/c)

    }

    if(x<0)//----

        printf("sqrt(%lf) = %.3lfi\n",a,b)

    else

        printf("sqrt(%lf) = %.3lf\n",a,b)

    getchar()

    return 0

}

---