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

      function zp = odefile(r,z)
      global qqq a1 a2
      zp = zeros(length(z),1);
      kf = a1*z(1)^a2;   kfp = a1*a2*z(1)^(a2-1);
      zp(1) = z(2);
      if r > 0
        zp(2) = -qqq/kf - (1/r + (kfp/kf)*z(2))*z(2);
      else
        zp(2) = -qqq/(2*kf);
      end
%
%   end of function
