This article presents a numerical solution for the unsteady slip velocity flow of a thin viscous
liquid film over a heated horizontal stretching surface in the presence of variable heat flux. Similarity
transformations are used to transform the governing equations to a set of coupled non-linear ordinary
differential equations. The obtained differential equations are solved approximately by the homotopy perturbation
method (HPM). The effects of various parameters governing the flow and heat transfer in this
study are discussed and presented graphically. Comparison of numerical results is made with the earlier
published results under limiting cases.