A recursion relationship for the modified Bessel function is used to develop recursive formulas for functions which are expressed as a Neumann series expansion. The method is applied to devise an algorithm for calculating the generalized

function. Comparisons with some alternative schemes are discussed and an error analysis is presented.