Некорректная работа функции pow

#include<iostream>
using namespace std;
#include<math.h>
#include<conio.h>
#include<Windows.h>
#define M_PI       3.14159265358979323846
 
void vertification(double& number);
double a, b, c, p, q, d, u, v, r, f, y, y2, y3, i1, i2, x1, x2, x3;
 
 
int main()
{
    int repeat = 1;
    do{
    cout << "x^3+ax^2+bx+c=0" << endl;
    SetConsoleOutputCP(1251);
    cout << "Введіть коефіцієнт а: ";
    vertification(a);
    cout << "Введіть коефіцієнт b: ";
    vertification(b);
    cout << "Введіть коефіцієнт с: ";
    vertification(c);
    cout << "Рівняння матиме вигляд: "<<endl;
    cout << "x^3";
    if (a < 0)
        cout << a<<"x^2";
    else
        cout << "+" << a << "x^2";
    if (b < 0)
        cout << b << "x";
    else
        cout << "+" << b << "x";
    if (c < 0)
        cout << c;
    else
        cout << "+" << c;
    cout << "=0" << endl;
 
    p = b - pow(a,2) / 3;
    q = 2 * pow(a,3) / 27 - a*b / 3 + c;
    d = pow(p, 3) / 27 + pow(q, 2) / 4;
    u = pow(-q / 2 + sqrt(d), 1.0 / 3.0);
    v = -p / (3 * u);
    r = pow(-pow(p, 3) / 27, 0.5);
    f = acos(-q / 2 * r);
 
    if (d == 0)
    {
        y = 3 * p / q - a / 3;
        y2 = -3 * q / (2 * p) - a / 3;
        y3 = y2;
        x1 = y;
        x2 = y2;
        x3 = y3;
    }
    if (d < 0)
    {
        y = 2 * fabs(pow(r, 1.0 / 3.0))*cos(f / 3) - a / 3;
        y2 = 2 * fabs(pow(r, 1.0 / 3.0))*cos((f + 2 * M_PI) / 3) - a / 3;
        y3 = 2 * fabs(pow(r, 1.0 / 3.0))*cos((f + 4 * M_PI) / 3) - a / 3;
        x1 = y;
        x2 = y2;
        x3 = y3;
    }
    if (d>0)
    {
        y = u + v - a / 3;
        y2 = -(u + v) / 2 - a / 3;
        i1 = pow(3, 2)*(u - v) / 2;
    }
    cout << "Корені рівняння:" << endl;
    if (d <= 0)
    {
        cout << x1<<endl;
        cout << x2<<endl;
        cout << x3<<endl;
    }
    if (d > 0)
    {
        cout << y-a / 3 << endl;;
        cout << y2 - a / 3 << "+" << "i*" << i1 << endl;
        cout << y2 - a / 3 << "-" << "i*" << i1 << endl;
    }
    if (d <= 0)
    {
        cout << pow(x1, 3) + a*pow(x1, 2) + b*x1 + c << endl;
    }
    cout << "Для продовження натисніть ENTER, для виходу з програми натисніть ESCAPE." << endl;
    while (repeat != 13 || repeat != 27)
    {
        repeat = _getch();
        if (repeat == 13 || repeat == 27)
            break;
        cout << endl;
    }if (repeat == 27)
        break;
    } while (repeat != 27);
    return 0;
}
 
void vertification(double& number) //определение функции vertification
{
    while ((cin >> number).fail() || cin.get() != '\n')//цикл который повторяется пока не будет присвоено значение переменной number
    {
        cin.clear();
        cin.sync();
        cout << "Помилка, введено некоректні дані. Введіть дані повторно." << endl;
    }
}

u = pow(-q / 2 + sqrt(d), 1.0 / 3.0);