We propose a direct solver for the three-dimensional Poisson equation with a variable coefficient, and an algorithm to directly solve the associated sparse linear systems that exploits the sparsity pattern of the coefficient matrix. Introducing some appropriate finite difference operators, we derive a second-order scheme for the solver, and then two suitable high-order compact schemes are also discussed. For a cube containing N nodes, the solver requires arithmetic operations and memory to store the necessary information. Its efficiency is illustrated with examples, and the numerical results are analysed.

A general birth-and-death queueing model is considered with N waiting positions (0 ≦ N ≦ ∞), s servers (1 ≦ s ≦ ∞) and a first-come-first-served queueing discipline. The first and second order moments of the steady-state waiting-time (excluding service) for a non-lost arriving customer are given. By setting the busy period equal to the time where at least one service is in progress, we obtain the first and second order moments of the length of a busy period and also the distribution of the number served during it, given an arbitrary number of customers present originally. Using a direct approach all expressions are given in explicit forms which, although being far from elegant, are suitable for evaluation on a computer.