sigmaZFun function

Sub-routine to get covariance of z-transformed variables