function c=fkuepper(dt,pt,n,nt)
% FKUEPPER generate a kupper wavelet
% input (dt,frequency,sides,samples)

pt=1/pt;
c=1:nt;
c=c*0;
nmax=pt/dt;

	for i=1:nmax,
	t=i*dt;
            c(i) = sin (n*pi*t/pt) - (n/(n+2)) * sin ((n+2)*pi*t/pt);

	end



