  We use a numerical method, the finite-mode approach, to study inhomogeneous
condensation in effective models for QCD in a general framework. Former
limitations of considering a specific ansatz for the spatial dependence of the
condensate are overcome. Different error sources are analyzed and strategies to
minimize or eliminate them are outlined. The analytically known results for
$1+1$ dimensional models (such as the Gross-Neveu model and extensions of it)
are correctly reproduced using the finite-mode approach. Moreover, the NJL
model in $3+1$ dimensions is investigated and its phase diagram is determined
with particular focus on the inhomogeneous phase at high density.
