First session finish!

Add the energy function of SG potential. The initial value can be set
at pericenter. Start coding the colliding part.
Add massFallback.c and start to beautify the code