  We discuss a numerical algorithm for solving nonlinear integro-differential
equations, and illustrate our findings for the particular case of Volterra type
equations. The algorithm combines a perturbation approach meant to render a
linearized version of the problem and a spectral method where unknown functions
are expanded in terms of Chebyshev polynomials (El-gendi's method). This
approach is shown to be suitable for the calculation of two-point Green
functions required in next to leading order studies of time-dependent quantum
field theory.
