A method is presented for estimating definite and indefinite integrals over a density function, such as a local density of states, defined by a three-term recurrence relation. This may be generated, for example, by the 'recursion method' applied to some Hamiltonian, and properties of the approximation are given, and the results derived, in that context.