In this study, we propose an extension of the modified Newton-Househölder methods to find multiple roots with unknown multiplicity of nonlinear equations. With four functional evaluations per iteration, the proposed method achieves an optimal eighth order of convergence. The higher the convergence order, the quicker we get to the root with a high accuracy. The numerical examples have shown that this scheme can compete with the existing methods. This scheme is also stable across all of the functions tested based on the graphical basins of attraction.
