def segment_area(seg):
r0 = seg[<NUM_LIT:0>][COLS.R]<EOL>r1 = seg[<NUM_LIT:1>][COLS.R]<EOL>h2 = point_dist2(seg[<NUM_LIT:0>], seg[<NUM_LIT:1>])<EOL>return math.pi * (r0 + r1) * math.sqrt((r0 - r1) ** <NUM_LIT:2> + h2)<EOL>
Compute the surface area of a segment.

    Approximated as a conical frustum. Does not include the surface area
    of the bounding circles.