%
%   Fuel_Temp_SH1_Eqns.M   Function file for Fuel_Temp.m
%
%   Given differential eqn.:   T'' + (1/r)T' + qqq/kf =  0
%
%   File prepared by J. R. White, UMass-Lowell (Sept. 2003)
%

      function zp = odefile(r,z)
      global qqq kf
      zp = zeros(length(z),1);
      zp(1) = z(2);
      if r > 0
        zp(2) = -qqq/kf - z(2)/r;
      else
        zp(2) = -qqq/(2*kf);
      end

%
%   end of function
%
