cov_joint function

Covariance for joint distribution