function [f,tau] = discoRgauss(N,beta,a,T)
%function [f,tau] = disco0gauss(N,beta,a,T) obtiene la aproximacion Fourier-Bessel del perfil gaussiano
%Para que la raiz la obtengamos en 1, que ser´a los que nos interesara gamma se calcula a partir de beta.

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

B=0;i=1;r=1;

for R=0:0.1:1.0

    for tau=0:0.01:T
        for n=1:1:N
            B = B + As(n)*besselJ(0,alphas(n)*R)*exp(-1*tau*alphas(n)^2);
        end
    
        f(r,i) = B;
        i=i+1;
        B=0;
    end
    r=r+1;
    i=1;
end
tau=0:0.01:T;    