b=2.0;
sbeps=5.e-6;
DftInt dftint = new DftInt();
for (i=0;i<N;i++) {
w=0.1*i;
dftint.dftint(this,a,b,w,cosint,sinint);
// System.out.printf(cosint << " " << sinint);
prefa=exp(-a)/(1.0+SQR(w));
prefb=exp(-b)/(1.0+SQR(w));
cwa=cos(w*a);
swa=sin(w*a);