The nonlinear stability problem in nonparallel boundary layer flow for two-dimensional disturbances was studied by using a newly presented method called Parabolic Stability Equations (PSE). A series of new modes generated by the nonlinear interaction of disturbance waves were tabulately analyzed, and the Mean Flow Distortion (MFD) was numerically given. The computational techniques developed, including the higher-order spectral method and the more effective algebraic mapping, increased greatly the numerical accuracy and the rate of convergence. With the predictor-corrector approach in the marching procedure, the normalization condition was satisfied, and the stability of numerical calculation could be ensured. With different initial amplitudes, the nonlinear stability of disturbance wave was studied. The results of examples show good agreement with the data given by the DNS using the full Navier-Stokes equations.