function [As,alphas] = cuadComplex(N)
%function [As,alphas] = cuadComplex(N) nos devolvera los coeficientes del desarrollo Fourier-Bessel del perfil cuadratico aproximado
%(1-R)^2, que nos permitir´a de alguna forma emular el comportamiento de una gaussiana.

    alphas = raizBesselJ(N,0.1);
    i=1;
    D=0;
    
for i=1:1:N
    A = (alphas(i)/2);
    B = 1/(2*(besselJ(1,alphas(i)))^2);
    D=0;
    
    for k=0:1:N
        D = D + (-1)^k*(coeficiente(k,A))^2/((k+1)*(k+2)*(k+1.5));
    end
    As(i)=B*D;
end

