adding stochastic component to Gordon-Schaefer model