#include #include #include #include #define EPSILON (1e-10) #define iszero(x) (fabs(x) < EPSILON) using namespace std; ios& fixed(ios &io) { io.setf(ios::fixed, ios::floatfield); return io; } int main(void) { ifstream cin("solar.in"); double a1, b1, t1, a2, b2, t; int i = 0; while(cin >> a1 >> b1 >> t1 >> a2 >> b2 >> t) { if(a1 == 0.0) break; double t2 = sqrt((t1 * t1 * a2 * a2 * a2) / (a1 * a1 * a1)); t = fmod(t, t2); double f0 = 2.0 * M_PI * a2 * t / t2; double u, u0 = 0.0, u1 = 2 * M_PI; while(true) { u = (u0 + u1) / 2.0; double f = a2 * u - sqrt(a2 * a2 - b2 * b2) * sin(u); if(iszero(f - f0)) break; if(f < f0) u0 = u; else u1 = u; } cout << "Solar System " << (++i) << ": " << fixed << setprecision(3) << a2 * cos(u) << " " << b2 * sin(u) << endl; } return 0; }