High-order schemes for neutron kinetics calculations, based on a local polynomial approximation