Hola. No sé cual es el error en el siguiente código, ¿alguien podría ayudarme?

//Problema 1: Programa que cree una matriz nxn, ingrese datos y luego imprima la suma de los elementos de la diagonal principal

#include<iostream>
#include<cstdlib>

using namespace std;

void llenarymostrarM(int *matriz, int f, int c);
int sumadiagp(int *matriz, int f, int c);

int main()
{
    int *mat, fil, col;
    cout<<"El numero de filas es: ";
    cin>>fil;
    cout<<"El numero de columnas es: ";
    cin>>col;
    mat=new int[fil, col];
    llenarymostrarM(mat, fil, col);
    cout<<endl;
    cout<<"La suma de los elementos de la diagonal principal es: "      <<sumadiagp(mat, fil, col)<<endl;
    delete []mat;
    system("pause");
}

void llenarymostrarM(int *matriz, int f, int c){
     int i,j;
     cout<<"Ingrese los elementos de la matriz: ";
     for(i=0; i<f; i++){
              for(j=0; j<c; j++){                                         
                       cin>>matriz[i,j];
                       }
     }
     cout<<endl;
     cout<<"La matriz es: ";
     for(i=0; i<f; i++){
              for(j=0; j<c; j++){
                       cout<<matriz[i,j]<<endl;
                       }
     }
}                           

int sumadiagp(int *matriz, int f, int c){
      int i,j, sdp=0;
      for(i=0; i<f; i++){
              for(j=0; j<c; j++){
                       sdp=sdp+matriz[i,i];
                               }
              }
      return sdp;
}

Por ejemplo cuando pongo una matriz 2x2: 1, 2, 3, 4, me sale que la matriz es: 3 4 3 4 y la suma de la diagonal principal es 14 

Muchas gracias por la ayuda. Creo que este el código final además le agregue otro problema jeje.

#include<iostream>
#include<cstdlib>

using namespace std;

void llenarymostrarM(int **matriz, int o);
int mayor(int **matriz, int o);
int sumadiagp(int **matriz, int o);

int main()
{
    int **mat, orden;
    cout<<"El orden de la matriz es: ";
    cin>>orden;
    mat=new int*[orden];
    llenarymostrarM(mat, orden);
    cout<<endl;
    cout<<"El mayor elemento es: "<<mayor(mat, orden)<<endl;
    cout<<"La suma de los elementos de la diagonal principal es: "<<sumadiagp(mat, orden)<<endl;
    delete []mat;
    system("pause");
}

void llenarymostrarM(int **matriz, int o){
     int i,j;
     cout<<"Ingrese los elementos de la matriz: ";
     for(i=0; i<o; i++){
              for(j=0; j<o; j++){                                        
                       cin>>matriz[j];
                       }
     }
     cout<<endl;
     cout<<"La matriz es: ";
     for(i=0; i<o; i++){
              for(j=0; j<o; j++){
                       cout<<matriz[j]<<endl;
                       }
     }
}                            

int mayor(int **matriz, int o){
    int i,j, mayor=matriz[0][0];
    for(i=1; i<o; i++){
              for(j=0; j<o; j++){
                       if(matriz[j]>mayor){
                                      mayor=matriz[j];
                                      }
              }
    }
    return mayor;
}

int sumadiagp(int **matriz, int o){
      int i,j, sdp=0;
      for(i=0; i<o; i++){
              for(j=0; j<o; j++){
                       if(i==j)
                       sdp=sdp+matriz[j];
                               }
              }
      return sdp;
}

Con respecto al return 0, ¿es siempre obligatorio?, nunca lo he usado en otros códigos.  Donde dice matriz[j] en realidad debe ser matriz i[j], no sé porque se ve asi. Por último, cuando termino de ejecutar el programa sale que debe cerrarse 

Muchas gracias por la ayuda, recién por cuestiones de tiempo he podido entrar.