plot(f/(60 + ((0.177+f*0.0253)^2)*f*(0.71/0.9/0.85)), f=21..33)