In this paper, an analytical method and a partial differential equation based solution to control temperature distribution for functionally graded (FG) plates is introduced. For the rectangular FG plate under consideration, it is assumed that the material properties such as thermal conductivity, density, and specific heat capacity vary in the width direction, and the governing heat conduction equation of the plate is a second-order partial differential equation. Using Lyapunov’s theorem, it is shown that by applying controlled heat flux through the boundary of the domain, the temperature distribution of the plate will approach a desired steady-state distribution. Numerical simulation is provided to verify the effectiveness of the proposed method such that by applying the boundary transient heat flux, in-domain distributed temperature converges to its desired steady-state temperature.