Квадратная матрица

#include <iostream>
using namespace std;

#define maxN 100

int main()
{
    double matrix[maxN][maxN]={0};
    int n=0;

    cout<<"Enter n"<<" ";
    cin>>n;
   
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {  
            cout<<"mas["<<i<<","<<j<<"]=";
            cin>>matrix[j];
        }

    double sum=0;
    for(int k=0;k<n;k++)
        sum=sum+matrix[k][k];

    cout<<"Sum="<<sum<<endl;

    return 0;
}