VMS function

Covariance matrix of M and S