This study focuses on two-dimensional flows of airfoils, looking for the wake position for a uniform flow over an airfoil. The flow is assumed to satisfy the potential theory. First of all, the Boundary Integral Equation is applied to solve the velocity potential on the boundary of the airfoil. Once the strength of the velocity potential is solved, it is substituted into the analytical formula to find the flow field velocity components at the collocation points. The constraint of normal velocity components to be zero are then used to determine the location of the wake. The Joukowski airfoil and NACA airfoils are used as test cases.