This paper present a novel system that simulates Time-of-Flight MR Angiography (TOF-MRA). In the current stage, pie focus on steady flow of blood. We first create a 3D blood vessel model obeying the branching rules of real vascular systems. In the model, laminar flows inside blood vessels are modeled using nested cylinders with a different velocity in each cylinder shell. Pseudo TOF-MRA images are generated by voxelizing the model. The voxelization is based on flow-related enhancement which is the basis of TOF-MRA. Slices of different resolutions and thicknesses can be generated with different settings for maximum velocity of blood flow as well as for the orientation of the phase and frequency encoding axes.