A finite-difference method is presented for solving pseudo-two-dimensional boundary-value problems. The sparse and nearly block tridiagonal properties of the matrices generated by using the finite-difference method for problems of this type are fully utilized and maintained, which yields a method that is highly efficient in the use of storage space and computation. An example shows that the central process unit time required by the method is significantly less than that required by an alternative method.