In this study, stability analysis of boundary layers on airfoils is performed by using parabolized stability equations(PSE). Boundary layer edge conditions are obtained by compressible inviscid flow calculations. Mean velocity and temperature profiles of the laminar boundary layer are obtained by solving compressible boundary layer equations in generalized curvilinear coordinates with fourth order accuracy in the wall normal direction. Laminar mean flow profiles are used as input data for PSE to investigate growth rates of disturbances and stability characteristics. For the cases of boundary layer on NACA0012 and HSNLF(1)-0213 airfoils at Mach number 0.5, growth rates with respect to disturbance frequencies and profiles of disturbance amplitude are investigated. The effect of angle of attack on stability characteristics are examined at both upper and lower surfaces. The neutral stability curves, effect of Mach number and effect of airfoil section shapes are also analyzed.