function f=f_2(x)
   f=5000./(2.*x.^2.*(x-1)).*(sin(20.0./x)).^2;
end