function f = fun(xx)
  es=0.002;
f=es*(xx^4-10000)+xx;
end