% unitsquare.m sets up the unit square US = [[0 0]; [0 1]; [1 1]; [1 0]; [0 0]]';