double f3 ( double x ) 
{
  double value;

  value = pow ( x, 4 ) + 2 * pow ( x, 3 ) - 14 * pow ( x, 2 ) + 2 * x + 1;

  return value;
}