A simulation algorithm of step or hole coverage for non-linear surface reaction on Chemical Vapor Deposition (CVD) is developed. The high speed simulation code for coverage is constructed by combining this algorithm with simple Monte Carlo simulation method. This code contained non-linear surface reaction could quantitatively explain the step coverages under high and low raw material (WF6: Tungsten (VI) fluoride) concentration on Tungsten thermal CVD.