Transition prediction of three-dimensional boundary layer is numerically investigated by the use of a modified e(exp N) method. Main idea is the corporation of a complex ray theory into the e(exp N) method so as to uniquely determine a wave number, which corresponds to a most amplified disturbance wave. From the calculation for a flow around a yawed circular cylinder with sweep angle 50 deg, it is found that the results agree well with a theoretical one obtained by N. Itoh (1986) excepting the N factor.