Fix the formulas for standard deviations
Implement rand_gamma
Another one
Use new function
Use dp
use old function
Fix dp
Now it fails all the time
Improve the formulas