function [As,alphas] = cuadComplex(N)
%function [As,alphas] = cuadComplex(N) nos devolvera los coeficientes del desarrollo Fourier-Bessel del perfil cuadratico aproximado
%(1-R)^2, que nos permitir´a de alguna forma emular el comportamiento de una gaussiana.
alphas = raizBesselJ(N,0.1);
i=1;
D=0;
for i=1:1:N
A = (alphas(i)/2);
B = 1/(2*(besselJ(1,alphas(i)))^2);
D=0;
for k=0:1:N
D = D + (-1)^k*(coeficiente(k,A))^2/((k+1)*(k+2)*(k+1.5));
end
As(i)=B*D;
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
