Abstract. In this work, we give a detailed convergence and stability analysis for the author's generalized Richardson extrapolation process GREP (m) as this is being applied to linearly convergent or divergent infinite sequences {A n }, whereThe quantity we would like to compute is A, whether it is the limit or antilimit of {A n }. Such sequences arise, for example, as partial sums of power series and of Fourier series of functions that have algebraic and/or logarithmic branch singularities. Specifically, we define the GREP (m) approximation A (m,j) n to A, with n = (n 1 , . . . , n m ) and α > 0, via the linear systemsn k , whereβ ki are additional unknowns. We study the convergence and stability properties of A (m,j)is a true convergence acceleration method for the sequences considered. In addition, we show that GREP (m) is stable for the case being studied, and we also quantify its stability properties. The results of this work are the first ones pertaining to GREP (m) with m > 1.