Почему не инициализируются переменные

#include <iostream>
#include <math.h>
using namespace std;
#define pi 3.14159
int main ( )
{
float a, b, c, d, r, s, t, p, q, D, ro, fi, x1, x2, x3, u, v, h, g;
cout<<" a = "; cin >>a;
cout<<" b = "; cin >>b;
cout<<" c = "; cin >>c;
cout<<" d = "; cin >>d;
r=b/a; s=c/a; t=d/a;
p=(3*s -r * r ) / 3; q=2* r* r * r /27 - r * s/3+t;
D=(p /3) * ( p /3) * ( p /3) +(q /2) * ( q /2);
if (D<0)
{
ro=sqrt ( ( float )( -p* p* p/27) );
fi= acos(-q /(2 * ro ));
fi=pi/2 - atan ( fi / sqrt (1 - fi * fi ) );
x1=2*pow( ro, ( float ) 1/3) * cos ( fi /3)- r /3;
x2=2*pow( ro, ( float ) 1/3) * cos ( fi /3+2* pi /3)- r /3;
x3=2*pow( ro, ( float ) 1/3) * cos ( fi /3+4* pi /3)- r /3;
cout<<" \n x1 = "<<x1<<" \t x2 = "<<x2;
cout<<" \t x3 = "<<x3<<" \n ";
}
else
{
if ( -q/2+sqrt (D) >0) u=pow(( - q/2+sqrt (D) ),( float )
1/3);
else
if ( -q/2+sqrt (D) <0) u=-pow( fabs( -q/2+sqrt (D) ),( float )
1/3);
else u=0;
if (-q/2 - sqrt (D) >0) v=pow(( -q/2 - sqrt (D) ),( float )
1/3);
else
if ( -q/2 - sqrt (D) <0) v=-pow( fabs( -q/2 - sqrt (D) ),(
float ) 1/3);
else v=0;
x1=u+v-r /3;
h= -(u+v)/2 - r /3;
g=(u-v) /2 -sqrt (( float ) 3);
cout<<" n x1 = "<<x1;
if (x2>=0)
{
cout<<x1<<" + "<<x2<<" i \t ";
cout<<x1<<" -"<<x2<<" i \n ";
}
else
{
    cout<<x1<<" -"<<fabs ( x2 )<<" i \t ";
    cout<<x1<<" + "<<fabs ( x2 )<<" i \n ";
    }
    }
    if (g>=0)
    {
    cout<<" \t x2 = "<<h<<" + "<<g<<" i ";
    cout<<" \t x3 = "<<h<<" -"<<g<<" i \n ";
    }
    else
    {
    cout<<" \t x2 = "<<h<<" -"<<fabs (g)<<" i ";
    cout<<" \t x2 = "<<h<<" + "<<fabs (g)<<" i ";
    }
    return 0;
    }

float D{}, x2{}, g{};

#include <iostream>
#include <cmath>
 
const float pi=3.14159;
int main ( )
{
    float a{}; std::cout<<" a = "; std::cin >>a;
    float b{}; std::cout<<" b = "; std::cin >>b;
    float c{}; std::cout<<" c = "; std::cin >>c;
    float d{}; std::cout<<" d = "; std::cin >>d;
 
    
    const float r=b/a;
    const float s=c/a;
    const float t=d/a;
    const float p=(3*s -r * r ) / 3;
    const float q=2* r* r * r /27 - r * s/3+t;
    const float D=p*p*p/3 +q*q/2;
    
    if (D<0)
    {
        const float ro=sqrt(-p* p* p/27);
        
        float fi= acos(-q /(2 * ro ));
        fi=pi/2 - atan ( fi / sqrt (1 - fi * fi ) );
        
        const float x1=2*cbrt(ro) * cos ( fi /3 + 0       )- r /3;
        const float x2=2*cbrt(ro) * cos ( fi /3 + 2* pi /3)- r /3;
        const float x3=2*cbrt(ro) * cos ( fi /3 + 4* pi /3)- r /3;
       
        std::cout
        <<" \n x1 = "<<x1
        <<" \t x2 = "<<x2
        <<" \t x3 = "<<x3
        <<" \n ";
    }
    else
    {
        float  u{};
        if      ( -q/2+sqrt (D) >0){ u=+cbrt(     (-q/2+sqrt(D) ) );}
        else if ( -q/2+sqrt (D) <0){ u=-cbrt( fabs(-q/2+sqrt(D) ) );}
        else {u=0;}
 
        float v{};
        if      (-q/2 - sqrt (D) >0){ v=+cbrt(     (-q/2 - sqrt (D) ) );}
        else if (-q/2 - sqrt (D) <0){ v=-cbrt( fabs(-q/2 - sqrt (D) ) );}
        else {v =0;}
        
        const float x1=u+v-r /3;
        const float h=-(u+v)/2 - r /3;
        const float g=(u-v) /2 -sqrt (3.);
        std::cout<<" n x1 = "<<x1;
        
        if (x2>=0)
        {
            std::cout<<x1<<" + "<<x2<<" i \t ";
            std::cout<<x1<<" - "<<x2<<" i \n ";
        }
        else
        {
            std::cout<<x1<<" - "<<fabs ( x2 )<<" i \t ";
            std::cout<<x1<<" + "<<fabs ( x2 )<<" i \n ";
        }
    }
        
    if (g>=0)
    {
        std::cout<<" \t x2 = "<<h<<" + "<<g<<" i ";
        std::cout<<" \t x3 = "<<h<<" - "<<g<<" i \n ";
    }
    else
    {
        std::cout<<" \t x2 = "<<h<<" - "<<fabs (g)<<" i ";
        std::cout<<" \t x2 = "<<h<<" + "<<fabs (g)<<" i ";
    }
    
    return 0;
}