Vicinal Si (001) surfaces are interesting because they are good substrates for the growth of III-V semiconductors. Spots in RHEED patterns from vicinal surfaces are split due to scattering from ordered step arrays and this splitting can be used to determine the misorientation angle, using kinematic arguments. Kinematic theory is generally regarded to be inadequate for the calculation of RHEED intensities; however, only a few dynamical RHEED simulations have been attempted for vicinal surfaces. The multislice formulation of Cowley and Moodie with a recently developed edge patching method was used to calculate RHEED patterns from vicinal Si (001) surfaces. The calculated patterns are qualitatively similar to published experimental results and the positions of the split spots quantitatively agree with kinematic calculations.RHEED patterns were calculated for unreconstructed (bulk terminated) Si (001) surfaces misoriented towards [110] ,with an energy of 15 keV, at an incident angle of 36.63 mrad ([004] bragg condition), and a beam azimuth of [110] (perpendicular to the step edges) and the incident beam pointed down the step staircase.