[r,l,g,c]=RLGCExpQ22(f);

function [r,l,g,c]=RLGCExpQ22(f);
r0c=252.5;
ac=0.2643;
l0=763.2e-6;
linf=575.7e-6;
fm=386e3;
b=0.794;
g0=26.7e-9;
ge=0.8805;
cinf=44.1e-9;
c0=75.1e-9;
ce=0.245;
r=sqrt(sqrt(r0c^4+ac.*f.^2));
l=(l0+linf*(f./fm).^b)./(1+(f./fm).^b);
g=g0.*f.^ge;
c=cinf+c0.*f.^(-ce);