import sys

read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines

def solve(T, p, w, v, gx, gy):
    T = int(T)
    p, w, v, gx, gy = map(float, (p, w, v, gx, gy))
    r = complex(1 + v, w)
    power = [r**i for i in range(T + 1)]
    g = complex(gx, gy) - power[T]
    c = 1 / sum(abs(z) ** 2 for z in power[:T])
    for i in range(T):
        a = c * g * power[T-1-i].conjugate()
        yield (a.real, a.imag)

"""
def check(p, w, v, gx, gy, uvs):
    z = 1
    r = complex(1 + v, w)
    for dx, dy in uvs:
        z *= r
        z += complex(dx, dy)
        p -= dx**2 + dy**2
    return abs(z - complex(gx, gy)) < 1e-6 and p > 0


uvs = ((1, 0), (1, 0), (1, 0), (0, 1), (0, 1))
check(5.5, 0.5, -0.5, -1.125, 2.125, uvs)
"""

N = int(readline())
probs = readlines()

for prob in probs:
    uvs = solve(*prob.split())
    print('\n'.join('{} {}'.format(x, y) for x, y in uvs))