Динамический массив, квадратная матрица

int n, i ,j;
cout<<"Enter size : ";
cin>>n;
int **a=new int *[n];
for (i=0; i<n; i++)
a[i]=new int[n];
cout<<"Enter matrix\n";
for (i=0; i<n; i++)
     for (j=0; j<n; j++)
          cin>>a[i][j];

#include <iostream>
#include <iomanip>
 
using namespace std;
 
 
typedef int**ppint;
 
 
void Mem_matr(ppint &a, int n)
{
    int i;
  int **a=new int *[n];
for (i=0; i<n; i++)
a[i]=new int[n];
 
}
 
        void In_mas(ppint &a, int n) // Ввод массива. Первый индекс является необязательным т.к.  функции передаётся ссылка на первый элеемнт массива(в данном случа сам массив)
        {
 
            for(int i=0;i<n;i++)
            {
                
                for(int j=0;j<n;j++)
                {
                    cout <<"\n a["<<i<<"]["<<j<<"]=";
                    cin>>a[i][j];
                }
            }
         }
 
                void Out_mas(ppint &a, int n) //Вывод массива. последующие индексы нужны для того чтобы компилятор мог корректно выделить место под элементы массива
                {
                    for (int i=0;i<n;i++)
                    {
                        for (int j=0;j<n;j++)
                        {
 
                            cout<< a[i][j]<< " ";
                        }
                        cout << endl;
                    }
                }
 
                 void field(ppint&b, ppint &a, int n) //функция вычисления закрашенной области
                 {
                        for (int i=0; i<n; i++)
                        {
                            for (int j=0; j<n; j++)
                            {
                                b[i][j] = 0;
 
                                for (int x=i; x<n; x++) // for (x=n/2;x<n;x++)
                                {
                                    for (int y=0; y<=j; y++) // for(y=0;y<n/2;y++)
                                    {
                                        b[i][j]=b[i][j]+a[x][y];
                                    }
                                }
                            }
                        }
                 }
 
                void Out_mas2 (ppint&b, int n)
                {
                for (int i=0;i<n;i++)
                {
                    for (int j=0;j<n;j++)
                    {
                        cout<<setw(3)<< b[i][j]<< " ";
                    }
                    cout << endl;
                }
                cout << endl << endl;
                }
 
                void Del(ppint&a,ppint&b, int n)
                {
                    for (int i = 0; i < n; i++)
                        delete []a[i];
                    delete []a;
                    delete []b;
                    a=NULL;
                    b=NULL;
                }
                int main()
                {
                    setlocale(LC_ALL, "Russian");
                    int n;
                    ppint a;
                    ppint b;
                    char c;
 
                    do{
 
                  cout << "Введите размерность матрицы:" << endl;
                 cin>>n;
                
                  Mem_matr(a,n);
                   cout << "Введите ваши числа:" << endl;
                   In_mas(a,n);
                   cout << "Полученый массив : " << endl;
                   Out_mas(a,n);
 
                   field(b,a,n);
 
                   cout << "Результат" << endl;
                   Out_mas2(b,n);
                   Del(a,b,n);
 
                  cout << "Хотите продолжить? Нажмите y если да." << endl;
                  cin >> c;
                 }while(c=='y');
 
              }

a=new int*[n]