def CAT(N, rho, k):
from numpy import zeros, arange<EOL>cat = zeros(len(rho))<EOL>for p in arange(<NUM_LIT:1>, len(rho)+<NUM_LIT:1>):<EOL><INDENT>rho_p = float(N)/(N-p)*rho[p-<NUM_LIT:1>]<EOL>s = <NUM_LIT:0><EOL>for j in range(<NUM_LIT:1>, p+<NUM_LIT:1>):<EOL><INDENT>rho_j = float(N)/(N-j)*rho[j-<NUM_LIT:1>]<EOL>s = s + <NUM_LIT:1.>/rho_j<EOL><DEDENT>cat[p-<NUM_LIT:1>] = s/float(N) - <NUM_LIT:1.>/rho_p<EOL><DEDENT>return cat<EOL>
r"""Criterion Autoregressive Transfer Function :

    .. math::  CAT(k) = \frac{1}{N} \sum_{i=1}^k \frac{1}{\rho_i} - \frac{\rho_i}{\rho_k}

    .. todo:: validation