Density functional theory (DFT) has emerged as one of the most versatile and lucrative approaches in electronic structure calculations of many-electron systems in past four decades. Here we give an account of the development of a variational DFT method for atoms and molecules completely in a Cartesian grid. The non-relativistic Kohn-Sham equation is solved by using an LCAO-MO ansatz. Atom-centered localized basis set, electron density, molecular orbitals, twobody potentials are directly constructed on the grid. We adopt a Fourier convolution method for classical Coulomb potentials by making an Ewald-type decomposition technique in terms of shortand long-range interactions. It produces quite accurate and competitive results for various properties of interest, such as component energy, total energy, ionization energy, potential energy curve, atomization energy, etc. Both local and non-local functionals are employed for pseudopotential as well as full calculations. While most results are offered in a uniform grid, initial exploratory attempts are made in a non-uniform grid, which can significantly reduce the computational overhead.This offers a practical, viable alternative to atom-centered grid-based implementations, currently exploited by the majority of programs available world-wide. . † AKR dedicates this chapter to his kind-hearted uncle Dr. Pranjalendu Roy, on the occasion of his 67th birthday. Dr. Roy is a noted zoologist and social activist, who has inspired the author greatly. 1 3 . The second term corresponds to the external potential energy containing the nuclear-electron attraction, and the last expression refers to the classical electrostatic Hartree repulsion energy. While the sheer simplicity of replacing the troublesome manyelectron SE by a single equation in terms of ρ(r) alone is manifestly appealing, underlying