sizeofsys -> double   (sizeofsys = pow ( numberofpart / rho , 1.0 / dimofsys ) )