for (j=0;j<M;j++) {
jj=M*(i%2)+j;
data[j]=10.0*exp(-1.*SQR(jj-M)/SQR(N/8))*
(cos(2.0*pi*32*jj/N)+cos(2.0*pi*64*jj/N));
}
sp.adddataseg(data,hann);
}
spec=sp.spectrum();
freq=sp.frequencies();
// for (i=0;i<M+1;i++)
// System.out.printf(setw(20) << freq[i] << " " << setw(20) << spec[i]);