Аварийное завершение программы (Process terminated with status -1073741510). Что делать?

#include "math.h"
#include "iostream"
#include "algorithm"
#include <iomanip>
#include <vector>
using namespace std;
float m,r,y;
int i,j;
 
 
void VvodMas(float **A, int const N)
{
  cout << "Исходный массив: "<< endl;
    for (i=1; i<= N; i++){
        for (j=1; j<= N; j++){
            A[i][j]=i/j-log(i*j);
            cout <<  setw(9) << A[i][j] << "  ";
        }
     cout << endl;
     }
 
}
 
void Vek(float **A, int const N, std::vector<float> &X)
{
  cout << "Исходный вектор: "<< endl;
    for (i=1; i<= N; i++){
        m= A[i][1]*fabs(A[1][i]);
        for (j=2; j<= N; j++){
            max(m,A[i][j]*fabs(A[j][i]));
        }
        X.push_back(m);
        cout << setw(9) << "X["<< i <<"]=" << X[i] ;
    }
   cout << endl;
}
 
void SortMas(float **A, int const N)
{
  cout << "Отсортированный массив: "<< endl;
    for (i=1; i<= N-1; i++){
        for (j=1; j<= N-1; j++){
            if (A[j][j]> A[j+1][j+1]){
                         r=A[j][j];
                         A[j][j]= A[j+1][j+1];
                         A[j+1][j+1]= r;
                         }
        }
    }
    for (i=1; i<= N; i++){
        for (j=1; j<= N; j++){
            cout << setw(9) << A[i][j] << "  " ;
        }
        cout << endl;
    }
}
void ProizMas(std::vector<float> X, int const N)
{
    y=1;
    cout << "Произведение массива: "<< endl;
    for (i=1; i<= N-1; i++){
        for (j=1; j<= N-1; j++){
            y=y*X[j]/(j*i);
        }
    }
    cout << "y=" << y << endl;
 
}
int main()
{
setlocale(LC_CTYPE, "rus");
cout << " Введите длинну квадратной матрицы:\n";
int N;
cin >> N;
float **A = new float *[N];
for (int i=1; i<=N; i++)
A[i] = new float [N];
 
std::vector<float> X{};
VvodMas(A,N);
Vek(A,N, X);
SortMas(A,N);
ProizMas(X,N);
 
}

#include "math.h"
#include "iostream"
#include "algorithm"
#include <iomanip>
#include <vector>
using namespace std;
float m,r,y;
int i,j;
 
 
void VvodMas(float **A, int const N)
{
  cout << "Исходный массив: "<< endl;
    for (i=0; i< N-1; i++){
        for (j=0; j< N-1; j++){
            A[i][j]=i/j-log(i*j);
            cout <<  setw(9) << A[i][j] << "  ";
        }
     cout << endl;
     }
 
}
 
void Vek(float **A, int const N, std::vector<float> &X)
{
  cout << "Исходный вектор: "<< endl;
    for (i=0; i< N-1; i++){
        m= A[i][1]*fabs(A[1][i]);
        for (j=0; j< N-1; j++){
            max(m,A[i][j]*fabs(A[j][i]));
        }
        X.push_back(m);
        cout << setw(9) << "X["<< i <<"]=" << X[i] ;
    }
   cout << endl;
}
 
void SortMas(float **A, int const N)
{
  cout << "Отсортированный массив: "<< endl;
    for (i=0; i< N-1; i++){
        for (j=0; j< N-1; j++){
            if (A[j][j]> A[j+1][j+1]){
                         r=A[j][j];
                         A[j][j]= A[j+1][j+1];
                         A[j+1][j+1]= r;
                         }
        }
    }
    for (i=0; i< N-1; i++){
        for (j=0; j< N-1; j++){
            cout << setw(9) << A[i][j] << "  " ;
        }
        cout << endl;
    }
}
void ProizMas(std::vector<float> X, int const N)
{
    y=1;
    cout << "Произведение массива: "<< endl;
    for (i=0; i< N-1; i++){
        for (j=0; j< N-1; j++){
            y=y*X[j]/(j*i);
        }
    }
    cout << "y=" << y << endl;
 
}
int main()
{
setlocale(LC_CTYPE, "rus");
cout << " Введите длинну квадратной матрицы:\n";
int N;
cin >> N;
float **A = new float *[N];
for (int i=0; i<N-1; i++)
A[i] = new float [N];
 
std::vector<float> X{};
VvodMas(A,N);
Vek(A,N, X);
SortMas(A,N);
ProizMas(X,N);
 
}

A[i][j]=(float)(i+1)/(j+1)-log((i+1)*(j+1));

#include <iostream>
#include <iomanip>
#include <cmath>
using namespace std;
 
void Set(double **&A, int n)
{
    for (int i = 0; i < n; i++)
        for (int j = 0; j < n; j++)
        A[i][j]=(double)(i+1)/(j+1)-log((i+1)*(j+1));
}
 
double *CreateX(double **&A, int n)
{    
    double *X = new double[n];
    
    double m;
    for (int i = 0; i < n; i++)
    {
        X[i]=A[i][0]*fabs(A[0][i]);
        for (int j = 1; j < n; j++){
            m=A[i][j]*fabs(A[j][i]);
        if(m>X[i]) X[i]=m;
        }
    }
    return X;
}
 
void SortMainDiagonal(double **&A, int n)
{    
    for (int i = 0; i < n-1; i++)    
        for (int j = 0; j < n-i-1; j++)
            if (A[j][j]> A[j+1][j+1])
            swap(A[j][j],A[j+1][j+1]);  
}
 
void Print2D (double **&A, int n)
{
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)       
        cout << setw(10) << setprecision(6) << fixed << A[i][j] << " ";        
    cout << "\n";
    }
}
 
void Print (double *&X, int n)
{
    for (int i = 0; i < n; i++)             
        cout << setw(10) << setprecision(6) << fixed << X[i] << " ";        
    cout << "\n";    
}
 
void Free2D (double **&A, int n)
{
    for (int i = 0; i < n; i++)
    delete[]A[i];
    delete[]A;
}
 
void Free(double *&X)
{
    delete[]X;
} 
 
int main()
{
    int n;
    cout << "n="; cin >> n;
    
      double **A = new double*[n];
    for (int i = 0; i < n; i++)
      A[i] = new double[n];
    
    Set(A,n);
    cout << "Matrix A:\n";
    Print2D(A,n);
    
    double *X=CreateX(A,n);
    
    cout << "Vector X:\n";
    Print(X,n);
    SortMainDiagonal(A,n);
    cout << "Modified matrix A:\n";
    Print2D(A,n);
    Free2D(A,n);
    Free(X);
system("pause");
return 0;
}