%
%   Fuel_Temp_Exact_Int.M  Function to evaluate integral of temperature  
%                               profile for variable kf case
%                              (called from Fuel_Temp_Exact)
%
%   File prepared by J. R. White, UMass-Lowell (Oct. 2004)
%
      function int = eqnfile(r)
      global Ts qqq Rf a1 a2
      c = 1+a2;
      T = (Ts^c + (c/a1)*(qqq/4)*(Rf*Rf - r.*r)).^(1/c);
      int = 2*pi*r.*T;
%
%   end of function
