double f5 ( double x )
{
  double value;

  value = pow ( x, 4 ) + 2.0 * pow ( x, 2 ) + x + 3.0;

  return value;
}