In this study a numerical solution for the problem of pulsating flow in rigid tubes is described. The method applies to the case of known flow rate waveform, as opposed to Womersley solution where the pressure gradient was the known quantity. The solution provides the pressure gradient and wall shear stress waveforms as well as the instantaneous velocity profiles. Results show that the method can be used to study the blood flow characteristics in large arteries.