The distribution of neutron population in nuclear reactor is described by using transport equations. One of possible approximations of neutron transport equation is given by the neutron diffusion equation. The paper presents numerical solution method of one group neutron diffusion equation with one group of delayed neutrons.