def ai(vp,rho):
    """
    Computes de acoustic impedance

    Parameters
    ----------
    vp : array
        P-velocity.
    rho : array
        Density.

    Returns
    -------
    ai : array
        Acoustic impedance.
    """
    
    ai = vp*rho
    
    return (ai)