Two dimensional parabolic stability equations (PSE) are numerically solved using expansions in orthogonal functions in the normal direction.The Chebyshev polynomials approximation,which is a very useful form of orthogonal expansions, is applied to solving parabolic stability equations. It is shown that results of great accuracy are effectively obtained.The availability of using Chebyshev approximations in parabolic stability equations is confirmed.