Równania

// Równanie kwadratowe
// (C)2019 mgr Jerzy Wałaszek
// Metody numeryczne
//---------------------------

#include <iostream>
#include <cmath>
#include <iomanip>

using namespace std;

// Tutaj definiujemy dane wejściowe
//----------------------------------
double eps = 1e-12;  // Dokładność porównania z zerem

// Funkcja rozwiązująca równanie kwadratowe
// a,b,c - współczynniki wielomianu
//-----------------------------------------
void qe(double a, double b, double c)
{
    double delta,x1,x2,xre,xim;

    cout << endl;
    cout << "Równanie: " << a << " * x^2 + " << b << " * x + " << c << " = 0" << endl;

    // Jeśli a = 0, to mamy równanie liniowe

    if(fabs(a) < eps)
    {
        cout << "Równanie redukuje się do równania liniowego." << endl;
        if(fabs(b) <= eps)
        {
            cout << "Równanie jest ";
            if(fabs(c) <= eps) cout << " nieoznaczone";
            else               cout << " sprzeczne";
        }
        else
        {
            x1 = - c / b;
            cout << "Wynik x = " << x1;
        }
    }
    else
    {
        // Obliczamy wyróżnik

        delta = b * b - 4 * a * c;

        // W zależności od wartości wyróżnika mamy 3 różne przypadki

        if(fabs(delta) < eps) // 1. delta = 0
        {
            x1 = - b / 2 / a;
            cout << "PIERWIASTEK PODWÓJNY:" << endl
                 << "x1 = x2 = " << x1;
        }
        else if (delta < 0)   // 2. delta < 0
        {
            xre = - b / 2 / a;          // Część rzeczywista
            xim = sqrt(-delta) / 2 / a; // Część urojona
            cout << "PIERWIASTKI ZESPOLONE:" << endl
                 << "x1 = " << xre << " + " << xim << "i" << endl
                 << "x2 = " << xre << " + " << -xim << "i";
        }
        else                  // 3. delta > 0
        {
            delta = sqrt(delta);
            x1 = (- b - delta) / 2 / a;
            x2 = (- b + delta) / 2 / a;
            cout << "PIERWIASTKI RZECZYWISTE:" << endl
                 << "x1 = " << x1 << endl
                 << "x2 = " << x2;
        }
    }
    cout << endl;
}

// Program główny
//---------------

int main()
{
    setlocale(LC_ALL,"");

    cout << fixed;
    cout << "Równania kwadratowe" << endl
         << "-------------------" << endl;

    qe(3,12,15);
    qe(9,12,4);
    qe(1,-1,-2);

    cout << endl;

    return 0;
}