The Fejer [1955] pulse interaction experiment is simulated numerically by solving, in an iterative fashion, the electron energy balance and magnetoionic propagation equations over a height-time "matrix" representing a disturbed ionosphere. The complex interdependence of physical variables is more accurately represented by including "selfdistortion" effects on the disturbing pulse, and allowing electron collision frequency and the energy loss coefficient (G) to vary with electron energy and time. Provision is made for nonrectangular disturbing pulses as well as finite-length probing pulses. Phase as well as amplitude interaction effects are calculated. Absorption and phase shift computations are based on the Appleton-Hartree or Sen-Wyller complete (arbitrary angle of propagation) formulations, but the insertion of more sophisticated developments in interaction theory, including electron density perturbations, should be possible without appreciable increase in program complexity. Representative calculations demonstrate the effects of variations in experimental parameters , perturbations of ionospheric conditions, and various interpretational assumptions. The numerical approach offers advantages of simplicity and flexibility to those engaged in design and interpretation of wave interaction experiments. Synthesis of D-region electron density and collision frequency profiles, using the technique "in reverse," should prove superior to presently used trial-and-error methods.