Calcolo media valori diagonale destra matrice

#include <iostream>
#include <ctime> // srand
#include <cstdlib> // srand
#include <math.h> // round
#include <iomanip> // setprecision

using namespace std;

int main()
{
    int righe,colonne,elementi,i=0,j=0,media=0,k; //i -> righe, j -> colonne, k -> altezza lato della tangente
    double tangente; //qui verrà memorizzata la tangente dell'angolo beta o gamma
    cout << "Digitare la grandezza della matrice ed il numero degli elementi\n";
    cout << "Righe :";
    cin >> righe;
    cout << "Colonne :";
    cin >> colonne;
    cout << "Numero di elementi :";
    cin >> elementi;
    int matrice[righe][colonne];
    srand(time(0));
    while (i<righe)
    {
        j=0;
        while (j<colonne)
        {
            matrice[i][j]=rand()%elementi;
            cout << matrice[i][j] << " ";
            j++;
        }
        cout << "\n";
        i++;
    }

    i=0; // 'i' adesso viene usato come contatore e 'j' come variabile temporanea;

    if (righe<colonne)
    {
        tangente=double(righe)/colonne; // tangente di gamma
        j=colonne;
        while (i<j)
        {
            k=round((j-i)*tangente)-1;
            media+=matrice[k][i];
            cout << "matrice[" << k << "][" << i << "]=" << matrice[k][i] << "\n";
            i++;
        }

    } else
    {
        tangente=double(colonne)/righe; // tangente di beta
        j=righe;
        while (i<j)
        {
            k=round((j-i)*tangente)-1;
            media+=matrice[i][k];
            cout << "matrice[" << i << "][" << k << "]=" << matrice[i][k] << "\n";
            i++;
        }
    }
    cout << "La media dei valori della diagonale destra della matrice e': " << setprecision(3) << double(media)/i;
    return 0;
}