The temporal subdomain method based on the Ritz-Galerkin method is investigated as a method for the solution of space-time dependent neutron dynamics equations. In the temporal subdomain finite element method, the time domain is divided into subdomains and within each subdomain the unknown coefficients of the time dependent trial functions are determined by making the residual of an appropriate functional orthogonal to the step function.