skew_diag_mat function

Simulate a skewed diagonal variance matrix