A first-order simulation algorithm for multivariable non-linear stochastic integral equations of the Ito type is derived, which utilizes the second-order properties of incremental Brownian motion. Various questions of stability and convergence of the algorithm are investigated. By way of application the algorithm is used in evaluating models and optimal control laws for biological watewater treatment plants of the activated sludge type.