如何解决一道c语言数学计算

Question

不得不说，这个问题困扰了我很久，求这个问题的具体代码，经过了实践的代码会得到悬赏

Answer 1

一元三次方程求解代码如下：

#include<stdio.h>
#include<math.h>
int main() {
        double a, b, p, q, x1, m, n, i, w, s1, s2, s3;
    printf("请输入一元三次方程标准形式ax^3+bx^2+px+q=0\n");
    printf("a b p q的值，用空格隔开，其中a的值为1，b的值为0\n");
    scanf("%lf%lf%lf%lf", &a, &b, &p, &q);
    m = (3 * a * p - b * b) / (3 * a * a);
    n = (27 * a * a * q - 9 * a * b * p + 2 * b * b * b) / (27 * a * a * a);
    w = -0.5;
    s1 = -b / (3 * a);
    if (-n / 2 + sqrt(n * n / 4 + m * m * m / 27) > 0) 
    {
        s2 = pow((-n / 2 + sqrt(n * n / 4 + m * m * m / 27)), 1.0 / 3);
    }
    else 
    {
        s2 = -pow(-(-n / 2 + sqrt(n * n / 4 + m * m * m / 27)), 1.0 / 3);
    }
    if (-n / 2 - sqrt(n * n / 4 + m * m * m / 27) > 0) 
    {
        s3 = pow((-n / 2 - sqrt(n * n / 4 + m * m * m / 27)), 1.0 / 3);
    }
    else 
    {
        s3 = -pow(-(-n / 2 - sqrt(n * n / 4 + m * m * m / 27)), 1.0 / 3);
    }
    x1 = s1 + s2 + s3;
    printf("x1的解为：%f\n", x1);
}