c中经典例题 用牛顿迭代法求实根

---

前言

c经典例题题目：用牛顿迭代法求根。方程为 ax³+bx²+cx+d=0，系数 a ， b ， c ， d的值依次为 1，2，3，4，由主函数输入。求 x 在 1 附近的一个实根。求出根后由主函数输出。

---

一、案例分析？

1.牛顿迭代法公式：x=x0-f(x0)/f’(x0)
2.首先确定其要在1附近求实根则确定x0=1
3.f‘（x）=3x² + 6x+3 第一次迭代=1-f(1)/f’(1)=1-10/12=0.1666666,以此类推求x的相似值。
4.直到|x-x0|<=10^-3

二、使用步骤

1.具体代码

代码如下（示例）：

#include<stdio.h>
#include<math.h>
int main()
{
    float solut(float a,float b,float c,float d);
    float a,b,c,d;
	printf("input a,b,c,d\n");
	scanf("%f,%f,%f,%f,\n",&a,&b,&c,&d);
	printf("x=%10.7f\n",solut(a,b,c,d));
	return 0;
}
float solut(float a,float b,float c,float d)
{
    float x=1,x0,f,f1;
    do{
        x0=x;
        f=((a*x0+b)*x0+c)*x0+d;
        f1=(3*a*x0+2*b)*x0+c;
        x=x0-f/f1;
    }
    while(fabs(x-x0)>=1e-3);
    return(x);
}

2.读入数据

代码如下（示例）：

1,2,3,4

---

结果

总结

首先我们先定义了一个名为solut的函数，形参a,b,c,d
在主函数main中格式化输入输出我们要输入和输出的结果
注意（输出结果时%10.7f是输出宽度为10精度为7的结果）
然后定义了一个do循环结束语为fabs(x-x0)>=1e-3(也就是|x-x0|>=10^-3)
循环中直接套用了牛顿迭代法的公式最后solut函数返回x。