#include <iostream>
#include "quadratic_roots.h"

void print_pair( std::pair<double, double> p ) {
    std::cout << "(" << p.first << ", " << p.second << ")" << std::endl << std::endl;
}

int main() {
    std::cout.precision( 17 );

    std::cout << "(-0.77697504995388685, 0.55958374560606076)" << std::endl;
    print_pair( quadratic_roots( 2.3, 0.5, -1.0 ) );

    std::cout << "-1e100, -1.0" << std::endl;
    print_pair( quadratic_roots( -1e400, -1e400, -1e300 ) );

    std::cout << "-1e150, -1e-150" << std::endl;
    print_pair( quadratic_roots( 1e-75,  1e75,  1e-75 ) );

    std::cout << "1e-150, 1e150" << std::endl;
    print_pair( quadratic_roots( 1e-75, -1e75,  1e-75 ) );

    std::cout << "-1e-150, 1e150" << std::endl;
    print_pair( quadratic_roots( 1e-75, -1e75, -1e-75 ) );

    std::cout << "-1e150, 1e-150" << std::endl;
    print_pair( quadratic_roots( 1e-75,  1e75, -1e-75 ) );

    return 0;
}