The problem of integrating the time-dependent Schrodinger equation (TDSE) describing the interaction of a polyelectronic atom with a laser pulse is treated by expanding the time-dependent wave function %(r, t) in terms of wave functions 4"E computed for discrete, autoionizing, and scattering states separately. The TDSE is transformed into a system of coupled first-order differential equations with time-dependent coeScients, whose number (in the thousands), necessary for convergence to be reliable, depends mainly on the degree of the contribution of the continuous spectrum, as a function of the frequency and strength of the field. This approach allows the systematic incorporation of the significant electronic structure, electron correlation, and spectral characteristics of each N-electron system under investigation. Furthermore, since the free-electron function is computed numerically in the polarized core potential of the remaining (N -1)-electron atom, properties such as the angular distribution and partial above-threshold ionization (ATI) of the photoelectron are directly computable. We present results from the application of our methods to H and Li . For the applications to H, which served as testing grounds for the method, the state-specific wave functions for discrete and continuum states were obtained numerically, for n and I up to 12 and 5, respectively, and for positive energies up to a= 34 eV with I up to 6. When comparisons with other time-independent and time-dependent results are possible, very good agreement is observed. On the other hand, our calculations do not confirm recent experimental results on absolute ionization rates for laser pulses of 248 nm. For Li, our results on ATI for photon energy tie=1.36 eV demonstrate the effects of initial-state electron correlation and of final-state fieldinduced coupling of open channels (the Li 1s 2s Sand 1s 2p P'), as a function of field intensity.