"Time evolution operator" U = exp(-i / hbar integral(H(x,t),t)) i hbar d(U,t) "Verify equation (1)" check(i hbar d(U,t) == H(x,t) U) "ok" Run