def ckm_wolfenstein(laC, A, rhobar, etabar):
rho = rhobar/(<NUM_LIT:1> - laC**<NUM_LIT:2>/<NUM_LIT>)<EOL>eta = etabar/(<NUM_LIT:1> - laC**<NUM_LIT:2>/<NUM_LIT>)<EOL>return np.array([[sqrt(<NUM_LIT:1> - laC**<NUM_LIT:2>)*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:6>*((-<NUM_LIT>)*eta + rho)*((<NUM_LIT>)*eta + rho)),<EOL>laC*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:6>*((-<NUM_LIT>)*eta + rho)*((<NUM_LIT>)*eta + rho)),<EOL>A*laC**<NUM_LIT:3>*((-<NUM_LIT>)*eta + rho)],<EOL>[-(laC*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:4>)) - A**<NUM_LIT:2>*laC**<NUM_LIT:5>*sqrt(<NUM_LIT:1> - laC**<NUM_LIT:2>)*((<NUM_LIT>)*eta + rho),<EOL>sqrt(<NUM_LIT:1> - laC**<NUM_LIT:2>)*sqrt(<NUM_LIT:1> -  A**<NUM_LIT:2>*laC**<NUM_LIT:4>) - A**<NUM_LIT:2>*laC**<NUM_LIT:6>*((<NUM_LIT>)*eta + rho),<EOL>A*laC**<NUM_LIT:2>*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:6>*((-<NUM_LIT>)*eta + rho)*((<NUM_LIT>)*eta + rho))],<EOL>[A*laC**<NUM_LIT:3> - A*laC**<NUM_LIT:3>*sqrt(<NUM_LIT:1> - laC**<NUM_LIT:2>)*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:4>)*((<NUM_LIT>)*eta + rho),<EOL>-(A*laC**<NUM_LIT:2>*sqrt(<NUM_LIT:1> - laC**<NUM_LIT:2>)) - A*laC**<NUM_LIT:4>*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:4>)*((<NUM_LIT>)*eta + rho),<EOL>sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:4>)*sqrt(<NUM_LIT:1> - A**<NUM_LIT:2>*laC**<NUM_LIT:6>*((-<NUM_LIT>)*eta + rho)*((<NUM_LIT>)*eta + rho))]])<EOL>
r"""CKM matrix in the Wolfenstein parametrization and standard phase
    convention.

    This function does not rely on an expansion in the Cabibbo angle but
    defines, to all orders in $\lambda$,

    - $\lambda = \sin\theta_{12}$
    - $A\lambda^2 = \sin\theta_{23}$
    - $A\lambda^3(\rho-i \eta) = \sin\theta_{13}e^{-i\delta}$

    where $\rho = \bar\rho/(1-\lambda^2/2)$ and
    $\eta = \bar\eta/(1-\lambda^2/2)$.

    Parameters
    ----------

    - `laC`: Wolfenstein parameter $\lambda$ (sine of Cabibbo angle)
    - `A`: Wolfenstein parameter $A$
    - `rhobar`: Wolfenstein parameter $\bar\rho = \rho(1-\lambda^2/2)$
    - `etabar`: Wolfenstein parameter $\bar\eta = \eta(1-\lambda^2/2)$