The initial-boundary value problem of two-dimensionalincompressible fluid flow in stream function form is considered.A prediction-correction Legendre spectral scheme is proposed, which iseasy to be performed.The numerical solution possesses the accuracy of second-order in time and higher order in space. Thenumerical experiments show the high accuracy of this approach.