Simulation by Empirical Orthogonal Functions