epsilon1 = ((0,0,0),(0,0,1),(0,-1,0)) epsilon2 = ((0,0,-1),(0,0,0),(1,0,0)) epsilon3 = ((0,1,0),(-1,0,0),(0,0,0)) Jx = -i epsilon1 Jy = -i epsilon2 Jz = -i epsilon3 U = ((-1/sqrt(2),i/sqrt(2),0),(0,0,1),(1/sqrt(2),i/sqrt(2),0)) Udag = transpose(conj(U)) "Verify U is unitary" check(inv(U) == Udag) "ok" Jxprime = dot(U,Jx,Udag) Jyprime = dot(U,Jy,Udag) Jzprime = dot(U,Jz,Udag) Jxprime Jyprime Jzprime Run