Approximation of the solution of certain nonlinear ODEs with linear complexity

EZEQUIEL DRATMAN

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

ELSEVIER SCIENCE BV

Lugar: Amsterdam; Año: 2010 vol. 233 p. 2339 - 2350

0377-0427

We study the positive stationary solutions of a standard finite-difference discretization of the semilinear heat equation with nonlinear Neumann boundary conditions. We prove that there exists a unique solution of such a discretization, which approximates the unique positive stationary solution of the "continuous´´ equation. Furthermore, in this case we exhibit an algorithm computing an epsilon-approximation of such a solution by means of a homotopy continuation method. The cost of our algorithm is linear in the number of nodes involved in the discretization and the logarithm of the number of digits of approximation required.