changed to const pass-by-reference signature
factorial and fact_div_fact functions live here
variable declaration fixed
take MoleculeSTS instead of Molecule
actually, it should take Particles as input
omegaint (IPL needs fixing and no ESA yet)