ind1 = np.ceil(np.cos(alpha - m * delphi) * (l2 - l1) + l1)