How to make the mean, standard deviation and mode statistics of an NxN matrix
As pointed out in a comment,  this is easily found by a search.
The scipy library offers functions for all of these, either for the whole array or for specified axes. scipy.stats.describe will give the mean and variance all at once, and standard deviation is easily computed from variance. For mode use scipy.stats.mode
Alternatively for mean and standard deviation, there are numpy.mean and numpy.std respectively.
I would like to know how to do the average, standard deviation and mode statistics of an NxN matrix.
Can anybody suggest how to do this in python?