Added partial support fot the 2d harmonic oscillator and inclued a
general interaction element solver.
Cleaning up the code style. Files, functions and variables not used has
been removed.