'Compute orthogonalized MA coefficient matrices.
For this purpose a matrix  P is used which fulfills
:math:`\Sigma_u = PP^\prime`. P defaults to the Cholesky
decomposition of :math:`\Sigma_u`
Parameters
maxn : int
Number of coefficient matrices to compute
P : ndarray (neqs x neqs), optional
Matrix such that :math:`\Sigma_u = PP\'`. Defaults to Cholesky
decomposition.
Returns
coefs : ndarray (maxn x neqs x neqs)'

def orth_ma_rep(self, maxn=10, P=None):

    return orth_ma_rep(self, maxn, P)
