double SofPhi(const int& x, const complex<double>& phix);
double rhoPhi(const int& x, const complex<double>& phixCandidate);
};
double SofPhi(const int& x, const complex<double>& phix);
double rhoPhi(const int& x, const complex<double>& phixCandidate);
};