PCMVarAtTime function

Calculate the variance covariance k x k matrix at time t, under a PCM model