Test of Spin operators

Let us check commutation relation [Sx(k),Sy(k)]=iSz(k)

> psi:=vec([1,0,0,0],[0,0,1,1]);

psi := vec([1, 0, 0, 0],[0, 0, 1, 1])

> Np(3,2,psi);

vec([1, 0, 0, 0],[0, 0, 1, 1])

> site:=1;

site := 1

> spinMtest:=array(1..dimH);

spinMtest := array(1 .. 36,[])

> for i from 1 to dimH do
spinMtest[i]:=Sy(site,Sz(site,Basis[i]))-Sz(site,Sy(site,Basis[i]))-I*Sx(site,Basis[i]);
od:
print(spinMtest);

vector([0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0...

> Sy(site,Sz(site,psi))-Sz(site,Sy(site,psi));
I*Sx(site,psi);

1/2*I*vec([0, 0, 0, 0],[1, 0, 1, 1])

1/2*I*vec([0, 0, 0, 0],[1, 0, 1, 1])

Test passed!!!

>