void Gra(int O, int T){
	Axes P = {Sun[O].xM, Sun[O].yM, 0, 0}; // 惯性/力 
	Axes G, GT; // 合力 和 对象施的力 
	double GT_; //对象力的大小 
	GT_ = g * SunWeight * SunWeight / (Dis(Sun[O].x, Sun[O].y, Sun[T].x, Sun[T].y) * Dis(Sun[O].x, Sun[O].y, Sun[T].x, Sun[T].y));
	GT.x = GT_ * sqrt((Sun[O].x - Sun[T].x) * (Sun[O].x - Sun[T].x)) / Dis(Sun[O].x, Sun[O].y, Sun[T].x, Sun[T].y) - Sun[O].x;
	GT.y = GT_ * sqrt((Sun[O].y - Sun[T].y) * (Sun[O].y - Sun[T].y)) / Dis(Sun[O].x, Sun[O].y, Sun[T].x, Sun[T].y) - Sun[O].y;
	G.x = P.x + GT.x;
	G.y = P.y + GT.y;
	Sun[O].xM = G.x;
	Sun[O].yM = G.y;
	std::cout<<Sun[O].xM<<" "<<Sun[O].yM<<std::endl;
}

我不知道我写的对不对QAQ