We consider an optimal control problem (OCP) for a partial differential equation (PDE) with random coefficients. The function is deterministic, distributed forcing term that minimizes expected quadratic regularized loss functional. For the numerical approximation of this PDE-constrained OCP, we replace expectation in objective functional by suitable quadrature formula and, eventually, discretiz...