We present a novel scheme for an unbiased, nonperturbative treatment of strongly correlated fermions. The proposed approach combines two of the most successful many-body methods, the dynamical mean field theory and the functional renormalization group. Physically, this allows for a systematic inclusion of nonlocal correlations via the functional renormalization group flow equations, after the local correlations are taken into account nonperturbatively by the dynamical mean field theory. To demonstrate the feasibility of the approach, we present numerical results for the two-dimensional Hubbard model at half filling.