Please see if this helps. It is giving pretty good results around 30. to poisson let maxv 52 let minv 1 let lamda 30 let z minv + random-poisson (lamda - minv) show z end Let me know how it worked out for you. Thank you.