Код решения кубического уравнения

#pragma once
#include <iostream>
#include <stdio.h>  
#include <conio.h>
#include <locale.h>
#include <math.h>
const double PI = 3.1415926535897932384626433832795;
using namespace std;
double a,b,c,d,p,q,s,f,x1,x2,x3,x2i,x3i;
int main()
{
setlocale(LC_ALL,"Russian");
m:cout<<"Введите коэффициенты кубического уравнения Ax^3+Bx^2+Cx+D"<<endl;
cout<<"A= "<<endl;
cin>>a;
cout<<"B= "<<endl;
cin>>b;
cout<<"C= "<<endl;
cin>>c;
cout<<"D= "<<endl;
cin>>d;
if (a==0)
{
    cout<<"Ошибка !"<<endl;
    goto m;
}
else
{
p=((3.*a*c-b*b)/(3.*a*a));
q=((2.*b*b-9.*a*b*c+27.*a*a*d));
s=(((q*q)/4.)+(p*p*p)/27.);
if (q<0)
{
    f=(atan(pow(-s,0.5)/(-q/2)));
}
else if (q>0)
{
    f=(atan(pow(-s,0.5)/(-q/2))+PI);
}
else
{
    f=(PI/2);
}
if (s<0)
{
    x1=(2.*pow((-p/3.),0.5)*cos(f/3.)-b/3.*a);
    x2=(2.*pow((-p/3.),0.5)*cos(f/3.+(2.*PI)/3.)-b/3.*a);
    x3=(2.*pow((-p/3.),0.5)*cos(f/3.+(2.*PI)/3.)-b/3.*a);
}
else if (s>0)
{
x1=(pow(-q/2.+pow(s,0.5),1./3.)+pow(-q/2.-(pow(s,0.5)),1./3.)-b/(3.*a));
    x2=(-0.5*(pow(-q/2.+pow(s,0.5),1./3.)+pow(-q/2.-(pow(s,0.5)),1./3.)-b/(3.*a)));
    x2i=((pow(3.,0.5)/2.)*(pow(-q/2.+(pow(s,0.5)),1./3.)-pow(-q/2.-(pow(s,0.5)),1./3.)));
    x3=(-0.5*(pow(-q/2+pow(s,0.5),1/3)+pow(-q/2-(pow(s,0.5)),1/3)-b/(3*a)));
    x3i=((pow(3.,0.5)/2.)*(pow(-q/2.+(pow(s,0.5)),1./3.)-pow(-q/2.-(pow(s,0.5)),1./3.)));
}
else
{
    x1=(2.*pow(-q/2.,1./3.)-b/(3.*a));
    x2=(-1.*pow(-q/2.,1./3.)-b/(3.*a));
    x3=(-1.*pow(-q/2.,1./3.)-b/(3.*a));
}
cout<<"Х1= "<<x1<<endl;
if (s>0)
{
    cout<<"X2= "<<x2<<"+"<<x2i<<"i"<<endl;
    cout<<"X3= "<<x3<<"-"<<x3i<<"i"<<endl;
}
else
{
cout<<"X2= "<<x2<<endl;
cout<<"X3= "<<x3<<endl;
}
}
system("pause");
}