simulate function

Simulate data from a structural equation model