import numpy as np
import pyfar as pf
from pyrato.analytic import rectangular_room_rigid_walls
L = np.array([8, 5, 3])/10
source_pos = np.array([5, 3, 1.2])/10
receiver_pos = np.array([1, 1, 1.2])/10
rir, _ = rectangular_room_rigid_walls(
    L, source_pos, receiver_pos,
    reverberation_time=1, max_freq=1e3, n_samples=2**16,
    speed_of_sound=343.9)
pf.plot.time_freq(rir)