We will perform a detailed study of the matter-ekpyrotic bouncing scenario in Loop Quantum Cosmology using the methods of the dynamical systems theory. We will show that when the background is driven by a single scalar field, at very late times, in the contracting phase, all orbits depict a matter dominated Universe, which evolves to an ekpyrotic phase. After the bounce the Universe enters in the expanding phase, where the orbits leave the ekpyrotic regime going to a kination (also named deflationary) regime. Moreover, this scenario supports the production of heavy massive particles conformally coupled with gravity, which reheats the universe at temperatures compatible with the nucleosynthesis bounds and also the production of massless particles non-conformally coupled with gravity leading to very high reheating temperatures but ensuring the nucleosynthesis success. Dealing with cosmological perturbations, these background dynamics produce a nearly scale invariant power spectrum for the modes that leave the Hubble radius, in the contracting phase, when the Universe is quasi-matter dominated, whose spectral index and corresponding running is compatible with the recent experimental data obtained by PLANCK's team. Big Bang singularity: for example violating the null energy condition in General Relativity by incorporating new forms of matter such as phantom [2] or quintom fields [3], Galileons [4] or phantom condensates [5], or by adding terms to Einstein-Hilbert action [6], but the simplest one is to go beyond General Relativity and consider holonomy corrected Loop Quantum Cosmology (LQC), where a Big Bounce replaces the Big Bang singularity [7]. In fact, other future singularities such as Type I (Big Rip) and Type III (Big Freeze) are also forbbiden in holonomy corrected LQC [8].On the other hand, it is well known that a matter domination period in the contracting phase is dual to the de Sitter regime in the expanding phase [9], which provides a flat power spectrum of cosmological perturbations in the matter bounce scenario. Moreover, an abrupt ekpyrotic phase transition is needed, in the contracting phase, in order to solve the famous Belinsky-Khalatnikov-Lifshitz (BKL) instability: the effective energy density of primordial anisotropy scales as a −6 in the contracting phase [10] and, more important, to produce enough particles which will be responsible for thermalizing the universe in the expanding phase [11,12]. This new scenario named matter-ekpyrotic scenario was introduced in [13], being developed within the two-field model in [14], while in [15] there is a numerical discussion on the primordial anisotropy issue (see [16] for a review). It was treated in the context of Loop Quantum Cosmology in [17], showing that it could be compatible with the new experimental data [18] provided by PLANCK's team (see also [19]).The main goal of the present work is to study mathematically, from the viewpoint of the dynamical systems theory, the matter-ekpyrotic scenario in the context of holonomy corrected LQC. More precisely...