+ + +

Subroutine pdel presented in Section 9.13.1 can be modified to include the derivative BC on the right-hand side simply by solving Eq. (9.121) for f(imax,j) (j = 2,jmax - 1). Similar procedures can be implemented on the other three sides of a rectangular domain.