We have previously developed a grid-free approach to implement local density functionals based entirely on matrix representations of operators known from conventional electronic structure methods (Y.C. Zheng and J. Almlöf, Chem. Phys. Lett., 214 (1993) 397). In this article we present a scheme to implement non-local gradient corrected functionals and analytical energy derivatives of density functional theory. The results obtained show that our grid-free DFT approach is indeed a viable alternative to these grid-based methods.