We introduce the notion of two-scale convergence for partial differential equations with random coefficients that gives a very efficient way of finding homogenized differential equations with random coefficients. For an application, we find the homogenized matrices for linear second order elliptic equations with random coefficients. We suggest a natural way of finding the two-scale limit of second order equations by considering the flux term.