Started to add module for harmonically trapped fermi gas.  Also polynomial approximations off polylog functions.