Rolling contact fatigue (RCF) is an important problem in the railway field. In rail, RCF cracks may be initiated by the near-surface plastic deformation and start by growing down at a shallow angle to the surface, so called co-planar crack growth, and may form branch cracks when they reach a certain length. This paper attempts to simulate the co-planar crack growth and branching in various conditions. The numerical simulation was based on the K-values obtained by BEM, Hourlier and Pineau's criterion and crack growth law due to the experiments.