function value = quartic ( x )
  value = 2.0 * x.^4 - 4.0 * x.^2 + x + 20.0;
  return
end

