Employs 3 mathematical models to solve problems of efficient and stable stacking palletization of interlocking boxes of both rectangular and square dimensions.