The distribution of a sum S of independent binomial random variables, each with different success probabilities, is discussed. An efficient algorithm is given to calculate the exact distribution by convolution. Two approximations are examined, one based on a method of Kolmogorov, and another based on fitting a distribution from the Pearson family. The Kolmogorov approximation is given as an algorithm, with a worked example. The Kolmogorov and Pearson approximations are compared for several given sets of binomials with different sample sizes and probabilities. Other methods of approximation are discussed and some compared numerically. The Kolmogorov approximation is found to be extremely accurate, and the Pearson curve approximation useful if extreme accuracy is not required.