This is the initial release of this package. The package has been restructured such that the organisation is tidier and more easy to extend. The current simulation supports the modelling of free dynamics plus x^4 non-linearity in the form of the harmonic potential, pulsed squeezing, and driving and cooling on resonance and at double the frequency.