Решение СЛАУ с комплексными коэффициентами

#include <iostream>
#include <fstream>
#include "complex.h"
#include "SLAU.h"
 
using namespace std;
 
int main ()
{
    int i, j, n;
    double L, C, E1, E2, *r, w, v;
    complex **a, *b, *X, *I;
    n = 3;
    ifstream f;
    f.open("/home/andrey/abc.txt");
    if (!f)
        cout << "Fail ne naiden";
    else
    {   
        f >> v >> L >> C >> E1 >> E2;
        w = 2*M_PI*v; cout << "f = " << v;
        cout << " Гц\nL = " << L << " Гн\nC = " << C*1e+6;
        cout << " РјР¤\nE1 = " << E1 << " Р’\nE2 = " << E2 << " Р’" << endl;
        a = new complex*[n];
        for (i = 0; i < n; i++)
            a[i] = new complex[n];
        b = new complex[n];
        I = new complex[n];     
        r = new double[n];
        for (i = 0; i < n; i++)
        {
            f >> r[i];
            cout << "r" << i+1 << " = " << r[i] << " РћРј" << endl;
        }
        X = new complex[n];
        X[0] = complex(r[0], w*L);
        X[1] = complex(r[1], -1/(w*C));
        X[2] = r[2];
        a[0][0] = 1; a[0][1] = -1; a[0][2] = 1; b[0] = 0;
        a[1][0] = X[0]; a[1][1] = X[1]; a[1][2] = 0; b[1] = E1+E2;
        a[2][0] = 0; a[2][1] = X[1]; a[2][2] = X[2]; b[2] = E2; 
        int otvet = SLAU(a, n, b, I);
        cout << endl;
        if (otvet)
            if (otvet == -2)
                cout << "Mnojestvo reshenij";
            else cout << "Net reshenija";
        else
        {
            for (i = 0; i < n; i++)
                cout << "I(" << i+1 << ") = " << I[i] << " Рђ\n";
            complex s;
            cout << "Проверка\n";
            for (i = 0; i < n; cout << s << endl, i++)
            {
                for (s = 0, j = 0; j < n; j++)
                    s += a[i][j]*I[j];
                s-=b[i];
            }
            
        }
    }
    return 0;
}