Вычисление корня квадратного из комплексного числа

#include<iostream>
#include<cmath>
using namespace std;
int main()
{
    setlocale(0,"");
    double x,y,phi,z,z1,z2,i,pi=3.14,X1,Y1,X2,Y2;
    cout<<"Введите x=";
    cin>>x;
    cout<<"Введите y=";
    cin>>y;
    phi=atan2(y,x);
    z=sqrt(x*x+y*y);
    cout<<"z1="<<sqrt(z)<<"*exp("<<phi/2<<"*i)""\n";
    cout<<"z2="<<sqrt(z)<<"*exp("<<phi/2<<"*i+i*"<<pi<<")\n";
    X1=z1*cos(phi);
    Y1=z1*sin(phi);
    X2=z2*cos(phi);
    Y2=z2*sin(phi);
    cout<<"X1=   "<<X1<<"Y1=   "<<Y1<<endl;
    cout<<"X2=   "<<X1<<"Y2=   "<<Y1<<endl;
    system("pause");
    return 0;
}

    cout<<"X1=   "<<X1<<"Y1=   "<<Y1<<endl;
    cout<<"X2=   "<<X1<<"Y2=   "<<Y1<<endl;

#include<iostream>
using namespace std;
int main()
{
    setlocale(0, "");
    double x, y, phi, z;
    cout << "Введите x=";    cin >> x;
    cout << "Введите y=";    cin >> y;
    phi = atan2(y, x);      z = sqrt(x*x + y*y);
    
    for (int k = 0; k <= 1; k++)
        cout << "z" << k + 1 << "=" << sqrt(z) << "*exp(" << phi / 2 << "+" << 2*k << "*pi)" << endl;
    
    system("pause");
    return 0;
}

#include "stdafx.h" // for Visual Studio
#include <iostream>
#include <cmath>  
using namespace std;
#define PI 3.14159265
 
int main()
{
    setlocale(0, "");
    double x, y, phi, r;
    cout << "Введите x=";    cin >> x;
    cout << "Введите y=";    cin >> y;
    phi = atan2(x, y);      r = sqrt(x*x + y*y);
    
    for (int k = 0; k <= 1; k++)
        cout << "z" << k + 1 << "=" << r << "*exp ^ i(" << phi / 2 << "+" << k*PI << ")\n";
    
    system("pause");
    return 0;
}

phi = atan2(y, x);

phi = atan2(y, x);// угол
z = sqrt(x*x + y*y); // модуль 
cout<<"z="<<x<<'+'<<y<<'i'<<endl; // вывод обычного представления
cout<<"modul="<<z<<" phi="<<phi<<endl; // вывод модуля и угла
 cout << "z=" <<z<< "*exp(i*" << phi <<")"<<endl; // вывод экспоненциального представления

#include <cmath>
#include <iostream>
using namespace std;
 
class CComplex
{
protected:
    bool bExp;
    double re;
    double im;
public:
    CComplex();
    double Abs ();
    double Angle();
    void ToExp();
    void ToAlg();
    void power(float p);
    void print();
 
    double &Re(){return re;}
    double &Im(){return im;}
};
 
CComplex::CComplex(){
    re = 0;
    im = 0;
    bExp = 0;
}
 
double CComplex::Abs(){
    if( bExp )
        return re;
    else
        return sqrt(re*re + im*im);
}
 
double CComplex::Angle(){
    if( bExp )
        return im;
    else
        return atan2(im, re);
}
 
void CComplex::ToExp(){
    if( !bExp ) {
        float abs = Abs();
        float phi = Angle();
        re = abs;
        im = phi;
        bExp = 1;
    }
}
 
void CComplex::ToAlg(){
    if( bExp )  {
    float abs = re;
    float phi = im;
    re = abs*cos(phi);
    im = abs*sin(phi);
    bExp = 0;
    }
}
 
void CComplex::power(float p){
    if( bExp )
    {
        re = pow(re, p);
        im = im*p;
    }
    else
    {
        ToExp();
        power(p);
        ToAlg();
    }
}
 
void CComplex::print(){
    if( bExp )
        cout<<"EXP FORM : "<<re<<"*exp(j*"<<im<<")"<<endl;
    else
        cout<<"ALG FORM : "<<re<<" + j*("<<im<<")"<<endl;
}
 
 
 
int main()
{
    CComplex val;
    val.Re() =  25;
    val.Im() = -63;
    val.print();
    cout<<"SQRT : "<<endl;
    val.power(0.5);
    val.print();
    val.ToExp();
    val.print();
    cin.get();
    return 0;
}