Ztrans (intop_in(ta.p_t)/intop_in(1)-intop_out(ta.p_t)/intop_in(1))/(intop_mid(w)/intop_mid(1))/(rho0*c0) Transfer impedance
iomega0 i*2*pi*freq Complex angular frequency
kv sqrt(-rho0*iomega0/mu0) Viscous wave number
kth sqrt(-iomega0*rho0*Cp0/kappa0) Thermal wave number
Gv -besselj(2,kv*a)/besselj(0,kv*a) Help variable 
Gth -besselj(2,kth*a)/besselj(0,kth*a) Help variable
Zc 1/sqrt((gamma0-(gamma0-1)*Gth)*Gv) Characteristic impedance
kc 2*pi*freq/c0*sqrt((gamma0-(gamma0-1)*Gth)/Gv) Complex wave number
Zlin_NOec 2*i*Zc*sin(kc*tp/2) Linear transfer impedance (no end correction)
EC 16*a/(3*pi) End correction
ECfit 0.8 End correction fitting factor
Zlin 2*i*Zc*sin(kc*(tp+ECfit*EC)/2) Linear transfer impedance (analytical)
Ztrans_rms (abs(intop_in(ta.p_t))/intop_in(1)-abs(intop_out(ta.p_t))/intop_in(1))/(abs(intop_mid(w))/intop_mid(1))/(rho0*c0) Transfer impedance from model (RMS)