from sympy import *
init_printing()

from sympy.physics.quantum import *
from sympy.physics.quantum.boson import *
from sympy.physics.quantum.fermion import *
from sympy.physics.quantum.operatorordering import *

a, b = BosonOp("a"), BosonOp("b")
c, d = FermionOp("c"), FermionOp("d")

Commutator(a, Dagger(a)).doit()

AntiCommutator(a, Dagger(a)).doit()

normal_ordered_form(AntiCommutator(a, Dagger(a)).doit())

Commutator(a, b).doit()

AntiCommutator(a, b).doit()

Commutator(a, b).doit(independent=True)

AntiCommutator(a, b).doit(independent=True)

AntiCommutator(c, Dagger(c)).doit()

Commutator(c, Dagger(c)).doit().expand()

Commutator(c, d).doit()

AntiCommutator(c, d).doit()

Commutator(c, d).doit(independent=True)

AntiCommutator(c, d).doit(independent=True)

Commutator(a, c).doit()

AntiCommutator(a, c).doit()

%reload_ext version_information

%version_information sympy