We propose an effective and useful numerical simulation scheme for the investigation of the ultra-fast laser pulses in tenuous plasmas. The accuracy of the method is tested by numerical examples. We check some special examples to investigate the laser envelope evolving and modulation in plasmas. Asymmetric two-peak modulation structure is found and its underlying physics is analyzed. The advantages and shortages of the method are also discussed.