This paper deals with the numerical solutions of the nonlinear heat transfer equation with inhomogeneous Neumann boundary conditions on a semiinfinite interval.Laguerre-Gauss-Radua nodes are used to construct the Nthdegree Lagrange interpolation function that approximates the solution to the nonlinear heat transfer equation on a semi-infinite interval.Efficient algorithms are implemented.Numerical results demonstrate efficiency and high accuracy of this approach.Especially,it is much easier to deal with the nonlinear heat transfer on a semi-infinite interval.The proposed method is also applicable to other nonUnear problems defined on certain unbounded domains.