Implement the Virasoro Algebra and it's representations. The Virasoro Algebra is an infinite dimensional Lie algebra with generators L_n, n::integer satisfying [L_m,L_n]=(m-n)*L_(m+n)+c/12*m*(m^2-1)*KroneckerDelta[m,-n] This is quite easy to implement with the "Physics" package: with(Physics); Setup(quantumop={L},algebrarule={ %Commutator=(m-n)*L[m+n]+c/12*m*(m^2-1)*KroneckerDelta[m,-n]});