This article investigates a family of approximate solutions for the fractional model (in the Liouville-Caputo sense) of the Ebola virus via an accurate numerical procedure (Chebyshev spectral collocation method). We reduce the proposed epidemiological model to a system of algebraic equations with the help of the properties of the Chebyshev polynomials of the third kind. Some theorems about the convergence analysis and the existence-uniqueness solution are stated. Finally, some numerical simulations are presented for different values of the fractional-order and the other parameters involved in the coefficients. We also note that we can apply the proposed method to solve other models.