By a novel exact diagonalization technique we show that bound pairs propagate
between repulsive Hubbard clusters in a superconducting fashion. The size of
the matrices that must be handled depends on the number of fermion
configurations {\em per spin}, which is of the order of the square root of the
overall size of the Hilbert space. We use CuO4 units connected by weak O-O
links to model interplanar coupling and c-axis superconductivity in Cuprates.
The numerical evidence on Cu2O8 and Cu3O12 prompts a new
analytic scheme describing the propagation of bound pairs and also the
superconducting flux quantization in a 3-d geometry.Comment: 5 pages, 3 figure