sim function

New generic for simulating matrix distributions