Goal

To simulate how the oscillator works.

Methods

1. Establish ODE equations based on Mass-action law;
2. Investigate reasonable parameter sets from previous researches;
3. Simulation;

Results

We solved these DDEs with R language. We also went one step further. We simulated the situation in which lag obeys a specific gaussian distribution, and the lag $\tau$changes in every certain interval. We hoped by running a random test, we could get closer to real life situation. The results are below.