A method is described for the computation of the back-scattered intensities of atomic beams, diffracted from the evanescent field generated outside an optical plate by internal counter-propagating laser beams. The method derives from a procedure developed for the similar but importantly differing problem of Low Energy Electron Diffraction, (Lynch and Smith, 1983). Modifications to that theory required for the present problem bring the equations closer to the RHEED solution proposed by Ichimiya (1983). Results from multi-slicing from the same narrow field depth (2 Aangstroems) in this strong field case show the success and also limitations of the program in its present form. Computation to greater depth in the strong field leads to numerical instabilities due to the incorporation of very large tunnelling terms. This requires the application of boundary conditions at each slice rather than the end of the multi-slice calculation. 7 refs., 3 figs.