%
%            function SIR_fun3
%

  function  Y = SIR_fun3(time,params)

  global N Y0

  ode_options = odeset('RelTol',1e-2);
  [t,Y] = ode45(@SIR_rhs3,time,Y0,ode_options,params);