записать абсолютное отклонение каждого элемента матрицы от её среднего арифметического

#include <iostream>
#include <stdlib.h>
#include <iomanip>
#include <ctime>
using namespace std;
main ()
{
    int M[4][5],N[20],i,j,c,a,k;
    float sr,sum;
    setlocale(LC_ALL,"Rus");
    time_t t;
    srand((unsigned)time(&t));
    cout<<"Мартрица М"<<endl;
    for (i=0; i<4; i++)
    {
        for (j=0; j<5; j++)
        {
            M[i][j]=rand()%101-30;
            cout<<setw(6)<<M[i][j];
        }
        cout<<endl;
    }
    sum=0.0;
    for (int i = 0; i <4; i++)
    {
        for (int j = 0; j < 5; j++)
        {
       sum+=M[i][j];
        }
    }
       c=4*5;
       sr=sum/c;
       cout<<"Сумма: "<<sum;
       cout << endl;
       cout<<"Среднее арифтическое: "<<sr;
      k=-1;
       for (i=0; i<4; i++)
      {
        for (j=0; j<5; j++)
        {
        a=0;
        a=sr-M[i][j];k++;
         N[k]=a;    
        }
      }
       cout << endl; 
       cout<<"Абсолютное отклонение каждого элемента: "<<a;
       cout << endl;
       cout << "Одномерный массив N [20]";
       for (int i=0; i<k; i++)
       {
        cout<<N[k];
       }
}

#include <iostream>
#include <cstdlib>
#include <iomanip>
#include <ctime>
#include <cmath>
using namespace std;
 
int main ()
{
    srand((unsigned)time(0));
    const int n=4, m=5, k=n*m;
    int M[n][m], p=0;
    double N[k], avg=0.;    
    
    cout << "Matrix M:\n";
    for (int i = 0; i < n; i++)
    {
        for (int j=0; j < m; j++)
        {
            M[i][j]=rand()%19-9;
            avg+=M[i][j];
            cout << setw(6) << M[i][j];
        }
        cout << "\n";
    }    
    
    avg/=k;
    cout << "average=" << avg << "\n";
    
    cout << "Array N:\n";
    for (int j = 0; j < m; j++)    
        for (int i = 0; i < n; i++)        
        N[p++]=fabs(avg-M[i][j]);
       
    for (int i = 0; i < k; i++)       
        cout << N[i] << "  ";
    cout << "\n";
    
system("pause");
return 0;
}