C语言 牛顿迭代法求根

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

#include <stdio.h>
#include <math.h>

// 牛顿迭代法
double newton_raphson(double a, double b, double c, double d) {
    double x = 1.0;  // 初始猜测
    double tolerance = 1e-6;
    double diff = 1.0;
    
    while (diff > tolerance) {
        double fx = a*x*x*x + b*x*x + c*x + d;
        double fpx = 3*a*x*x + 2*b*x + c;
        double new_x = x - fx / fpx;
        diff = fabs(new_x - x);
        x = new_x;
    }
    
    return x;
}

int main() {
    double a = 1, b = 2, c = 3, d = 4;
    double root = newton_raphson(a, b, c, d);
    
    printf("方程的一个实根为：%.6f\n", root);
    
    return 0;
}

解释说明：

牛顿迭代法：

• 函数 newton_raphson 使用牛顿迭代法求解三次方程 ax^3 + bx^2 + cx + d = 0 在 x = 1 附近的一个实根。
• 主函数输入方程的系数 a, b, c, d，调用 newton_raphson 函数并输出求得的根。