The exact expression for the maximum tension of a pendulum string is used to obtain a closed-form approximate expression for the solution of a simple pendulum in terms of elementary functions. This approximate solution has a rational harmonic expression and depends on an unknown function, which must be determined. This unknown function is expanded using the Padé approximant and two new parameters are introduced which are determined by means of a term-by-term comparison of the power series expansion for the approximate maximum tension with the corresponding series for the exact one. Using this approach, accurate approximate analytical expressions for the periodic solution are obtained. We also compared the Fourier series expansions of the approximate solutions and the exact ones. This allowed us to compare the coefficients for the different harmonics in these solutions. We also compared the approximate and exact solutions as a function of time for several oscillation amplitudes. Finally, in this procedure we used some of the approximate expressions for the simple pendulum frequency which can be found in the bibliography; however, the procedure can be applied using other approximate frequencies.