"Problem 16.1.1" ket00 = (1,0,0,0) ket10 = (0,1,0,0) ket01 = (0,0,1,0) ket11 = (0,0,0,1) b1 = outer(ket00,ket10) - outer(ket01,ket11) b2 = outer(ket00,ket01) + outer(ket10,ket11) b1dag = transpose(b1) b2dag = transpose(b2) b1 b2 I = ((1,0,0,0),(0,1,0,0),(0,0,1,0),(0,0,0,1)) "Verify commutation relations" check(dot(b1dag,b1) + dot(b1,b1dag) == I) check(dot(b2dag,b2) + dot(b2,b2dag) == I) check(dot(b1dag,b2dag) + dot(b2dag,b1dag) == 0) check(dot(b1dag,b1dag) + dot(b1dag,b1dag) == 0) -- equation (16.22) check(dot(b1,b1) + dot(b1,b1) == 0) check(dot(b1,b2) + dot(b2,b1) == 0) check(dot(b2,b2) + dot(b2,b2) == 0) -- equation (16.14) check(dot(b1dag,b1dag) + dot(b1dag,b1dag) == 0) check(dot(b1dag,b2dag) + dot(b2dag,b1dag) == 0) check(dot(b2dag,b2dag) + dot(b2dag,b2dag) == 0) -- equation (16.32) check(dot(b1dag,b1) + dot(b1,b1dag) == I) check(dot(b1dag,b2) + dot(b2,b1dag) == 0) check(dot(b2dag,b1) + dot(b1,b2dag) == 0) check(dot(b2dag,b2) + dot(b2,b2dag) == I) "ok" Run