Useful libraries:

*  repaired factorial
 *  matrix exponential function
 *  round function