def t60_arau(Sx, Sy, Sz, alpha, volume, c=SOUNDSPEED):
    a_x = -np.log(1 - alpha[0])
    a_y = -np.log(1 - alpha[1])
    a_z = -np.log(1 - alpha[2])
    St = np.sum(np.array([Sx, Sy, Sz]))
    A = St * a_x**(Sx/St) * a_y**(Sy/St) * a_z**(Sz/St)
    t60 = 4.0 * np.log(10.0**6.0) * volume / (c * A)
    return t60
Reverberation time according to Arau. [#arau]_
dataset/ETHPy150Open python-acoustics/python-acoustics/acoustics/room.py/t60_arau