Update ForceFieldFunctions.f90

Added finite difference tests to some harmonic derivatives. (needs to compile)