In this paper, the gear design program is presented. The profile of gears is created using classical mathematic formulations. In each gear, a kinematic joint is applied and one can define the 2D contact condition between gear pairs. Initial and boundary conditions such as force, torque, velocity, acceleration, etc. can be set. Thus, it is possible to analyze dynamic characteristics of gear pairs such as reaction moment and the variation of angular velocity. In order to find the optimal profile of gear pairs, two optimization methods based on design of experiments are inserted in the program; One is the Taguchi method and the other is the response surface analysis method. To verify the program, the rack & pinion gear is created and analyzed. Simulation results show that the developed program is useful and result data is reliable.