Computing expansion coefficients in orthogonal bases of ultraspherical polynomials (Q1893555)

The problem is the calculation of the coefficients in the approximation of a function by an expansion in Gegenbauer polynomials. The authors propose to replace the single integrals by double integrals. (This means that the integrand contains only elementary functions which are readily available). It is claimed that there are advantages in this. However the authors choose as basis for comparison two inappropriate methods which will clearly give poor results. The third method with which it is compared is the use of the appropriate Gaussian formulae. The results seem to show that the proposed method has little, if any improvement in accuracy over this. It may be that their method would be useful in situations when the Gaussian formulae were not available. There is no discussion of the calculation of values of the gamma function which occur in the formulae.