0

カルダノの方法を使用して 3 次関数の根を計算するプログラムを作成しようとしています

ここに私のコードがあります:

#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){
       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;
}

このコードでユーザー定義の ComplexNumber を実装するのを手伝ってくれる人はいますか? このリンクhttp://en.wikipedia.org/wiki/Cubic_functionを使用したいのですが、これを理解できません。

4

1 に答える 1

3

このコードでユーザー定義の ComplexNumber を実装するのを手伝ってくれる人はいますか?

ユーザー定義の ComplexNumber 型を実装しないでください。言語で提供されているものを使用してください。ただ#include <complex>。これにより、単純に複雑な変数を持つことができますcomplex<double> variable_name。

コードのコメント:

  1. std::sqrt(x)よりも使ったほうがいいですstd::pow(x, 0.5)。
  2. POSIX 準拠のマシンを使用している場合、数学ライブラリのヘッダーに cbrt(double) という立方根関数があります。(C++ ヘッダーにエクスポートされない場合があります。) これも より好ましいstd::pow(x,0.33333333333333333333333333333333)です。
于 2012-09-28T14:43:17.590 に答える