The Kuhfuß theory of convection is derived from the exact hydrodynamic equations by means of anelastic approximation and diffusion approximation. In the local and time-independent limit the result is identical to the MLT case. In the more general case the theory includes overshooting and convective mixing as well as time-dependence. We present test calculations in a simple box model and first steps toward the implementation in a stellar code.