function value = center ( f, x, h )

  value = ( f ( x + h ) - f ( x - h ) ) ./ h / 2.0;

  return
end

