We approximate the solution of initial boundary value problems for nonlinear parabolic equations. In space we discretize by finite element methods. The discretization in time is based on linear multistep schemes. One part of the equation is discretized implicitly and the other explicitly. The resulting schemes are stable, consistent and very efficient, since their implementation requires at each time step the solution of a linear system with the same matrix for all time levels. We derive optimal order error estimates. The abstract results are applied to the Kuramoto-Sivashinsky and the Cahn-Hilliard equations in one dimension, as well as to a class of reaction diffusion equations in R ν , {\mathbb {R}} ^{\nu }, ν = 2 , 3. \nu = 2, 3.