vl, vr = map(float, input().split())
d = float(input())
w = float(input())

# d0 == (w + vr) * t1
# t1 = d0 / (w + vr)
# d1 = d0 - (vl + vr) * t1 = d0 * (w - vl) / (w + vr)
# d1 == (w + vl) * t2
# t2 = d1 / (w + vl) = d0 * (w - vl) / (w + vr) / (w + vl) = t1 * (w - vl) / (w + vl)
# d2 = d1 - (vl + vr) * t2 = d1 * (w - vr) / (w + vl)
# d2 == (w * vr) * t3
# t3 = d2 / (w + vr) = d0 * (w - vl) / (w + vr)**2 / (w + vl) = t1 * (w - vl) / (w + vr) / (w + vl)
# d3 = d2 - (vl + vr) * t3 = d2 * (w - vl) / (w + vr)
# ans = w * (t1 + t2 + t3 + ...)
#     = w * d0 * (1 + (w - vl) * (1 + 1/(w + vl) + 1/(w+vr)/(w+vl) + ...)) / (w + vr)

i = 0
ans = 0.0
while d > 0.000001 and i < 10000000:
    t1 = d / (w + vr)
    ans += w * t1
    d = d - (vl + vr) * t1
    t2 = d / (w + vl)
    ans += w * t2
    d = d - (vl + vr) * t2
    i += 1

print(ans)