sos100=sumofsquares(100, 10, zerosum=TRUE, maxt=Inf),
sos200=sumofsquares(200, 10, zerosum=TRUE, maxt=Inf),
sos400=sumofsquares(400, 10, zerosum=TRUE, maxt=Inf),
sos800=sumofsquares(800, 10, zerosum=TRUE, maxt=Inf), and
sos1000=sumofsquares(100, 10, zerosum=TRUE, maxt=Inf)