cubic formula

i am trying to make a program that will compute the roots(real and complex number) of a cubic equation.
can anyone help me to develope this code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
#include<iostream>
#include<cmath>
using namespace std;

int main()
{
    double a, b, c, d, value;
    double f, g, h;
    double i, j, k, l, m, n, p, po;
    double r, s, t, u;
    double x1, x2, x2re, x2im, x3re, x3im, x3;
    
    cin >> value;
    for(int w=1; w <= value; w++){
       cin >> a >> b >> c >> d;
       cout << "CUBIC EQUATION : " << a << " x^3 + " << b << " x^2 + " << c <<" x + " << d << " = 0" << endl;
    
       f = ((3*c/a)-((b*b)/(a*a)))/3;
       g = ((2*(b*b*b)/(a*a*a))-(9*b*c/(a*a))+(27*d/a))/27;   
       h = ((g*g)/4)+((f*f*f)/27);
           
    if(f==0 && g==0 && h==0){     // all roots are real and equal
       x1 = pow((d/a),0.33333333333333333333333333333333);
       x2 = pow((d/a),0.33333333333333333333333333333333);
       x3 = pow((d/a),0.33333333333333333333333333333333);
       cout << "x = " << x1 << endl;
       cout << "x = " << x2 << endl;
       cout << "x = " << x3 << endl;
       }
    else if(h<=0){         // all 3 roots are real
       i = pow((((g*g)/4)-h),0.5);
       j = pow(i,0.33333333333333333333333333333333);
       k = acos((g/(2*i))*-1);
       l = j * -1;
       m = cos(k/3);
       n = sqrt(3) * sin(k/3);
       p = (b/(3*a))*-1;
       x1 = (2*j)*m-(b/(3*a));
       cout << "x = " << x1 << endl;
       x2 = l * (m+n) + p;
       cout << "x = " << x2 << endl;
       x3 = l * (m-n) + p;
       cout << "x = " << x3 << endl;
       }
    else if(h>0){        // only 1 root is real
       r = ((g/2)*-1)+pow(h,0.5);
       s = pow(r,0.33333333333333333333333333333333);
       t = ((g/2)*-1)-pow(h,0.5);
       u = pow((t),0.33333333333333333333333333333333);
       x1 = (s+u) - (b/(3*a));
       cout << "x = " << x1 << endl;
       x2re = (((s+u)*-1)/2) - (b/(3*a));
       x2im = -(s-u)*pow(3,0.5)/2;
       cout << "x = (" << x2re << "," << x2im << ")" << endl;
       x3re = (((s+u)*-1)/2) - (b/(3*a));
       x3im = (s-u)*pow(3,0.5)/2;
       cout << "x = (" << x3re << "," << x3im << ")" << endl;
       }
       }
    return 0;
}


using this link : http://en.wikipedia.org/wiki/Cubic_function#Cardano.27s_method

this should be the output

CUBIC EQUATION : -1 x^3 + 3 x^2 + 2 x + 7 = 0
x = 3.95367
x = (-0.476836,-1.24223)
x = (-0.476836,1.24223)


where my input is
-1 3 2 7
Last edited on
Topic archived. No new replies allowed.