  Based on the non-local light-ray operator technique, we develop an algorithm
for the computational calculation of evolution kernels for higher-twist
operators in leading order of the perturbation theory. We compute the evolution
kernel for the twist-3 operators in the flavour singlet channel. Our result
confirms the local anomalous dimensions computed by Bukhvostov, Kuraev, and
Lipatov as well as the non-local evolution kernels of Balitzky and Braun.
