Two dimensional and axisymmetric supersonic turbulent boundary layer with pressure gradient are calculated by using implicit finite difference method. Fluctuating terms in the equations are removed by using eddy viscosity and eddy conductivity concepts. Prandtl mixing length formulas are modified to account for pressure gradient and compressibility effects. The computational scheme used in the study is Keller Box Scheme which is efficient and accurate for this type of numerical calculation. Comparisons with published experimantal data show good agreements for boundary layers with favorable or mild adverse pressure gradient. And it is also found that the models of eddy viscosity and turbulent Prandtl number used in the study should be modified to account for the non-equiliblium flow effect of boundary layer developing rapidly with strong adverse pressure gradient.