An approach to numerically integrate the Landau-Lifshitz fluctuating hydrodynamic equations is outlined. The method is applied to one-dimensional systems obeying the nonlinear Fourier equation and the full hydrodynamic equations for a dilute gas. Static spatial correlation functions are obtained from computer-generated sample trajectories (time series). They are found to show the emergence of long-range behavior whenever a temperature gradient is applied. The results are in very good agreement with those obtained from solving the correlation equations directly.