I think this can be expressed as:

`lorenz_points = (arange(1, int(n) + 1) / n, sdata.cumsum() / Sn)`