added robust multinomial calculation used in Jie's model (but not in a class, darn)
Merge pull request #37 from seacode/jimdev

added robust multinomial calculation used in Jie's model (but not in a c...