def compute_join_curvature(px, py):
    from scipy.integrate import quad
    xp = polyder(px, 1)
    xpp = polyder(px, 2)
    yp = polyder(py, 1)
    ypp = polyder(py, 2)
    pn = polyadd(polymul(xp, ypp), polymul(yp, xpp))
    pd = polyadd(polymul(xp, xp), polymul(yp, yp))
    integrand = (lambda t: fabs((polyval(pn, t) / (polyval(pd, t) ** 1.5))))
    return quad(integrand, 0, 1)[0]
def compute_join_curvature(px, py):
    from scipy.integrate import quad
    xp = polyder(px, 1)
    xpp = polyder(px, 2)
    yp = polyder(py, 1)
    ypp = polyder(py, 2)
    pn = polyadd(polymul(xp, ypp), polymul(yp, xpp))
    pd = polyadd(polymul(xp, xp), polymul(yp, yp))
    integrand = (lambda t: fabs((polyval(pn, t) / (polyval(pd, t) ** 1.5))))
    return quad(integrand, 0, 1)[0]