Considers how, in a reasonably motivated manner, to define and then evaluate unique finite differences between individually divergent sums S= Sigma infinity f(n) and integrals I= integral infinity dnf(n). This is done by first replacing f(n) by f(n mod lambda )=f(n)g(n mod lambda ), where g is a cutoff function (g(n to infinity mod lambda )=0) which obeys the permanence conditions g(n mod lambda to infinity )=1, and then transforming (S-I) into some convenient explicit functional D(f(n mod lambda )) which admits limlambda to infinity D(f(n mod lambda ))=D(limlambda to infinity f(n mod lambda ))=D(f(n)), where the final expression converges. There ought but there seems not to exist a convenient yet reasonably general theory for identifying admissible classes of summands and cutoffs and for deriving D(f), even though physicists have long dealt with simple cases ad hoc. Three specialised prescriptions are presented for D(f); one for suitably analytic f and g; another for merely differentiable f and g; and the unconventional ' epsilon -averaging method'. The mutual compatibility of these methods is discussed. Explicit differences are worked out for logarithmic, power-law, and exponential summands.