In this study, a numerical approach has been introduced in the elastic stress solutions of hollow disks made of functionally graded materials (FGM) that are exposed to linearly increasing temperature dispersion. The modulus of elasticity and the coefficient of thermal expansion of the FGM disk is assumed to vary in radial direction in different forms, and it is further assumed that the Poisson's ratio is constant. It causes the differential equation that manages the behaviour of the object under different material properties and boundary conditions to be a variable coefficient equation. Except for some simple grade materials and boundary conditions, it is hardly possible to produce an analytical solution of such equations. In this case, the solution of the problems can only be found with numerical approaches. Complementary Functions Method (CFM) was used to solve the problem. Different material models were used from the written works and corresponding radial, tangential and equivalent stresses and radial displacements were calculated. Simple, effective and well-structured solution steps can be easily implemented for disks.