Hi, So to summarize beta_deviate(a,b,x) or beta_deviate(a,b,x,eps) works in all cases. I hope you can get your example working and document this a bit (wherever that is done?) Dave