function [f,tau] = cilindroZInf(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.
[Am,betam] = coefAs(N,beta,a);
alpha = 1.47*10^-7;
B=0;i=1;r=1;
R= 1;
for Z=0.1:0.1:a
for tau=0:0.01:T
for n=1:1:N
B = B + Am(n)*besselJ(0,betam(n)*R)*exp(-1*tau*betam(n)^2);
end
f(r,i) = B*erf(Z/sqrt(4*tau));
i=i+1;
B=0;
end
r=r+1;
i=1;
end
tau=0:0.01:T;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
