"Spin flip" B = B0 cos(omega t) (0,0,1) Sx = 1/2 hbar ((0,1),(1,0)) Sy = 1/2 hbar ((0,-i),(i,0)) Sz = 1/2 hbar ((1,0),(0,-1)) S = (Sx,Sy,Sz) H = e / m dot(B,S) H "Verify equation (1)" c1 = exp(-i e / (2 m omega) B0 sin(omega t)) / sqrt(2) c2 = exp(i e / (2 m omega) B0 sin(omega t)) / sqrt(2) s = (c1,c2) check(i hbar d(s,t) == dot(H,s)) "ok" "Verify equation (2)" check(dot(conj(s),Sx,s) == 1/2 hbar cos(e / (m omega) B0 sin(omega t))) "ok" Run