On the time-fractional Schr¨odinger equation: theoretical analysis and numerical solution by matrix Mittag–Leffler functions✩