A set n distinct numbers taken from the interval [1, n^2] form a magic series if their sum is the nth magic constant M_n = 1/2 n(n^2 + 1) (Kraitchik 1942, p. 143). If the sum of the kth powers of these numbers is the magic constant of degree k for all k element [1, p], then they are said to form a pth order multimagic series.