A generation matrix theory of full-sib mating is developed in which 13 mating "classes" are distinguished according to identity of genes in individuals mated and identity of genotypes as belonging to homozygous, parental, or offspring sets. The 13 times 13 matrix reveals some properties of the full-sib mating system not shown by previous work. The eigenvalues and a set of eigenvectors for the generation matrix, and the general solution for the frequencies of mating classes among descendants of an original mating of genotypes ab times cd, are given. The genotypic array of descendants in an arbitrary generation is also given. A new formula is derived for the coefficient of inbreeding in generation n + m in terms of coefficients of inbreeding in earlier generations. An algorithm is presented for calculating the probability of a given situation of identity of alleles carried by two individuals given only the indices of their own respective generations and the generation of their most recent common ancestor. The application of such probabilities to obtaining covariances between relatives in a full-sib mating system, under the assumptions of independence and non-interaction among loci, is illustrated. All results are shown to agree with previous work in special cases. All possible full sib, generation n - 1 parent-generation n + m offspring, and generation n uncle-generation n + m nephew covariances for 1 less than n + m less than or equal to 8 are obtained using the given algorithm.