wd.u.Rd
Williamson-Downs upper bound
wd.u(delt, y.seq, Y1t, ddid)
the value to obtain the DoTT for
possible values of y