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.