vcov.garma_model function

Covariance matrix