Updated the calculations of m[1] and m[2] using n instead of n-1
Separated moment.hpp and moment.cpp out of usr.cpp and usrt.cpp