iarspider, а Вы ненароком подкинте теор. отделу идейку, что общее решение не нужно, пусть поищут частное, форму замкнутых в пределах одного оборота орбит в виде X=A*sin(a)+C, Y=B*cos(a)+D*cos(2*a)+E, может и выйдет что путное
(Х - ось проходит через неподвижные тела, 0 в центре тяжести их)