Найти решение неравенства

#include <iostream>
 
using namespace std;
int main()
{
    setlocale(LC_ALL, "Russian");
    double a,b;
    cout<<"введите a и b"<<endl;
    cin>>a>>b;
    if(b>0)
        cout<<"решение не существует"<<endl;
    else
        cout<<-sqrt(-b)<<"<"<<"x"<<"<"<<sqrt(-b)<<endl;
    system("pause");
    return 0;
}

///////////////////////////////////////////////////////////////////////////////
//1.
///////////////////////////////////////////////////////////////////////////////
//Для произвольных a, b вычислить решение неравенства (x^2 + b)(ax - 1) <= 0.
///////////////////////////////////////////////////////////////////////////////
#include <algorithm>
#include <cmath>
#include <cstdlib>
#include <ctime>
#include <iostream>
#include <vector>
///////////////////////////////////////////////////////////////////////////////
typedef std::vector     < double    >   T_roots;
///////////////////////////////////////////////////////////////////////////////
double  get_a_root( double  a )
{
    return  1 / a;
}
///////////////////////////////////////////////////////////////////////////////
double  get_b_root( double  b )
{
    return  std::sqrt(-b);
}
///////////////////////////////////////////////////////////////////////////////
void    print_solution_for_three_roots
    (
        double  a,
        double  a_root,
        double  b_root
    )
{
    T_roots     roots   {
                            a_root,
                            -b_root,
                            b_root
                        };
 
    std::sort   (
                    roots.begin     (),
                    roots.end       ()
                );
 
    if( a > 0 )
    {
        std::cout   <<  "x <= "
                    <<  roots[0]
                    <<  " or x >= "
                    <<  roots[1]
                    <<  " and x <= "
                    <<  roots[2]
                    <<  std::endl;
    }
    else
    {
        std::cout   <<  "x >= "
                    <<  roots[0]
                    <<  " and x <= "
                    <<  roots[1]
                    <<  " or x >= "
                    <<  roots[2]
                    <<  std::endl;
    }//else
}
///////////////////////////////////////////////////////////////////////////////
void    print_solution_for_one_root
    (
        double  a,
        double  a_root
    )
{
    if( a > 0 )
    {
        std::cout   <<  "x <= "
                    <<  a_root
                    <<  std::endl;
    }
    else
    {
        std::cout   <<  "x >= "
                    <<  a_root
                    <<  std::endl;
    }//else
}
///////////////////////////////////////////////////////////////////////////////
void    print_solution_for_two_roots( double    b_root )
{
    std::cout   <<  "x <= "
                <<  -b_root
                <<  " or x >= "
                <<  b_root
                <<  std::endl;
}
///////////////////////////////////////////////////////////////////////////////
void    print_solution
    (
        double  a,
        double  b
    )
{
    double  a_root  =   a   ==  0
                            ?   0
                            :   get_a_root(a);
 
    double  b_root  =   b   <=  0
                            ?   get_b_root(b)
                            :   0;
 
    std::cout   <<  std::fixed;
 
    if( a != 0 )
    {
        std::cout   <<  "a_root = "     <<  a_root  <<  std::endl;
 
        if( b <= 0 )
        {
            std::cout   <<  "b_root = "     <<  b_root  <<  std::endl;
 
            print_solution_for_three_roots
                (
                    a,
                    a_root,
                    b_root
                );
        }
        else
        {
            print_solution_for_one_root
                (
                    a,
                    a_root
                );
        }//else
    }
    else
    {
        if( b < 0 )
        {
            std::cout   <<  "b_root = "     <<  b_root  <<  std::endl;
 
            print_solution_for_two_roots( b_root );
        }
        else
        {
            std::cout   <<  "x any number"
                        <<  std::endl;
        }
    }//else
}
///////////////////////////////////////////////////////////////////////////////
void    rand_set_a_and_b
    (
        double  &   a,
        double  &   b
    )
{
    double  a_root{};
 
    if  (
            rand() % 2
        )
    {
        b   =   -rand();
 
        switch( rand() % 3 )
        {
        case    0   :
            a   =   0;
            break;
 
        case    1   :
            a_root  =   get_b_root(b) * 2;
            a       =   1 / a_root;
 
        default :
            a_root  =   get_b_root(b) / 2;
            a   =   1 / a_root;
            break;
        }//switch
    }
    else
    {
        b   =   rand();
 
        a   =   rand() % 2
                    ?   0
                    :   rand();
    }//else
 
    if( rand() % 2 )
    {
        a   *=  -1;
    }
}
///////////////////////////////////////////////////////////////////////////////
int     main()
{
    srand(unsigned(time(0)));
 
    for(;;)
    {
        double  a{};
        double  b{};
 
        rand_set_a_and_b( a, b );
 
        std::cout   <<  "\n\n"
                    <<  "a = "  <<  a   <<  std::endl
                    <<  "b = "  <<  b   <<  std::endl
                    <<   "(x^2 + b)(ax - 1) <= 0 for x values:"
                    <<  std::endl;
 
        print_solution(a, b);
        system("pause");
    }//for
}