Numerical problems in the computation of ellipsoidal harmonics