// quadratic-equation.C --- 複素係数の 2 次方程式を解く。 // // g++ -o quadratic-equation quadratic-equation.C // // 注意: 素朴な (桁落ちの対策などしていない) アルゴリズムを使っている。 // // 2016/3/22 修正 // 2000年頃動いていたプログラムがコンパイル出来なかったので修正。 // iostream.h, complex.h をインクルードしていた → .h を取る // namespace を指定していなかった (必要なかった) → using namespace std; // 4, 2 などの int 型を使っていた → double 型にする // 最近は拡張子が大文字のCなのは珍しい気がする (.cppとかが多い?) けど。 #include <iostream> #include <complex> using namespace std; int main(void) { complex<double> a, b, c, D, x1, x2; cout << "複素係数の2次方程式 a x^2+b x+c=0 (a≠0) を解きます。" << endl; cout << " 複素数は ( ) でくくり、実部と虚部をカンマ , で区切って表す。" << endl; cout << " 例えば 1+2i は (1,2) と表わします。" << endl; cout << "入力してください。" << endl; cout << "a="; cin >> a; cout << "b="; cin >> b; cout << "c="; cin >> c; cout << "a=" << a << ", b=" << b << ", c=" << c << endl; D = sqrt(b * b - 4.0 * a * c); x1 = (-b + D) / (2.0 * a); x2 = (-b - D) / (2.0 * a); cout << "x1=" << x1 << endl; cout << "x2=" << x2 << endl; return 0; }