Gamma

使用牛顿迭代法求解 Gamma 分布的逆累积分布函数（分位数）。

BUG FIX: 原实现的循环条件 |cdf - p| <= eps Or i = 100 是反的 —— 它会在估计值已经足够精确（或恰好迭代到第 100 次）时继续循环，反而在尚未收敛时立即退出， 因此通常只做了一两次迭代就返回，分位数结果不可靠。这里改为标准的收敛判据 （未收敛且迭代次数未用尽时继续），并加入 pdf=0 与越界保护。