© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This paper presents a novel and accurate method to implement the Kutta condition in solving subsonic (subcritical) inviscid isentropic compressible flow over isolated airfoils using the stream function equation. The proposed method relies on body-fitted grid generation and solving the stream function equation for compressible flows in computational domain using finite-difference method. An expression is derived for implementing the Kutta condition for the airfoils with both finite angles and cusped trailing edges. A comparison of the results obtained from the proposed numerical method and the results from experimental and other numerical methods reveals that they are in excellent agreement, which confirms the accuracy and correctness of the proposed method.