By studying an appropriate Riemann-Hilbert (RH) problem, the inverse scattering of a higher-order nonlinear Schrödinger equation for inhomogeneous Heisenberg ferromagnetic system with zero boundary condition is calculated. The RH problem of reflection coefficient with multiple high-order poles is obtained. Meanwhile, the calculation formulas of bound state (BS) solitons and multiple BS solitons corresponding to one high-order pole and multiple high-order poles are also calculated respectively. Finally, the corresponding soliton solution model is calculated according to the corresponding formula. Simultaneously, we also obtain the interaction between BS soliton and multiple BS solitons, and between multiple BS solitons and multiple BS solitons. Keywords A higher-order nonlinear Schrödinger equation • Heisenberg ferromagnetism system • Riemann-Hilbert problem • Initial value problem • Multiple higher-order poles • Bound-state soliton solution.