This work presents the analytic solution of flow of a viscous fluid between two stretching disks with slip boundaries. Suitable similarity transformations are used to normalize the governing system. Optimal homotopy analysis method has been used as a solution technique by considering one parameter; two parameters; three parameters and one step optimal HAM approaches. The solutions obtained through all these approaches are compared. It is observed that, although the solution series (obtained through Optimal HAM) attains convergence at initial few terms but, on the other hand the computational resources,required to acquire solution, increases indefinitely as we increase the order of approximation.The current analytic solution is compared with a numerical technique and it is found that both solutions are perfectly matched. The effects of material parameters on different physical quantities are illustrated through graphs and tables.