function Dk = coefDk(N,C,D)
%function Dk = coefDk(N,C,D) calculara los coeficientes Dk del desarrollo Fourier Bessel del perfil gaussiano
D=0;Dk=0;

for k=0:1:N
    
    H = ((-1)^k)*coeficiente(k,D);
    Cl=0;
    
    for l=0:1:k
        Cl = Cl + coeficiente(l,C);
    end
    
    Dk = Dk + H*Cl;
end
