"Spin part 2" s = (cos(theta/2), sin(theta/2) exp(i phi)) theta = 1/3 pi s "Verify expectation value" Sz = 1/2 hbar ((1,0),(0,-1)) dot(conj(s),Sz,s) check(last == 1/4 hbar) "ok" Run