added a test to check that the center of mass of the model is located in zero when computing Roland_Sylvain integrals