  We consider tilings of a rectangle which is n units wide and m units long by
non-overlapping 1 X 1 squares and s X s squares. Bivariate generating functions
are computed with the Transfer Matrix Method for moderately large but fixed
widths n as a function of the parameter m and of the number of s X s squares in
the rectangle.
