We present a numerically exact calculation of rovibrational levels of a five-atom molecule. Two contracted basis Lanczos strategies are proposed. The first and preferred strategy is a two-stage contraction. Products of eigenfunctions of a four-dimensional (4D) stretch problem and eigenfunctions of 5D bend-rotation problems, one for each K, are used as basis functions for computing eigenfunctions and eigenvalues (for each K) of the Hamiltonian without the Coriolis coupling term, denoted H-0. Finally, energy levels of the full Hamiltonian are calculated in a basis of the eigenfunctions of H-0. The second strategy is a one-stage contraction in which energy levels of the full Hamiltonian are computed in the product contracted basis (without first computing eigenfunctions of H-0). The two-stage contraction strategy, albeit more complicated, has the crucial advantage that it is trivial to parallelize the calculation so that the CPU and memory costs are independent of J. For the one-stage contraction strategy the CPU and memory costs of the difficult part of the calculation scale linearly with J. We use the polar coordinates associated with orthogonal Radau vectors and spherical harmonic type rovibrational basis functions. A parity-adapted rovibrational basis suitable for a five-atom molecule is proposed and employed to obtain bend-rotation eigenfunctions in the first step of both contraction methods. The effectiveness of the two methods is demonstrated by calculating a large number of converged J=1 rovibrational levels of methane using a global potential energy surface. (C) 2004 American Institute of Physics.