function [f,x,g] = FBcuadComplex(N)
%function [f,x,g] = FBgaussian(N,beta,a) obtiene la aproximacion Fourier-Bessel del perfil gaussiano

[As,alphas] = cuadComplex(N);
k=1;

for x=0:0.01:1
    B=0;
    for i=1:1:N
        B = B + As(i)*besselJ(0,alphas(i)*x);
    end
    
    f(k) = B;
    g(k) = (1-x)^2;
    k=k+1;
end
x=0:0.01:1;    