This method computes tight lower & upper bounds for the renewal function. I
t is based on Riemann-Stieltjes integration, and provides bounds for solvin
g certain renewal equations used in the study of availability. An error ana
lysis is given for the numerical bounds when inter-renewal time distributio
ns are sufficiently smooth. Three examples are explored that demonstrate th
e accuracy of these computed numerical bounds.