An approximation method is developed which is applicable to the solving of integro-differential equations using a digital computer. This method is a generalization of the finite difference method in that information about the expected solution is added as trial solutions. This allows the possibility of significantly reducing the number of discrete mesh points necessary to achieve an accurate solution over that required by the conventional finite difference method. The method is demonstrated on some space-time reactor kinetics problems.