%
%          hiv_fun
%

  function Y = hiv_fun(time,params)

  global coef Y0 wt

  parameters = params;

  ode_options = odeset('RelTol',1e-6);
  [t,Y] = ode15s(@hiv_rhs_dram,time,Y0,ode_options,parameters);