script to calculate ewens sampling distribution