The incomplete beta function is defined as where Beta(p, q) is the beta function. Dutka (1981) gave a history of the development and numerical evaluation of this function. In this article, an algorithm for computing first and second derivatives of Ix,p,q with respect to p and q is described. The algorithm is useful, for example, when fitting parameters to a censored beta, truncated beta, or a truncated beta-binomial model.