Originally written to solve systems of integrodifferential equations via collocation method in an arbitrary number of variables. The current implementation of the method is based on a finite-dimensional orthogonal system of functions.

Code to solve the non-linear Balitsky-Kovchegov evolution equation in QCD