Update ForceFieldFunctions.f90

Added a number of functions for various derivatives of the harmonic energy. (not compiled)