Add basic code to compute quadratic optomechanics.  Duplicates NJP figures