Computer simulations of a Boeing 727 class aircraft landing in turbulence were developed by programming the longitudinal aircraft equations of motion into a digital computer with various input values of vertical and horizontal wind speeds. Turbulent wind data was fed to the computer in one-second intervals. The computer computed in one-second intervals the aircraft speed, altitude, horizontal distance traveled, rate-of-descent, pitch attitude, glide path angle (from edge of runway) and elevator angle. All computer runs were made in the 'stick-fixed' mode. The RMS values of altitude and velocity perturbations (from equilibrium) were found to be large when horizontal wind gusts had sinusoidal components at or near the phugoid (long period) frequency. Maximum RMS altitude deviations occurred when the vertical wind had sinusoidal components which were 1/10 to 1/5 of the phugoid frequency. When real wind data (obtained from NCAR Queen Air) were used as input winds good correlations were found to exist between RMS velocity perturbations and both horizontal and vertical wind shears.