Numerical instability in B&D models

Abstract

When computing the expected value of time till extinction of a Birth and Death process, the usual textbook approach results in an extreme case of numerical ill-conditioning, which prevents us from getting accurate answers beyond the first few low-lying states; in this brief note we present a potential solution, together with a novel derivation of related formulas.