forgot to write the unit
primitive code for n-gamma distrimination