A FORTRAN program is given to calculate the run length distribution (RLD), the average run length (ARL), and the standard deviation of the run length (SDRL) for residual control charts used to monitor autocorrelated process output. RLD, ARL, and SDRL values are calculated for processes that can be modeled by pure autoregressive models of orderp(AR(p)), pure moving-average models of order 1 (MA(1)), and mixed autoregressive moving-average models of orderspand 1 (ARMA(p,1)), given that the assignable cause to be detected is a step shift in the process mean.