%load_ext sympyprinting

from sympy import sqrt, symbols, Rational
from sympy import expand, Eq, Symbol, simplify, exp, sin
from sympy.physics.quantum import *
from sympy.physics.quantum.qubit import *
from sympy.physics.quantum.gate import *
from sympy.physics.quantum.grover import *
from sympy.physics.quantum.qft import QFT, IQFT, Fourier
from sympy.physics.quantum.circuitplot import circuit_plot

phi, psi = Ket('phi'), Ket('psi')
alpha, beta = symbols('alpha beta', complex=True)

state = alpha*psi + beta*phi; state

ip = Dagger(state)*state; ip

qapply(expand(ip))

A = Operator('A')
B = Operator('B')
C = Operator('C')

A*B == B*A

expand((A+B)**2)

comm = Commutator(A,B); comm

comm.doit()

comm = Commutator(A*B,B+C); comm

comm.expand(commutator=True)

_.doit().expand()

Dagger(_)