Many of the equations encountered in fluid dynamics and plasma physics are cubically nonlinear. The properties of a relatively simple statistical theory of cubically nonlinear equations, analogous to Kraichnan’s direct interaction approximation for quadratically nonlinear equations, are explored. Its predictions are compared with the numerical solution of a set of equations called the random‐coupling model; good agreement is found.