Mathematical model is developed for unsteady boundary layer flow of Carreau fluid. Numerical solutions are constructed for momentum and energy equations. A comparison with the available published literature in limiting cases is performed. Impact of various pertinent parameters on the velocity and temperature profiles are discussed graphically.