function y = golden ( x )

  y = x - 1.0 - 1 ./ x;

  return
end

