Many extensive groundwater systems comprise several aquifers hydraulically connected to one another through layers which are usually semiconfining. The aquifers respond conjunctively to stresses imposed on the system. Two types of models are developed in this report. The first is a digital computer model simulating groundwater flow in multiaquifer systems. The aquifers can be isotropic or anisotropic with boundary conditions of the first kind (Dirichlet type) or of the second kind (Neumann type). The second are analytical models of one-dimensional flow in systems comprising three interacting aquifers which represent idealized textbook type artesian basins. The analytical solutions are programmed on a digital computer for evaluation and some typical results, using parameters selected from the Roswell Basin, New Mexico, are presented.