This function computes and returns properly-normalized Slepian tapers.
These tapers are normalized so that

(16.20.331) |

`num_points:`Input. The number of points in the taper.`nwin:`Input. The number of tapers computed.`lam:`Output. Upon return,`lam[0..nwin-1]`contains the eigenvalues of the tapers. Note that .`nwdt:`Input. The (total sample time) (frequency resolution bandwidth) product.`tapers:`Output. Upon return:`tapers[0..num_points-1]`contains the first taper,`tapers[num_points..2*num_points-1]`contains the second taper, and so on.`tapsum:`Output. On return`tapsum[0]`contains the sum of the`num_points`values of the first taper,`tapsum[1]`contains the sum of the values of the second taper, and so on. Note that because the odd-index Slepian taper functions are odd,`tapsum[1,3,5,...]`would vanish if it were not for round-off and other numerical error.

This function will print a warning message if the condition is not satisfied (see Section ).

- Author: Adapted from the original code (Lees and Park) by Bruce Allen (ballen@dirac.phys.uwm.edu) and Adrian Ottewill (ottewill@relativity.ucd.ie).
- Comments: There are a number of techniques for calculating the Slepian tapers. We have not extensively tested these routines, but they appear to work well. They make use of the standard EISPACK routines, translated from FORTRAN into C using f2c.