py = linspace(y[0] - rstep, y[-1] + rstep, points)