Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems, in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap. The L-leap method can better satisfy the leap condition. Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.