A thermal buckling analysis of functionally graded thick rectangular plates according to von Karman non-linear theory is presented. The material properties of the functionally graded plate, except for the Poisson's ratio, were assumed to be graded in the thickness direction, according to a power-law distribution, in terms of the volume fractions of the metal and ceramic constituents. Formulations of equilibrium and stability equations are derived using the high order shear deformation theory based on different types of shape functions. Analytical method for determination of the critical buckling temperature for uniform increase of temperature, linear and non-linear change of temperature across thickness of a plate is developed. Numerical results were obtained in МATLAB software using combinations of symbolic and numeric values. The paper presents comparative results of critical buckling temperature for different types of shape functions. The accuracy of the formulation presented is verified by comparing to results available from the literature.