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

[As,alphas] = coefAs(N,beta,a);
k=1;

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