main()
{
int i, k, n, N, kl, kh, nOmega;
float fl, fh, fs, fs1, Q, w, dw;
float xRe[M], xIm[M], Xre[M], Xim[M];
float hRe[ML], hIm[ML], Hdre[ML], Hdim[ML];
float yRe[M], yIm[M], Yre[M], Yim[M];
float Freq[M], HejwRe[M], HejwIm[M], Hejw[M];
float Hre, Him, Alpha, Alpha1, Sita[M];
float f1, f2, t;
clrscr();
N = 22 ;
fl = 500.0 ; fh = 700.0 ; fs = 2200.0 ;
f1 = 50.0 ; f2 = 600.0 ;
1