The Shanks transformation is a powerful nonlinear extrapolation method that is used to accelerate the convergence of slowly converging, and even diverging, sequences {A n }. It generates a two-dimensional array of approximations A ( j) n to the limit or anti-limit of {A n } defined as solutions of the linear systemswhereβ k are additional unknowns. In this work, we study the convergence and stability properties of A ( j) n , as j → ∞ with n fixed, derived from general linear sequencesi=0 β ki n γ k −i as n → ∞, where ζ k = 1 are distinct and |ζ 1 | = · · · = |ζ m | = θ, and γ k = 0, 1, 2, . . .. Here A is the limit or the anti-limit of {A n }. Such sequences arise, for example, as partial sums of Fourier series of functions that have finite jump discontinuities and/or algebraic branch singularities. We show that definitive results are obtained with those values of n for which the integer programming problems max s 1 ,...,s m m k=1 ( γ k )s k − s k (s k − 1) , subject to s 1 ≥ 0, . . . , s m ≥ 0 and m k=1 s k = n, 123 726 A. Sidi have unique (integer) solutions for s 1 , . . . , s m . A special case of our convergence result concerns the situation in which γ 1 = · · · = γ m = α and n = mν with ν = 1, 2, . . . , for which the integer programming problems above have unique solutions, and it reads A