  We implement an algorithm for the computation of Schouten bracket of weakly
nonlocal Hamiltonian operators in three different computer algebra systems:
Maple, Reduce and Mathematica. This class of Hamiltonian operators encompass
almost all the examples coming from the theory of (1+1)-integrable evolutionary
PDEs
