The problem of a strongly driven one‐dimensional electrostatic plasma is treated by the method of expansion of the distribution function in terms of Hermite polynomials in velocity space. Numerical results for hydrodynamic variables are similar to those found in simulation calculations. The shapes of the ion and electron distribution functions are significantly altered from Maxwellian form during heating.