function [da, x, dax] = kb_deap(f,l,n)

x = [-f*n/2:f*n/2-1]/(n);
b = (f-1/2)*pi*l;

sqa = sqrt(pi*pi*l*l*x.*x-b*b);
dax = sin(sqa)./(sqa);
dax = dax/dax(f*n/2);
da = dax'*dax;

