工程計算MATLAB code to calculate the reorthogonalized sine tapers
input: N = the length of the time series data to be tapered
p = the number of tapers requested
I = the gap structure a vector of length N
I(t) = 1 if there is data at time t, t=1, ..., N
I(t) = 0 if there is a gap at time t
output: X = N-by-p vector of the reorthogonalized sine taper
computes the eigenvalues of a symmetric tridiagonal
* matrix T. The user may ask for all eigenvalues, all eigenvalues
in the half-open interval (VL, VU], or the IL-th through IU-th eigenvalues.