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)
# d2 = d1 - (vl + vr) * t2 = d1 * (w - vr) / (w + vl)
# d2 == (w * vr) * t3
# t3 = d2 / (w + vr)
# d3 = d2 - (vl + vr) * t3 = d2 * (w - vl) / (w + vr)
# ans = w * (t1 + t2 + t3 + ...)
#     = w * ((d0 + d2 + ...) / (w + vr) + (d1 + d3 + ...) / (w + vl))

ans = 0.0
while d > 0.000001:
    t1 = d / (w + vr)
    ans += w * t1
    d = d - (vl + vr) * t1
    t2 = d / (w + vl)
    ans += w * t2
    d = d - (vl + vr) * t2

print(ans)