f_num withsol('sol1',freq,setind(lambda,1)) "First eigenfrequency"
f_r real(f_num) "Resonance frequency (of 1st mode)"
omega_r 2*pi*f_r "Resonance angular frequency"
alpha_r imag(f_num)*2*pi "Attenuation coefficient"
Q_r omega_r/(2*alpha_r) "Q factor"
df_r f_r/Q_r "Resonance half power width"
P_simple 20*log10(abs(omega_r^2/((ta.iomega)^2+2*alpha_r*ta.iomega+omega_r^2)))-20*log10(abs(omega_r/(2*i*alpha_r))) "Simple harmonic oscillator model fit"