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");
}