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]);
> Np(3,2,psi);
> site:=1;
> spinMtest:=array(1..dimH);
>
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);
>
Sy(site,Sz(site,psi))-Sz(site,Sy(site,psi));
I*Sx(site,psi);
Test passed!!!
>