Present paper addresses the numerical study of Sisko fluid model over stretching cylinder with variable thermal conductivity. The governing equations are simplified by incorporating the boundary layer approximations. After employing suitable similarity transformations partial differential equations are reduced to ordinary differential equations. To obtain numerical solution shooting method in conjunction with Runge-Kutta-Fehlberg method is used. For the analysis of model, variations due to different physical parameters involved in momentum and heat equations are reflected through graphs. Also, the effects of physical parameters on skin-friction coefficient and Nusselt number are represented through graphs as well as tables.
