In this paper we present properties of an algorithm to determine the maximum likelihood estimators of the covariance matrix when two processes jointly affect the observations. Additionally, one process is partially modeled by a compound symmetry structure. We perform a simulation study of the properties of an iteratively determined estimator of the covariance matrix.