Код аппроксимации данных линейной функции

    double[][] matrix = {
        {1, 2, 3, 4, 5},//x
        {5.1, 4.4, 3.2, 2.7, 2.55},};//y
 
    double[][] matrix1;
    public int XX = 5, YY = 2;
    public int l = 15, k, m, j;
    double d1 = 0;
 
    double[] data = new double[15];
 
    double[] a = new double[XX];
    double[] b = new double[XX];
    double[] d = new double[XX];
 
    public static void main(String[] args) {
        МНК Zadacha5 = new МНК();
        Zadacha5.Initialization();
        Zadacha5.Start();
        Zadacha5.Start2();
        Zadacha5.Start3();
 
    }
 
    void Initialization() {
        matrix1 = new double[YY][XX];
 
    }
 
    public void Start() {
        double x1, y1;
        for (j = 0; j < XX; j++) {
            switch (j) {
                case 0: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[0][i] = matrix[0][i];
                        matrix1[1][i] = matrix[1][i];
                    }
                    break;
                }
                case 1: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[1][i] = matrix[0][i] * matrix[1][i];
                    }
                    break;
                }
                case 2: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[1][i] = 1 / matrix[1][i];
                    }
                    break;
                }
                case 3: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[1][i] = matrix[0][i] / matrix[1][i];
                    }
                    break;
                }
                case 4: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[1][i] = Math.log10(matrix[1][i]);
                    }
                    break;
                }
                case 5: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[0][i] = Math.log10(matrix[0][i]);
                        matrix1[1][i] = matrix[1][i];
                    }
                    break;
                }
                case 6: {
                    for (int i = 0; i < XX; i++) {
                        matrix1[1][i] = Math.log10(matrix[1][i]);
                    }
                    break;
                }
            }
            paramenties();
 
            d1 = d[0];
            k = 0;
            for (int i = 0; i < XX; i++) {
                if (d[i] < d1) {
                    d1 = d[i];
                    k = i;
                }
            }
            
            DecimalFormat dF = new DecimalFormat("#.###");
            System.out.println();
            for (int i = 0; i < YY; i++) {
                for (int j = 0; j < XX; j++) {
                    System.out.print(" " + dF.format(matrix1[i][j]));
                }
                System.out.println();
            }
            System.out.println();
 
 
            System.out.println("Формула N " + k + "\n d = " + dF.format(d1) + 
                    "\n a = " + dF.format(a[k]) + "\n b = " + dF.format(b[k]) + "\n");
        }
    }
 
    public void paramenties() {
        double k0 = 0, b0 = 0, a1 = 0, b1 = 0, a2 = 0, b2 = 0, f2 = 0;
        for (int i = 0; i < XX; i++) {
            a1 = a1 + matrix1[1][i];
            b1 = b1 + matrix1[0][i];
            a2 = a2 + matrix1[1][i] * matrix1[1][i];
            b2 = b2 + matrix1[0][i] * matrix1[0][i];
        }
        d1 = XX * a2 - a1 * a1;
        k0 = (XX * b2 - a1 * b1) / d1;
        b0 = (b1 * a2 - a1 * b2) / d1;
 
        for (int i = 0; i < XX; i++) {
            f2 = f2 + matrix1[0][i] * matrix1[0][i];
            d1 = d1 + (matrix1[0][i] - k0 * matrix1[1][i] - b0)
                    * (matrix1[0][i] - k0 * matrix1[1][i] - b0);
        }
        d[j] = Math.sqrt(d1 / f2);
        a[j] = k0;
        b[j] = b0;
    }
 
    public void Start2() {
        double mx = 0;//summ X
        double my = 0;//summ Y
        double a2 = 0;
        double a11 = 0;
        double kk = 0, bb = 0;
        for (int i = 0; i < XX; i++) {
            mx = mx + matrix[0][i];
            my = my + matrix[1][i];
        }
 
        mx = mx / XX;
        my = my / XX;
 
        for (int i = 0; i < XX; i++) {
            for (int j = 0; j < XX; j++) {
                a2 = a2 + matrix[0][i] * matrix[0][j];//vector X^T * X
            }
        }
 
        for (int i = 0; i < XX; i++) {
            for (int j = 0; j < XX; j++) {
                a11 = a11 + matrix[0][i] * matrix[1][j];//vector X^T * Y
            }
        }
        a2 = a2 / XX;
        a11 = a11 / XX;
        kk = (a11 - mx * my) / (a2 - Math.pow(mx, 2));
        bb = my - kk * mx;
        System.out.println(" k= " + kk + " b= " + bb);
    }
 
    public void Start3() {
        double mx = 0;//summ X
        double my = 0;//summ Y
        double sumx2 = 0;
        double sumx_y = 0;
        for (int i = 0; i < XX; i++) {
            mx = mx + matrix[0][i];
            my = my + matrix[1][i];
            sumx2 = sumx2 + Math.pow(matrix[0][i], 2);
            sumx_y = sumx_y + matrix[0][i] * matrix[1][i];
        }
        double determinant = sumx2 * XX - mx * mx;
        if (determinant == 0) {
            System.out.println("Определитель равен нулю. "
                    + "Обратной матрицы не сущетсвует.");
            System.exit(0);
        }
 
        double[][] obratnay = new double[2][2];
        obratnay[0][0] = XX * (1 / determinant);
        obratnay[0][1] = -1 * mx * (1 / determinant);
        obratnay[1][0] = -1 * mx * (1 / determinant);
        obratnay[1][1] = sumx2 * (1 / determinant);
 
        double b = obratnay[0][0] * sumx_y + obratnay[0][1] * my;
        double a = obratnay[1][0] * sumx_y + obratnay[1][1] * my;
 
        System.out.println(" k= " + a + " b= " + b);
    }
 
}