Calculate necessary covariances for ARMA(p,q) process up to kappa(KappaSize,KappaSize)

procedure ARMAKappa(const Data: TVec; const Phi: TVec; const Theta: TVec; const Cov: TMtx; const KappaSize: Integer); overload;