Two or more distinct materials are combined into a single functionally graded material (FGM) where the microstructural composition and properties change gradually. Thermal post-buckling behavior of uniform slender FGM beams is investigated independently using the classical Rayleigh-Ritz (RR) formulation and the versatile Finite Element Analysis (FEA) formulation developed in this paper. The von-Karman strain-displacement relations are used to account for moderately large deflections of FGM beams. Bending-extension coupling arising due to heterogeneity of material through the thickness is included. Simply supported and clamped beams with axially immovable ends are considered in the present study. Post-buckling load versus deflection curves and buckled mode shapes obtained from both the RR and FEA formulations for different volume fraction exponents show an excellent agreement with the available literature results for simply supported ends. Response of the FGM beam with clamped ends is studied for the first time and the results from both the RR and FEA formulations show a very good agreement. Though the response of the FGM beam could have been studied more accurately by FEA formulation alone, the authors aim to apply the RR formulation is to find an approximate closed form post-buckling solutions for the FGM beams. Further, the use of the RR formulation clearly demonstrates the effect of bending-extension coupling on the post-buckling response of the FGM beams.