function [t,T,R,tasa] = tasaRgauss(N,beta,a,T)
%function T = tasaRgauss(N,beta,a,T) nos va a permitir representar informacion del campo de temperaturas resultadnte de resolver
%la ecuacion de calor para el primer modelo de aproximacion al resultado. Para simplificar las variables se van a adimensionalizar
% Tint: Temperatura del interior del disco
% Text: Temperatura del exterior del disco
% N: Numero de sumandos para el desarrollo en serie
%Las nuevas variables de trabajo van a variar de la siguiente forma.
% R: variara entre 0 y 1
% theta: Variara entre 0 y 2*pi
%Obtengamos en primer lugar las N raices de la funcion Jo[x] que nos interesan
[As,alphas] = coefAs(N,beta,a);
%Calculemos entonces en el punto R=0, por simplicidad
B=0;i=1;k=1;r=1;
t = [0:0.01:T];
'Calculando....'
for R=0:0.01:1
for tau=0:0.01:T
for n=1:1:N
B = B + As(n)*besselJ(0,alphas(n)*R)*exp(-1*tau*alphas(n)^2);
end
T(r,i) = B;
i=i+1;
B=0;
end
tau1 = buscaTemperatura(0.8*T(r,1),T(r,:),t);
tau2 = buscaTemperatura(0.2*T(r,1),T(r,:),t);
tasa(k) = 0.6*T(r,1)/(tau2-tau1);
k=k+1;
r=r+1;
i=1;
end
R=0:0.01:1;
'Modelo generado...'
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
