We study the effects of Coulomb repulsive interactions on a Majorana Benalcazar–Bernevig–Huges (MBBH) model. The MBBH model belongs to the class of second-order topological superconductors (HOTSC2), featuring robust Majorana corner modes. We consider an interacting strip of four chains of length L and perform a density matrix renormalization group (DMRG) numerical simulation based on a tensor-network approach. Study of the non-local fermionic correlations and the degenerate entanglement spectrum indicates that the topological phases are robust in the presence of interactions, even in the strongly interacting regime.