Табулирование кусочно-непрерывной функции

#include <iostream>
#include <cmath>
#include <cstdlib>
 
int main()
{
    const double a =  -1.0;     //согласно условию
    const double X1 = -1.0;     //
    const double X2 = 2.0;      //
    const double deltaX = 0.5;  //
 
    double f = 0;    
    double X = X1;
    
    std::cout << "Table of function values: " << std::endl;
    while (X <= X2)
    {
        if (abs(a*X) > 1.0)
        {
            f = 1/(a*X) - 2 * pow((a*X),3) + 3*log(abs(X));    
        }
        else if (abs(a*X) < 1.0)
        {
            f = sqrt(a*a + 4*X*X)*log(abs(a)) + a / (sqrt(a*a + X));
        }
        else
        {
            f = pow(X, 3)/3 - pow(a, 3)/2;
        }
        
        std::cout << "f(" << X << ") = " << f << std::endl;
        
        X+=deltaX;
    }
    
    
    return 0;
}