vl, vr = map(float, raw_input().split())
d1 = float(raw_input())
w = float(raw_input())

EPS = 0.000000001
a = 0
for i in xrange(10**7):
    t1 = d1/(vr+w)
    a1 = w * t1
    if a1 < EPS:
        break
    a += a1
    d2 = d1-(vl+vr) * t1
    t2 = d2 / (vl+w)
    a2 = w * t2
    if a2 < EPS:
        break
    a += a2
    d1 = d2-(vl+vr)*t2

print a