  A program dedicated to the numerical solution of the evolution equations for
twist-three multiparton correlation functions is presented. The solutions are
obtained by direct integration on a discretized momentum fraction grid. Both
flavor nonsinglet and flavor singlet evolution (in both $C$-parity sectors) can
be addressed. Physical applications include single spin asymmetries and the
subleading twist contribution to the polarized structure function $g_2$. An
arbitrary input is accepted for the initial distributions.
