The splitoperator method solving the wavefunction evolution of atomic systems in ultrashort laser fields is introduced. It gives more accurate simulation results in comparison with the standard difference expression. Much information can be extracted from the timeresolved wavefunction. The ionization probability and highorder harmonic spectrum of a simple atomic model is calculated．