A discrete ordinates code is developed with which to compute the beam spread function (BSF) without invoking the small-angle scattering approximation or performing Monte Carlo calculations. The computed BSF is used to predict the response of a detector versus its distance to the origin of a highly collimated beam, its angle with respect to the beam, and the two local angles that specify the detector orientation. Numerical results have been obtained for water models that simulate a clear ocean, a coastal ocean, and a turbid harbor. Six orders of magnitude or more change in the detector response caused by scattered photons can be predicted for different detector locations while simultaneously obtaining small changes for different detector orientations. This capability is useful for assessment of the sensitivity of the detector response to the interpretation of time-independent underwater imaging systems or visibility models.