ABSTRACT

The last two terms in Eq. (12.109) involve the derivative at the boundary points. For Dirichlet boundary conditions, these terms are not needed. For Neumann boundary conditions, these two terms introduce the derivative boundary conditions at the boundaries of the global solution domain. Substituting Eq. (12.109) into Eq. (12.108) yields

/(y(x)) = J:(-y'~ + Qyff.; - Fff.;)dx +Ybff.;(b) - Y~ff.;(a) = 0 (12.110) In terms of the global approximate solution y(x) and the discretized global solution

domain illustrated in Figure 12.4, Eq. (12.110) can be written as follows: /(y(x)) = /(l)(y(x)) + /(2)(y(x)) + ... + /(i-I)(y(X)) + /(i)(y(x)) + ...