  Poisson brackets between conserved quantities are a fundamental tool in the
theory of integrable systems. The subclass of weakly nonlocal Poisson brackets
occurs in many significant integrable systems. Proving that a weakly nonlocal
differential operator defines a Poisson bracket can be challenging. We propose
a computational approach to this problem through the identification of such
operators with superfunctions on supermanifolds.
