Nonlinear Parabolized Stability Equations(NSPE) can be effectively used to study more throughly the transition process. NPSE can efficiently analyze the stability of a nonlinear region in transition process with low computational cost compared to Direct Numerical Simulation(DNS). In this study, NPSE in general coordinate system is formulated and a computer code to solve numerically the equations is developed. Benchmark problems for incompressible and compressible boundary layers over a flat plate are analyzed to validate the present code. It is confirmed that the NPSE methodology constructed in this study is an efficient and effective tool for nonlinear stability analysis.