Calculating Mardia's coefficient of each sample of k dimensions and representing in a vector

Question

I have generated a sample of 4 multivariate populations where each population has n=5 and k=2. My R version is 4.0.3. I want to calculate mardia\'s coefficients of skewness