function y = f(x)
  y = x + x^2 * sin(x);
