摘要

By exploiting the implicit structure of Chandrasekhar H-equation, we devise a structured Newton's method to compute the solution, with physical meaning. We also consider the critical case of this equation, in which the Jacobian matrix at the solution is singular. By incorporating a shift technique, we show that our structured Newton's method can still converge quadratically to the solution. Numerical experiments are listed to indicate that the new method is a feasible and effective solver for Chandrasekhar H-equation.