Copyright © 2011 Mohammed Yusuf Waziri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The problems of radiative transfer give rise to interesting integral equations that must be faced with efficient numerical solver. Very often the integral equations are discretized to large-scale nonlinear equations and solved by Newton's-like methods. Generally, these kind of methods require the computation and storage of the Jacobian matrix or its approximation. In this paper, we present a new approach that was based on approximating the Jacobian inverse into a diagonal matrix by means of variational technique. Numerical results on well-known benchmarks integral equations involved in the radiative transfer authenticate the reliability and efficiency of the approach. The fact that the proposed method can solve the integral equations without function derivative and matrix storage can be considered as a clear advantage over some other variants of Newton's method.