function [f,tau,tasa] = tasa0gauss(N,beta,a,T)
%function [f,tau] = tasa0gauss(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.
%gamma = 1/(exp(1/beta)-1)
gamma=0;
[As,alphas] = coefAs(N,beta,a);
%R=0;
k=1;
for tau=0:0.01:T
B=0;
for n=1:1:N
%B = B + As(n)*besselJ(0,alphas(n)*R)*exp(-1*tau*alphas(n)^2);
B = B + As(n)*exp(-1*tau*alphas(n)^2);
end
f(k) = B;
k=k+1;
end
tau=0:0.01:T;
tau1 = buscaTemperatura(0.8*f(1),f,tau);
tau2 = buscaTemperatura(0.2*f(1),f,tau);
tasa = 0.6*f(1)/(tau2-tau1);
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
