Do we need this anymore?  Since we are using boost anyway for gamma(2n,b) in the calculations for the truncated case, we could probably just always use the solver to find the right value of b.  As discussed in issue #325, the existing calculation isn't super accurate for small n.