THE NUMERICAL METHOD FOR SOLVING THE RATE EQUATIONS OF ELECTRON'S AVERAGE OCCUPATION PROBABILITIES BY USING MULTIPLE-TIME-SCALE PERTURBATION THEORY