In order to simulate gene regulatory oscillators more effectively, Runge-Kutta RK integrators are adapted to the limit-cycle structure of the system. Taking into account the oscillatory feature of the gene regulatory oscillators, phase-fitted and amplification-fitted Runge-Kutta FRK methods are designed. New FRK methods with phase-fitted and amplification-fitted updated are also considered. The...