Определить сумму элементов матрицы, расположенных на главной диагонали

#include <iostream> 
 
using namespace std; 
 
int main() 
{ 
    int a[4][4],
    int sum = 0, i, n, j; 
 
 
 for (i = 0; i < n; i++) 
for (j = 0; j < m; j++) 
 if: (i == j)
 { 
   sum +=a[i,j]
 }

#include <iostream>
 
using namespace std;
 
int main()
{
    int a[4][4];
    int S=0;
    int i,j;
    
    for(i=1;i<=4;i++)
    {
    for(j=1;j<=4;j++)
    {
                    cin>>a[i][j]; 
                     }
}
 
 
  for(i=1;i<=4;i++)
    {
    for(j=1;j<=4;j++)
    {
                    cout<<a[i][j]<<" "; 
                     }
                     cout<<endl;
}
 
  for(i=1;i<=4;i++)
    {
    for(j=1;j<=4;j++)
    {
                    if(i==j) {S=S+a[i][j];} 
                     }
}
 
cout<<S<<endl;
 
}

#include <iostream>
#include <iomanip>
 
int sumOfMainDiagonal(const int [][4], int, int);
int productOfAdverseDiagonal(const int [][4], int, int, int);
 
int main()
{
    int matr[4][4];
 
    for(int i = 0; i < 4; i++)
        for(int j = 0; j < 4; j++)
        {
            std::cout << "matr[" << i << "][" << j << "]: ";
            std::cin >> matr[i][j];
        }
 
    std::cout << "Исходная матрица: " << std::endl;
 
    for(int i = 0; i < 4; i++)
    {
        for(int j = 0; j < 4; j++)
            std::cout << std::setw(5) << matr[i][j];
 
        std::cout << std::endl;
    }
 
    std::cout << "Сумма эл-ов на главной диагонали: "
                  << sumOfMainDiagonal(matr,0,3) 
                  << std::endl;
    
        std::cout << "Произведение эл-ов на побочной диагонали: "
                  << productOfAdverseDiagonal(matr,3,0,3) 
                  << std::endl;
 
    return 0;
}
 
int sumOfMainDiagonal(const int matr[][4], int step, int end)
{
    if(step == end)
        return matr[step][step];
    else
        return matr[step][step] + sumOfMainDiagonal(matr,step+1,end);
}
 
int productOfAdverseDiagonal(const int matr[][4], int row, int column, int step)
{
    if(!step)
        return matr[row][column];
    else
        return matr[row][column] * productOfAdverseDiagonal(matr,row-1,column+1,step-1);
}