function [R,T] = discoRcuad(N, h, tfin)
%function T = discocuad(N, h, tfin) 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
j0 = raizBesselJ(N,0.01);
%Calculemos entonces en el punto R=0, por simplicidad
B=0;i=1;k=1;r=1;
'Calculando....'
for R=0:h:1
for tau=0:0.01:tfin
for n=1:1:N
B = B + (8/(j0(n)^3*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
r=r+1;
i=1;
end
R = 0:h:1;
'Modelo generado...'
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
