Abstract In this paper, the nonlinear Volterra series expansion is extended and used to describe certain types of nonautonomous differential equations related to the inverse scattering problem in nuclear physics. The nonautonomous Volterra series expansion lets us determine a dynamic, polynomial approximation of the variable phase approximation (VPA), which is used to determine the phase shifts from nuclear potentials through first-order nonlinear differential equations. By using the first-order Volterra expansion, a robust approximation is formulated to the inverse scattering problem for weak potentials and/or high energies. The method is then extended with the help of radial basis function neural networks by applying a nonlinear transformation on the measured phase shifts to be able to model the scattering system with a linear approximation given by the first-order Volterra expansion. The method is applied to describe the ${}^1S_0$ NN potentials in neutron+proton scattering below 200 MeV laboratory kinetic energies, giving physically sensible potentials and below $1\%$ averaged relative error between the recalculated and the measured phase shifts.