added calculations for harmonic oscillator, using various implicit schemes.