function [As,alphas] = coefAs(N,beta,a)
%function As = coefAs(N,beta,a) con esta funcion trataremos de calcular los coeficientes del desarrollo de Fourier-Bessel del
%perfil cuadratico con el que estamos trabajando para un determinado numero de puntos, para una beta determinada y en un
%intervalo [0,a] que nosotros elijamos arbitrariamente grande.

    alphas = raizBesselJ(N,0.1);
    alphas = alphas/a;
    i=1;
    
    C = (a^2)/beta;
    
for i=1:1:N
    A = (alphas(i)/2)^2;
    D = A*beta;
    B = beta/(a*besselJ(1,a*alphas(i)))^2;
    As(i) = B*(exp(-D) - exp(-C)*coefDk(N,C,D));
end
