An accurate and effective numerical method is presented to solve the Maxwell-Bloch equations,which describe the optical pulse propagation and interaction with homogeneously broadened two-level medium.The convergence and stability of this numerical method are also proved by comparing with the analytic solutions derived under special condition,and the method maintains its erroneous exponents and is applicable.Simulations of arbitrary conditions are discussed by employing this new method,and the evolutions of pulse and Bloch vectors for different detunings,relaxation times and initial input pulse power are analyzed.The established numerical method can be used to solve the Maxwell-Bloch equations and their corrections.