The principal results of this
paper are the following: Every integral 2 × 2 matrix is the sum of at most 3 integral
squares, and this is best possible. Every integral n×n matrix with n > 2 is the sum
of at most k integral squares, where k = 7 if n is even, and k = 9 if n is odd. Every
n × n matrix over GF(2) is the sum of at most 2 matrix squares, and this is best
possible.