#ES(): the four entangled mixed states Sing,T1,T2,T3
ES:=proc()
[[0,1,-1,0]/sqrt(2),
 [0,1,1,0]/sqrt(2),
 [1,0,0,1]/sqrt(2),
 [1,0,0,-1]/sqrt(2)
]:
end: