function y=erb(x)
y=24.7*(4.37e-3*x+1);