PCMVar function

Expected variance-covariance matrix for each couple of tips (i,j)