A new technique to model ultrashort optical pulses is proposed and verified. The technique uses Pade approximant to account for the fast pulse propagational variations. Numerical parameters of the technique have been tested and it was shown that the method is simple, very stable, and accurate in modeli ng ultrashort optical pulses in long propagation interaction
