СЛАУ метод Зейделя

public class Zeidel extends IOStream
{
    final private int n = 3;
    final private double[][] koef = new double[n][n];
    private double[] otvet = new double[n];
    private double eps = 0.0001; 
    void InputKoef()
    {
        PrintLine("Введите матрицу коэфициентов:");
        int i, j;
        for(i = 0; i < n; ++i)
            for(j = 0; j < n; ++j)
                koef[i][j] = ReadDouble();
        PrintLine("Введите ответы:");          
        for(i = 0; i < n; ++i)
            otvet[i] = ReadDouble();    
    }
    void Out()
    {
        int i, j;
        PrintLine("Система уравнений:");
        for(i = 0; i < n; ++i, PrintLine())
            for(j = 0; j < n; ++j)
                Print(koef[i][j] + "x" + (j + 1) + 
                      (j == n - 1 ? " = " + otvet[i] : " + "));
    }
    private double[] x = new double[n];
    private boolean Equal(double x0, double x1, double x2)
    {
        return
            (((koef[0][0] * x0 + koef[0][1] * x1 + koef[0][2] * x2) - otvet[0] <= eps)&&
            ((koef[1][0] * x0 + koef[1][1] * x1 + koef[1][2] * x2) - otvet[1] <= eps) &&
            ((koef[2][0] * x0 + koef[2][1] * x1 + koef[2][2] * x2) - otvet[2] <= eps));
    }
    
    double[] z = new double[n];
 
    private static double Round(double chislo, int sign, int fract)
    {
        return (double)(int)(chislo * Math.pow(10, sign + fract)) / Math.pow(10, fract);
    }
    void Calculate()
    {
        int i, j;
        for(i = 0; i < n; ++i)
        {
            otvet[i] /= koef[i][i];
            z[i] = koef[i][i];
            for(j = 0; j < n; ++j)
                koef[i][j] /= z[i];
        }
        x[1] = x[2] = 0;
        i ^= i;
        do
        {
            x[0] = otvet[0] - koef[0][1] * x[1] - koef[0][2] * x[2];
            x[1] = otvet[1] - koef[1][0] * x[0] - koef[1][2] * x[2];
            x[2] = otvet[2] - koef[2][0] * x[0] - koef[2][1] * x[1];
            PrintLine("\ni = " + ++i +"\tx1 = " + Round(x[0], 4, 5) + "\tx2 = " + Round(x[1], 4, 5) + "\tx3 = " + Round(x[2], 4, 5));
        }while(!Equal(x[0], x[1], x[2]));
        PrintLine("\nx1 = " + Round(x[0], 4, 5) + "\tx2 = " + Round(x[1], 4, 5) + "\tx3 = " + Round(x[2], 4, 5));
    }
    void Verification()
    {
        int i, j;
        PrintLine("\n\nСистема уравнений:");
        for(i = 0; i < n; ++i, PrintLine())
            for(j = 0; j < n; ++j)
                Print(koef[i][j]*z[i] + " * " + Round(x[j], 4, 5) + 
                      (j == n - 1 ? " = " + otvet[i]*z[i] : " + "));
    }
}

public class Main extends IOStream
{
    public static void main(String[] args) 
    {
        jar = false;
        Zeidel z = new Zeidel();
        z.InputKoef();
        z.Out();
        z.Calculate();
        z.Verification();
    }
}