def compute_join_length(px, py, tlow=0.0, thigh=1.0):
    from scipy.integrate import quad
    xp = polyder(px, 1)
    yp = polyder(py, 1)
    xp2 = polymul(xp, xp)
    yp2 = polymul(yp, yp)
    p = polyadd(xp2, yp2)
    integrand = (lambda t: sqrt(polyval(p, t)))
    return quad(integrand, tlow, thigh)[0]
def compute_join_length(px, py, tlow=0.0, thigh=1.0):
    from scipy.integrate import quad
    xp = polyder(px, 1)
    yp = polyder(py, 1)
    xp2 = polymul(xp, xp)
    yp2 = polymul(yp, yp)
    p = polyadd(xp2, yp2)
    integrand = (lambda t: sqrt(polyval(p, t)))
    return quad(integrand, tlow, thigh)[0]