We consider a data-driven method, which combines Koopman operator theory with Extended Dynamic Mode Decomposition. We apply this method to the hypergeometric equation which is the Fuchsian equation with three regular singular points. The space of solutions at any of its singular points is a two-dimensional linear vector space on the field of reals when the independent variable is restricted to take values in the real axis and the unknown function is restricted to be a real-valued function of a real variable. A basis of the linear vector space of solutions is spanned by the hypergeometric function and its products with appropriate powers of the independent variable or the logarithmic function depending on the roots of the indicial equation of the hypergeometric equation. With our work, we obtain a new representation of the fundamental solutions of the hypergeometric equation and relate them to the spectral analysis of the finite approximation of the Koopman operator associated with the hypergeometric equation. We expect that the usefulness of our results will come more to the fore when we extend our study into the complex domain.