D. Slapian (1977): "Program to compute discrete prolate spheroidal wave functions, sequences, and eigenvalues." Technical Report TM 77-1218-8, Bell Telephone Laboratories.