Sloppily transpiled from some python code written for a paper that used some wick's theorem. Here for posterity and in case someone else might use it.

Writeup on wick's theorem at some point. For now look here: https://en.wikipedia.org/wiki/Second_quantization#Creation_and_annihilation_operators and here https://en.wikipedia.org/wiki/Normal_order.

More iffy on the bosonic results since I did not test it on those.