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.01);
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
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
