The Debye-Waller factor is calculated classically for the Mossbauer effect of an atom which is subject to a very anharmonic force from its neighbours in an otherwise harmonic lattice. The required integration over configuration space can be performed analytically over the coordinates of all atoms in the crystal but the active atom and its neighbours. The remaining integration, which, for a face-centered cubic lattice involves still 39 coordinates, is done by a Monte Carlo method. The numerical work assumes an FCC lattice with nearest-neighbour interactions, and a hard-core repulsion between the active atom and its neighbours, the core radius being one-half the cube edge. The Debye-Waller exponent, W, varies with wavevector k and temperature T approximately as k2(A+BT), A, B, being constants, and A by no means small. For fixed k, this resembles the behaviour predicted from simpler models, but no Einstein model with a spherically symmetric potential can reproduce the dependence on both variables, thus showing that it is essential to take into account the many-body aspects of the problem.