function d = der(xx)
  es=0.002;
d=es*(4*xx^3)+1;
end