function T = disco0cuad(N)
%function T = discocuad(N) 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;
    
    'Calculando....'
    
    for tau=0:0.01:2
            for n=1:1:N
            B = B + (8/(j0(n)^3*BesselJ(1,j0(n))))*exp(-1*tau*j0(n)^2);
        end
        T(i) = B;
        i=i+1;
        B=0;
    end
    
    'Modelo generado...'
    
