Метод Якоби для нахождения собственных значений

pair<int, int> max(double **arr, int n) {
    pair<int, int> ij; 
    double res = 0;
    for (int i = 0; i < n; i++) {
        for (int j = i+1; j < n; j++)
            if (fabs(arr[i][j]) > fabs(res)) {
                res = arr[i][j];
                ij.first = i;
                ij.second = j;
            }
    }
    return ij;
}

int sign(double a) {
    if (a > 0) return 1;
    if (a < 0) return -1;
    return 0;
}

void Jacobi(double **arr, int n, int i0, int j0) {
    double cos = 0, sin = 0;
 
    if (arr[i0][i0] == arr[j0][j0])
        cos = sin = sqrt(2) / 2;
    else {
        double x = -2 * arr[i0][j0], y = arr[i0][i0] - arr[j0][j0];
        cos = sqrt(0.5 * (1 + (fabs(y) / sqrt(x*x + y*y))));
        sin = (sign(x*y)*fabs(x)) / (2 * cos*sqrt(x*x + y*y));
    }
 
    for (int m = 0; m < n; m++) {
        if (m != i0 && m != j0) {
            double ai0m, aj0m;
            if (m > i0) ai0m = arr[i0][m];
            else ai0m = arr[m][i0];
            if (m > j0) aj0m = arr[j0][m];
            else aj0m = arr[m][j0];
 
            if (m > i0) arr[i0][m] = cos*ai0m - sin*aj0m;
            else arr[m][i0] = cos*ai0m - sin*aj0m;
            if (m > j0) arr[j0][m] = sin*ai0m + cos*aj0m;
            else arr[m][j0] = sin*ai0m + cos*aj0m;
        }
    }
 
    arr[i0][i0] = cos*cos*arr[i0][i0] - 2 * sin*cos*arr[i0][j0] + sin*sin*arr[j0][j0];
    arr[j0][j0] = sin*sin*arr[i0][i0] + 2 * sin*cos*arr[i0][j0] + cos*cos*arr[j0][j0];
    arr[i0][j0] = (cos*cos - sin*sin)*arr[i0][j0] + cos*sin*(arr[i0][i0] - arr[j0][j0]); // on Wiki there is a mistake
}

void algorithm(double **arr, int n) {
    pair<int, int> ij = max(arr, n);
    while (fabs(arr[ij.first][ij.second]) > eps){
        ij = max(arr, n);
        Jacobi(arr, n, ij.first, ij.second);
        printArray(arr, n);
    };
}