We study the Lanczos method for solving symmetric linear systems with multiple right-hand sides. First, we propose a numerical method of implementing the Lanczos method, which can provide all approximations to the solution vectors of the remaining linear systems. We also seek possible application of this algorithm for solving the linear systems that occur in continuation problems. Sample numerical results are reported. Copyright (C) 2002 John Wiley Sons, Ltd