"Problem 15.7.1" ahat(f) = sqrt(1/2) (f xi + d(f,xi)) ahatdag(f) = sqrt(1/2) (f xi - d(f,xi)) E(f) = i sqrt(hbar omega / 2 / epsilon0) (ahat(f) - ahatdag(f)) B(f) = sqrt(hbar omega mu0 / 2) (ahat(f) + ahatdag(f)) Psi = psi(t,x,y,xi) C = E(B(Psi)) - B(E(Psi)) C check(C == i hbar sqrt(mu0 / epsilon0) omega Psi) "ok" Run