Two methods are presented to solve transonic airfoil flow problems. The first method is based on the integral equation solution of the full-potential equation in terms of the velocity field. A shock capturing-shock fitting scheme has been developed. In the shock-fitting part of the scheme, shock panels are introduced at the shock location. The shock panels are fitted by using the Rankine-Hugoniot relations. The second method is a coupling of the integral solution of the full-potential equation with the pseudotime integration of Euler equations, which are used in a small embedded region around the shock. This scheme is caUed the integral equation-embedded Euler scheme. The two methods are applied to NACA 0012 and 64A010A airfoils over a wide range of Mach numbers, and the results are in good agreement with the experimental data and other computational results. The schemes converge within a number of iterations that is one order of magnitude less than the finite-difference schemes.