function zombie_ode_test02 ( )

%*****************************************************************************80
%
%% zombie_ode_test02 runs a specific zombie test case.
%
%  Licensing:
%
%    This code is distributed under the MIT license.
%
%  Modified:
%
%    22 April 2020
%
%  Author:
%
%    John Burkardt
%
  fprintf ( 1, '\n' );
  fprintf ( 1, 'zombie_ode_test02:\n' );
  fprintf ( 1, '  Demonstrate the eradication model.\n' );

% alpha = 0.005;
  alpha = 0.008;
  beta = 0.0095;
% zeta = 0.0001;
  zeta = 0.02;
  delta = 0.0001;
  k = 0.08;
  t = 20.0;
  dt = 0.01;

  fprintf ( 1, '\n' );
  fprintf ( 1, '  ALPHA = %g, death rate of zombies in human/zombie encounter.\n', alpha );
  fprintf ( 1, '  BETA =  %g, infection rate for human/zombie encounter.\n', beta );
  fprintf ( 1, '  ZETA =  %g, dead-human-zombie-resurrection rate.\n', zeta );
  fprintf ( 1, '  DELTA = %g, Non-zombie-related human death rate.\n', delta );
  fprintf ( 1, '  K     = %g, zombie kill rate.\n', k );
  fprintf ( 1, '\n' );
  fprintf ( 1, '  T =     %g, length of time period to be studied.\n', t );
  fprintf ( 1, '  DT =    %g, stepsize for Euler method.\n', dt );

  erad ( alpha, beta, zeta, delta, k, t, dt );

  return
end
