Add basic quadratic optomechanics code that reproduces NJP paper