def interval_intersection_width(a, b, c, d):
return max(<NUM_LIT:0>, min(b, d) - max(a, c))<EOL>
returns the width of the intersection of intervals [a,b] and [c,d]
    (thinking of these as intervals on the real number line)