[en] Inhomogeneity terms in the expansion of the kinetic energy density are included and the Euler-Lagrange equations solved. Shell effects may be incorporated in a simple way. The study of spherical shapes of large systems is given as an illustration of the method proposed