The magnetization profile and the related magnetic small-angle neutron scattering cross section of a single spherical nanoparticle with Néel surface anisotropy are analytically investigated. A Hamiltonian is employed that comprises the isotropic exchange interaction, an external magnetic field, a uniaxial magnetocrystalline anisotropy in the core of the particle and the Néel anisotropy at the surface. Using a perturbation approach, the determination of the magnetization profile can be reduced to a Helmholtz equation with Neumann boundary condition, whose solution is represented by an infinite series in terms of spherical harmonics and spherical Bessel functions. From the resulting infinite series expansion, the Fourier transform, which is algebraically related to the magnetic small-angle neutron scattering cross section, is analytically calculated. The approximate analytical solution for the spin structure is compared with the numerical solution using the Landau–Lifshitz equation, which accounts for the full nonlinearity of the problem. The signature of the Néel surface anisotropy can be identified in the magnetic neutron scattering observables, but its effect is relatively small, even for large values of the surface anisotropy constant.