Calcul matriciel en c++: Somme de la première diagonale

using namespace std;
#include<stdio.h>
#include<stdlib.h>
#include <iostream>
#include <time.h>

int somme_diagonale(int **M,int n){
    int somme=0;
    for(int i = 0 ; i < n ; i++){
              somme += M[i][i];
    }
    return somme;  
}
void affichage(int** Mat, int N, int M){
     for(int i = 0; i < N; i++) {
        for(int j = 0; j < M; j++) {
            cout << Mat[i][j] ;
            cout << "\t";
        }
        cout << "\n";
     }
     cout << "\n";
}
main()
{
      int n=3;
      int** M  = new int* [ n ];
      for (int i=0; i < n; i++)
      M[i] = new int[ n ];
      srand(time(NULL));
      for(int i=0; i<n; i++)
      for(int j=0; j<n; j++)
M[i][j] =rand()%10;
      affichage(M,n,n);
      printf("somme de la première diagonale : %d\n",somme_diagonale(M,n));
      system("pause");
}