We review the derivation of a phase field crystal (PFC) model from classical density functional theory (DFT). Through a gradient flow of the Helmholtz free energy functional and appropriate approximations of the correlation functions, higher order nonlinear equations are derived for the evolution of a time averaged density. The equation is solved by finite elements using a semi-implicit time discretization.