Ошибка: "Необработанное исключение: 0xC0000005: Нарушение прав доступа при записи"

#include <iostream>
#include <conio.h>
#include "mpi.h"
#include <ctime>
#include <time.h>
#include <math.h>
#include <fstream>
 
using namespace std;
 
#define LID 7077
#define _CRT_SECURE_NO_WARNINGS
 
int main(int argc, char* argv[])
{
    int myrank, ranksize;
    double time;
        
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &myrank);
    MPI_Comm_size(MPI_COMM_WORLD, &ranksize);
    MPI_Status Stat;
            
    if (myrank == 0) {
        int n = 0;
        cout << "Vvedite n:";
        cin >> n;
        double **mat = new double*[n];
        double **ans = new double*[n];
 
        cout << endl << "Matrix:" << endl;
        for(int i=0;i<n;i++)
        {
            mat[i] = new double[n];
            for(int j=0;j<n;j++)
            {
                mat[i][j] = rand() % 100;
                if(i==j) {
                    if(mat[i][j] == 0) {
                        mat[i][j] = 1;
                    }
                }
                cout << mat[i][j] << " ";
            }
            cout << endl;
        }
        double *c = new double[n];
        double *res = new double[n];
        
        time = MPI_Wtime();
 
        for (int i = 1; i < ranksize; i++)
        {
            MPI_Send(&n, 1, MPI_INT, i, LID*i + 47, MPI_COMM_WORLD);
        }
 
        for (int i = 0; i<n; i++) {
            for(int j=0;j<n;j++) {
                c[j] = mat[j][i];
            }
            MPI_Send(&i, 1, MPI_INT, i % (ranksize - 1) + 1, LID*(i % (ranksize - 1) + 1) + i / (ranksize - 1), MPI_COMM_WORLD);
            MPI_Send(c, n, MPI_DOUBLE, i % (ranksize - 1) + 1, LID*(i % (ranksize - 1) + 1) + i / (ranksize - 1), MPI_COMM_WORLD);
 
        }
        for (int i = 0; i<n; i++) {         
            MPI_Recv(ans[i], n, MPI_DOUBLE, i % (ranksize - 1) + 1, LID*(i % (ranksize - 1) + 1) + i / (ranksize - 1), MPI_COMM_WORLD, &Stat);
            for(int j=0;j<n;j++) {
                cout << res[j] << " ";
                ans[j][i] = res[j];
            }
        }
        cout << endl << "Obratnaja matrica:" << endl;
        for (int i = 0; i < n; i++)
        {
            for(int j=0;j<n;j++) {
                cout << ans[i][j] << " ";
            }
            cout << endl;
        }
        cout << endl;
        cout << "Ranksize: " << ranksize << endl;
 
        time = MPI_Wtime() - time;
        cout << "Time: " << time;
        double **pr = new double*[n];
        cout << endl << "Proverka:" << endl;
        for (int i=0;i<n;i++) {
            double s= 0;
            for (int j=0;j<n;j++) {
                for (int k=0;k<n;k++) {
                    s += mat[i][k]*ans[k][j];
                }
                pr[i][j] = floor(s+0.5);
            }       
        }
        for (int i=0;i<n;i++) {
            for (int j=0;j<n;j++) {
                cout << pr[i][j] << " ";
            }
            cout << endl;
        }
 
    }
 
    else {
 
        int nNew = 0, iNew = 0;
        MPI_Recv(&nNew, 1, MPI_INT, 0, LID*myrank + 47, MPI_COMM_WORLD, &Stat);
        
        double *vector = new double[nNew];
        double *res = new double[nNew];
 
        for (int i = 0; i< nNew / (ranksize - 1) + (nNew % (ranksize - 1) >= myrank ? 1 : 0); i++) {
            MPI_Recv(&iNew, 1, MPI_INT, 0, LID*myrank + i, MPI_COMM_WORLD, &Stat);
            MPI_Recv(vector, nNew, MPI_DOUBLE, 0, LID*myrank + i, MPI_COMM_WORLD, &Stat);
            for(int j=0;j<nNew;j++) {
                if(j==iNew) {
                    res[j] = 1.0/vector[j];
                    cout << res[j] << " ";
                }
                else {
                    res[j] = -vector[j]/vector[iNew];
                    cout << res[j] << " ";
                }
            }
            cout << endl;
            MPI_Send(res, nNew, MPI_DOUBLE, 0, LID*myrank + i, MPI_COMM_WORLD);
        }
    }
    
    MPI_Finalize();
    system("PAUSE");
    return 0;
}

MPI_Recv(res, n, MPI_DOUBLE, i % (ranksize - 1) + 1, LID*(i % (ranksize - 1) + 1) + i / (ranksize - 1), MPI_COMM_WORLD, &Stat);
            for(int j=0;j<n;j++) {
                cout << res[j] << " ";
                ans[j][i] = res[j];     //здесь не работает
            }