simulate_data function

Simulate design matrix