ind2 = np.floor(np.cos(alpha - m * delphi) * (l2 - l1) + l1 + diameter)