def METHOD_NAME(skeleton, ds):
    lengths = np.linalg.norm(skeleton[1:] - skeleton[:-1], axis=-1)
    dists = np.concatenate([np.array([0]), np.cumsum(lengths)])
    ts = inverse_interpolate(dists, ds * dists[-1]) 
    return ts / (len(skeleton) - 1)