Funzione con array bidimensionali

Scrivere una funzione C++ che, data una matrice
quadrata di numeri double, restituisce il valore più grande, il valore più piccolo
e la media di tutti gli elementi della matrice.
ecco il mio codice ma mi dà sempre lo stesso risultato per la media
ho fatto solo il codice per la media e vorrei risolvere questo prima di passare agli altri

#include<iostream>
using namespace std;
const int N=4;
void funzione(double [N][N]);
int main()
{
    double a[N][N]={{3.9,50.80,6.7,8},{2.23,3.456,45,6},{4,8,9,19.405},{3,5,6,7.8}};
    funzione(a);
    return 0;
}

void funzione(double x[N][N])
{
    int i,j;
    double somma=0;
    double media=0;
    for(i=0;i<N;i++)
    {
        for(j=0;j<N;j++)
        {
            somma+=x[i][j];
            j++;
        }
        i++;
    }
media=somma/(N*N);
cout<<media;

}