An integral method, using Coles combined inner and outer law as the velocity profile, is developed for calculation of turbulent boundary layers with blowing or suction and pressure gradients. The resulting ordinary differential equations are solved numerically for the distribution of skin friction coefficient and integral thickness along the surface. Comparisons of predicted skin friction coefficients with experimental data are made for a wide range of blowing and suction rates and for various pressure gradients, including adverse, zero and a strong favorable gradient. In addition to good agreement with experimental data for constant blowing fractions F, the method is also successfully tested on cases where the blowing fraction is variable with position. Predictions, in general, exhibit satisfactory agreement with the data. The integral method predictions are comparable to, or better than, a number of finite difference procedures in a limited number of cases where comparisons were made.