The ellipsoidal wave equation, which is the most general equation derived by separation of the Helmholtz equation in confocal coordinates, presents unusual computational difficulties, and its solutions, despite their importance for physical applications, have not hitherto been effectively computed. This paper describes a successful technique, which involves the solution of a four-term recursion and the simultaneous handling of two eigenparameters.