  We show how to find upper and lower bounds to the median of a gamma
distribution, over the entire range of shape parameter $k > 0$, that are the
tightest possible bounds of the form $2^{-1/k} (A + Bk)$, with closed-form
parameters $A$ and $B$. The lower bound of this form that is best at high $k$
stays between 48 and 50 percentile, while the uniquely best upper bound stays
between 50 and 55 percentile. We show how to form even tighter bounds by
interpolating between these bounds, yielding closed-form expressions that more
tightly bound the median. Good closed-form approximations between the bounds
are also found, including one that is exact at $k = 1$ and stays between 49.97
and 50.03 percentile.
