Fix calculation of variances

[goerz]

Variances are currently calculated by averaging over all trajectories. This is wrong (see http://michaelgoerz.net/notes/calculating-expectation-values-and-variances-for-quantum-trajectories.html)

Instead, we should keep track of the expectation value for each operator, and the operator squared. The variance can the easily be calculated on the fly for plotting.