We extend Zeilberger's fast algorithm for definite hypergeometric summation to non-hypergeometric holonomic sequences. The algorithm generalizes to differential and q-cases as well. Its theoretical justification is based on a description by linear operators and on the theory of holonomy.