This study aims to present an accelerated derivative-free method for solving systems of nonlinear equations using a double direction approach. The approach approximates the Jacobian using a suitably formed diag-onal matrix by applying the acceleration parameter. Moreover, a norm descent line search is employed in the scheme to compute the optimal step length. Under the primary conditions, the proposed method’s global con-vergence is proved. Numerical results are recorded in this paper using a set of large-scale test problems. Moreover, the new method is successfully used to address the problem of Chandrasekhar’s integral equation problem appearing in radiative heat transfer. This method outperforms the existing Newton and inexact double step length methods.