function [R, tasa] = tasaRconst(N,h, tfin)
%function [R,tasa] = tasaRconst(N, h, tfin) nos va a permitir representar informacion del campo de temperaturas resultadnte de resolver
%la ecuacion de calor para el modelo uniforme 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
% h: Precisión en la variable temporal
% tfin: tiempo adminesional en el que se estudiará la evolución de la temperatura
%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
j0 = raizBesselJ(N,0.01);
%Calculemos entonces en el punto R=1, por simplicidad
B=0;i=1;k=1;r=1;
t = [0:0.01:tfin];
'Calculando....'
for R=0:h:1
for tau=0:0.01:tfin
for n=1:1:N
B = B + (2/(j0(n)*BesselJ(1,j0(n))))*BesselJ(0,j0(n)*R)*exp(-1*tau*j0(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:h:1;
'Modelo generado...'
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
