Вызвано исключение по адресу 0x013C1891 в Jacobi.exe: 0xC0000005: нарушение прав доступа при чтении по адресу

#include <math.h>
const double eps = 0.001; ///< желаемая точность 
 
/// N - размерность матрицы; A[N][N] - матрица коэффициентов, F[N] - столбец свободных членов,
/// X[N] - начальное приближение, ответ записывается также в X[N];
void Jacobi(int N, double** A, double* F, double* X)
{
    double* TempX = new double[N];
    double norm; // норма, определяемая как наибольшая разность компонент столбца иксов соседних итераций.
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
    delete[] TempX;
}
 
int main()
{
    const int N = 3;
    double A[N][N] = { { 5, -4, -1 }, { 4, 1, -2 }, { 3, 1, -5 } };
    double F[N] = { -2, 8, 10 };
    double X[N] = { 0 };
    Jacobi(N, (double**)A, F, X);
}

int main()
{
    const int N = 3;
    double B[N][N] = { { 5, -4, -1 }, { 4, 1, -2 }, { 3, 1, -5 } };
    double F[N] = { -2, 8, 10 };
    double X[N] = { 0 };
    double **A=new double*[N];
    for (int i=0;i<N;++i) {
        A[i]=new double[N];
        for (int j=0;j<N;++j)
            A[i][j]=B[i][j];
    }
    Jacobi(N, A, F, X); //B так передавать нельзя, а A - можно 
    return 0;
}

void Jacobi(int N, double A[][3], double* F, double* X)

#include <math.h>
#include <vector>
 
void Jacobi(std::vector<std::vector<double>> A, std::vector<double> F, std::vector<double> X, double eps)
{
    size_t N = A.size();
    std::vector<double> TempX(N, 0);
    double norm(0.0);
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
}
 
int main()
{
    std::vector<std::vector<double>> A{ { 5, -4, -1 },{ 4, 1, -2 },{ 3, 1, -5 } };
    std::vector<double> F{ -2, 8, 10 };
    std::vector<double> X(3, 0);
    Jacobi(A, F, X, 0.001);
}

#include <math.h>
#include <vector>
 
void Jacobi(const std::vector<std::vector<double>> &A, const std::vector<double> &F, std::vector<double> &X, double eps)
{
    size_t N = A.size();
    std::vector<double> TempX(N, 0);
    double norm(0.0);
 
    do {
        for (int i = 0; i < N; i++) {
            TempX[i] = F[i];
            for (int g = 0; g < N; g++) {
                if (i != g)
                    TempX[i] -= A[i][g] * X[g];
            }
            TempX[i] /= A[i][i];
        }
        norm = fabs(X[0] - TempX[0]);
        for (int h = 0; h < N; h++) {
            if (fabs(X[h] - TempX[h]) > norm)
                norm = fabs(X[h] - TempX[h]);
            X[h] = TempX[h];
        }
    } while (norm > eps);
}
 
int main()
{
    std::vector<std::vector<double>> A{ { 5, -4, -1 },{ 4, 1, -2 },{ 3, 1, -5 } };
    std::vector<double> F{ -2, 8, 10 };
    std::vector<double> X(3, 0);
    Jacobi(A, F, X, 0.001);
}