rand_diag_mat function

Simulate a random diagonal variance matrix