In this paper, a nonstandard finite difference scheme is employed to approximate the solution of a nonlinear second order differential equation obtained from the model of heat transfer in extended surfaces. The proposed method for the system enjoys the nonlocal approximation of nonlinear terms and requires neither any iterative procedure nor the computation of Jacobian unlike the standard finite difference method. Accuracy of the method and conditions for generating expected positive solution are also discussed in detail.