s3o6=sqrt(3)/6
a=[1/4 1/4-s3o6; 1/4+s3o6 1/4]
b=[1/2,1/2]
c=[1/2-s3o6,1/2+s3o6]
fm(k,l)=b[k]*a[k,l]+b[l]*a[l,k]-b[k]*b[k]