Abstract: A method which Combines the numerical flux of high accurate TVD scheme with multistage method in time approaching is developed to calculate the problems of 2-D steady flow. Some special treatments are given to accelerate convergency. The numerical results of transonic flow over NAGA 0012 airfoil are obtained by this method. They show that this method is high accurate, of high resolut-ion nd time-saving.