I am using the MVN::mvn() function to test for univariate aswell as multivariate normality. With it, the functions delivers nicely arranged histograms to one big plot.