A new analytical solution of the Fractional Neutron Point Kinetics Equations is developed in the present work, considering multi-groups of precursors of delayed neutrons as well as a constant reactivity. Such solution is obtained in a short way, as a direct generalization of the integer case, being only necessary to use an elementary change of variable as well as the Heaviside's Expansion Formula. One of the most relevant novelties of the developed solution is that it has the same structure and mathematical form of the solution of the integer case. This symmetry property is an important advantage over other reported solutions because it allows extending, in a straightforward way, computational algorithms that are used for the integer case to the fractional one. A MATLAB code implementation of the developed solution is also provided, as well as numerical experiments, where the proposed analytical solution is applied to several cases and scenarios, observing a sub-diffusive behavior for step reactivities, a mixed profile for ramp simulations, as well as a displacement of the peaks in the case of feedback reactivities.