simulate_cosinor function

Simulate data from a cosinor model