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