J. Henderson et Hb. Thompson, Difference equations associated with fully nonlinear boundary value problems for second order ordinary differential equations, J DIF EQ AP, 7(2), 2001, pp. 297-321
We study the continuous problem y"=f(x,y,y'), xc[0,1], 0=G((y(0),y(1)),(y'(
0), y'(1))), and its discrete approximation (y(k+1)-2y(k)+y(k-1))/h(2) =f(t
(k), y(k), v(k)), k = 1,..., n-1, 0 = G((y(0), y(n)), (v(1), v(n))), where
f and G = (g(0), g(1)) are continuous and fully nonlinear, h = 1/n, v(k) =
(y(k) - y(k-1))/h, for k =1,..., n, and t(k) = kh, for k = 0,...,n. We assu
me there exist strict lower and strict upper solutions and impose additiona
l conditions on f and G which are known to yield a priori bounds on, and to
guarantee the existence of solutions of the continuous problem. We show th
at the discrete approximation also has solutions which approximate solution
s of the continuous problem and converge to the solution of the continuous
problem when it is unique, as the grid size goes to 0. Homotopy methods can
be used to compute the solution of the discrete approximation. Our results
were motivated by those of Gaines.