We discuss the possible extension of the bosonic classical field theory simulations to include fermions. This problem has been addressed in terms of the inhomogeneous mean-field approximation by Aarts and Smit. By performing a stochastic integration of an equivalent set of equations we can extend the original 1+1 dimensional calculations so that they become feasible in higher dimensions. We test the scheme in 2+1 dimensions and discuss some classical applications with fermions for the first time, such as the decay of oscillons.