function fpx = fp2( x )
  fpx = zeros(2,1);
  fpx(1) = 4 * x(1) - 4.2 * x(1).^3 + 6.0 * x(1).^5 + x(2);
  fpx(2) = x(1) + 2.0 * x(2);
  return
end

