vcov.joint function

Extract the variance-covariance matrix from a joint fit.