The Dirac equation is studied for a sufficiently large class of Lemaįtre—Tolman—Bondi cosmological models. While the angular equation (whose solution is known) separates directly, the spatial and temporal dependence de-couples only after a suitable separation procedure. The separated time equation is integrated by series. The separated spatial equation still depends on an arbitrary function relative to the integration of the cosmological model.