The magnetic forces between paramagnetic squares cannot be calculated using a classic dipolar model because the magnetic field distribution is not uniform within square particles. Here, we present the calculation of magnetic forces and torques on paramagnetic squares in a uniform 2-D magnetic field using a Laplace's equation solver. With these calculations, we simulate the variations in equilibrium configurations as a function of number of interacting squares. For example, a single square orients with its diagonal directed to the external field while a system of multiple squares will assemble into chain-like structures with their edges directed to the external field. Unlike chains of spherical magnetic particles, that easily stagger themselves to aggregate, chains consisting of magnetic squares are unable to aggregate due to interchain repulsion.