We develop a computational approximation to the intensity of a Gibbs spatial point process having interactions of any order. Limit theorems from stochastic geometry, and small-sample probabilities estimated once and for all by an extensive simulation study, are combined with scaling properties to form an approximation to the moment generating function of the sufficient statistic under a Poisson process. The approximate intensity is obtained as the solution of a self-consistency equation.