function value = goat ( x )

  value = - 0.5 * x .* sqrt ( 400 - x.^2 ) ...
    + ( x.^2 - 200 ) .* acos ( x / 20.0 ) ...
    + 100 * pi / 2;

  return
end