A method to solve the Liouville equation for an ensemble of two‐dimensional viscous flows that are driven by random forcing with arbitrary statistics is outlined. By appropriate transformations of both the dependent variable (probability distribution) and independent variables, and by expansion of the solution in eigenfunctions of the separable part of the Liouville operator, it is found possible to reduce the problem to that of solving a simultaneous system of nonhomogeneous linear algebraic equations. The equilibrium kinetic energy and energy‐transfer spectra can be calculated directly from the equilibrium probability distribution.