Bailey, Cheng, and Kipnis [3] developed a method for constructing trend-free run orders of factorial experiments called the generalized fold-over method (GFM). In this paper, we use the GFM of constructing run orders of factorial experiments to give a systematic method of constructing magic squares of higher order.