A new method for computing laminar-turbulent transition and turbulence in compressible boundary layers is proposed. It is especially useful for computation of laminar-turbulent transition and turbulence starting from small-amplitude disturbances. The laminar stage, up to the beginning of the breakdown in laminar-turbulent transition, is computed by parabolized stability equations (PSE). The direct numerical simulation (DNS) method is used to compute the transition process and turbulent flow, for which the inflow condition is provided by using the disturbances obtained by PSE method up to that stage. In the two test cases including a subsonic and a supersonic boundary layer, the transition locations and the turbulent flow obtained with this method agree well with those obtained by using only DNS method for the whole process. The computational cost of the proposed method is much less than using only DNS method.