# Bore and reed defition
D = mv.ReedDynamics(wr=2*np.pi*1500., qr=0.4, kr=8e8)
Ze = mac.Cylinder(L=.5, r=7e-3, radiates=False, nbmodes=10)