POMA is a very handy tool for Product Operator Formalism calculation for NMR pulse program development.
It works with Mathematica, and various operator spaces (Ix, Iy, Iz,
or I±, Iz) are available.
The original web page of
POMA in ETH by Peter Günter
You can also find useful description and POMA examples of Product Operator Fomalism in
Columbia U.