Middle-square method