f=[0.636 0.716 0.796 0.82 0.875 0.907 0.955 0.986 1.034 1.114 
    1.194 1.273 1.353 1.384 1.432 1.512 1.591 1.671 1.702 1.751 1.832 1.91];