The simplified spherical harmonics(SPN) method is utilized to discretize the angular variables of neutron transport time-space kinetics equation.The finite element method,fully implicit scheme and direct analytical time integration method are used to deal with the spatial,time variables and delayed neutron precursor equation,respectively.According to the model,a computer program is developed to solve multi-dimensional time-space kinetics equation in unstructured-meshes.The numerical results show that this method can be used to perform the kinetics calculation in complex environment.