I am trying to define some AntiCommutator-Relation to work with operators depending on a large number of creation- and annihilation-operators in the fermionic case.

Let "b" be the creation and "a" be the annihilation-operator. Then i am looking, for example, for "a[i]a[j]b[k]a[l]a[m]"...

Now I define the relations:

################

restart

with(Physics)

Setup(quantumop={a,b},algebrarule={%AntiCommutator(a[i],a[j])=KroneckerDelta[i,j...