  We present a C program to compute by successive approximations the regular
order reduction of a large class of ordinary differential equations, which
includes evolution equations in electrodynamics and gravitation. The code may
also find the regular order reduction of delay-differential equations.
