This paper deals with the numerical solution of the nonlinear heat transfer problem in a multilayered plate. Kansa's meshless method is used for the solution of this problem. In this approach, the unknown temperatures in layers are approximated by the linear combination of radial basis functions, while the governing equation and the boundary conditions are imposed directly at the collocation points. The multiquadrics [MQ] are used as the radial basis functions. In the presented method the radial basis functions contains a free parameter C, called the shape parameter. Usually, in the application of radial basis functions, this parameter is chosen arbitrarily depending on the author's experience. In the presented paper, special attention is paid to the optimal choice of the shape parameter for the radial basis functions. This optimal value of the shape parameter is obtained using a formula given by other authors for solution of the linear case.