middle-square method simulation Cplus