v[i] = calc_volume(r[i + 1], d) - calc_volume(r[i], d)