function y=HzToErbRate(x)


y=(21.4*log10(4.37e-3*x+1));
