A novel analytical design method for highly selective digital optimal equiripple comb FIR filters is presented. The equiripple comb FIR filters are optimal in the Chebyshev sense. The number of notch bands, the width of the notch bands and the attenuation in the pass bands can be independently specified. The degree formula and the differential equation for the generating polynomial of the filter is presented. Based on the differential equation, a fast simple algebraic recursive procedure for the evaluation of the impulse response of the filter is described. Its arithmetic robustness outperforms by far the known analytical design method. Highly selective equiripple comb FIR filters with thousands of coefficients can be designed. One example demonstrates the efficiency of the filter design.