A method for solving the inverse Sturm–Liouville problem on a finite interval is proposed. It is based on a Fourier–Legendre series representation of the integral transmutation kernel. Substitution of the representation into the Gel’fand–Levitan equation leads to a linear algebraic system of equations and consequently to a simple algorithm for recovering the potential. Numerical illustrations are presented.