This research article established an optimal global computational algorithm of indices for control systems matrices for single-delay autonomous linear neutral equations. The development of the index matrices exploited the continuity of these matrices for positive time periods, the step-wise backward continuation recursive procedure, change of variables, theory of linear difference equations and careful observations of emerging patterns
