This article addresses the nonlinear radiation effect on two-dimensional oblique stagnation point flow in a porous medium. Constitutive equations of viscoelastic fluid are employed in the mathematical development of the relevant problem. The resulting nonlinear analysis is computed using Chebyshev Spectral Newton Iterative Scheme. A comparative study of present results with that of previous studies have been made in a limiting sense and shown through tabular values. Excellent agreement is noted which clearly shows that the used numerical scheme is stable and the results are highly accurate. Impact of sundry variables in the flow quantities of interest are discussed. It is observed that shearing parameter helps to increase the fluid velocity. Thermal boundary layer thickness can be controlled due to small value of radiation parameter and surface heating parameter. It is also noted that with the increase in value of porosity parameter K, the velocity increases but the momentum boundary layer thickness decreases in the region of stagnation point. Moreover, the streamlines are plotted to predict the flow pattern.