def sound_exposure(pressure, fs, axis=-1):
    """
    Sound exposure :math:`E_{T}`.

    :param pressure: Instantaneous sound pressure :math:`p`.
    :param fs: Sample frequency :math:`f_s`.
    :param axis: Axis.

    .. math:: E_T = \\int_{t_1}^{t_2} p^2 (t) \\mathrm{d}t

    """
    return (pressure**2.0/fs).sum(axis=axis)
dataset/ETHPy150Open python-acoustics/python-acoustics/acoustics/standards/iso_tr_25417_2007.py/sound_exposure/original