The physical problem under consideration is the boundary layer problem of an incompressible, laminar flow, taking place over a flat plate in the presence of a pressure gradient and radiation. For the mathematical formulation of the problem, the partial differential equations of continuity, energy, and momentum are taken into consideration with the boundary layer simplifications. Using the dimensionless Falkner–Skan transformation, a nonlinear, nonhomogeneous, coupled system of partial differential equations (PDEs) is obtained, which is solved via the homotopy analysis method. The obtained analytical solution describes radiation and pressure gradient effects on the boundary layer flow. These analytical results reveal that the adverse or favorable pressure gradient influences the dimensionless velocity and the dimensionless temperature of the boundary layer. An adverse pressure gradient causes significant changes on the dimensionless wall shear parameter and the dimensionless wall heat-transfer parameter. Thermal radiation influences the thermal boundary layer. The analytical results are in very good agreement with the corresponding numerical ones obtained using a modification of the Keller’s-box method.