"Dirac equation 1" alphax = ((0,0,0,1),(0,0,1,0),(0,1,0,0),(1,0,0,0)) alphay = ((0,0,0,-i),(0,0,i,0),(0,-i,0,0),(i,0,0,0)) alphaz = ((0,0,1,0),(0,0,0,-1),(1,0,0,0),(0,-1,0,0)) beta = ((1,0,0,0),(0,1,0,0),(0,0,-1,0),(0,0,0,-1)) gammat = ((1,0,0,0),(0,1,0,0),(0,0,-1,0),(0,0,0,-1)) gammax = ((0,0,0,1),(0,0,1,0),(0,-1,0,0),(-1,0,0,0)) gammay = ((0,0,0,-i),(0,0,i,0),(0,i,0,0),(-i,0,0,0)) gammaz = ((0,0,1,0),(0,0,0,-1),(-1,0,0,0),(0,1,0,0)) check(dot(beta,alphax) == gammax) check(dot(beta,alphay) == gammay) check(dot(beta,alphaz) == gammaz) "ok" Run