HBUOJ--牛顿迭代法

HBUOJ–牛顿迭代法
牛顿迭代法
如何快速了解牛顿迭代法，请看👇
如何通俗易懂的讲解牛顿迭代法

迭代关系式：

用牛顿迭代法求根。方程为ax^3 + bx ^2+cx+d=0，系数a、b、c、d由主函数输入，求x在1附近的一个实根。求出根后，由主函数输出。
程序的运行示例如下：
输入方程的系数a、b、c、d：
1 2 3 4
方程是： 1.00x^3 + 2.00x ^2+ 3.00x+ 4.00=0
X= 0.0000176

输入数据格式："%f,%f,%f,%f"
输出数据格式："\n方程是：%5.2fx^3 +%5.2fx ^2+%5.2fx+%5.2f=0"
“\nX=%10.7f\n”

#include <stdio.h>
#include <math.h>
 
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-5);
 
    return x;
}                             
 
int main()
{                             
    float a, b, c, d;
 
    printf("\n输入方程的系数a、b、c、d：\n");
    scanf("%f,%f,%f,%f", &a, &b, &c, &d);
    printf("\n方程是：%5.2fx^3+%5.2fx^2+%5.2fx+%5.2f=0", a, b, c, d);
    printf("\nX=%10.7f\n", solut(a, b, c, d));
 
    return 0;
}