pargen function

Parameter simulation