An implicit numerical integration scheme is proposed in the present paper for the solution of the transient flows. The numerical method is based on the spectral Chebyshev collocation technique in the direction normal to the disk and forward marching in time. Besides being free of the numerical oscillations caused by the discontinuities between the initial and boundary conditions inherent to the classical finite difference techniques, the devised technique benefits from the advantages of being robust, unconditionally stable, highly accurate and straightforward to implement. Unlike to the finite difference methods, it is also compact in the sense that it resolves the flow field without necessitating further transformations. The numerical algorithm developed here is applied to the classical time-dependent von Karman swirling flow due to a porous rotating disk impulsively set into motion which progresses into the well-known steady state after a long time. The energy equation is also treated by the method and the physical parameters of paramount interest as such the radial and tangential skin-friction coefficients, the torque and the rate of heat transfer from the disk surface are numerically calculated that are shown to approach their steady state counterparts.