\vec g(\vec r_1+\vec vt_1,t_1)=G\int_V\frac{\rho(\vec r\,'_1+\vec vt_1,t_1)((\vec r\,'_1+\vec vt_1)-(\vec r_1+\vec vt_1))}{|(\vec r\,'_1+\vec vt_1)-(\vec r+\vec vt_1)|^3}\,d\vec r\,'_1