This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudosimilarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.