sim_data function

Simulate data