This file contains the information of and the instructions for running the programs
used for the following paper submitted to RED:
Title of the paper: "Knowledge, Germs and Output"
RED manuscript number: RED-D-20-00358
Author: Shouyong Shi
RED editor in charge: Loukas Karabarbounis

General information of the programs:
    Software and operating system used: matlab 2021a; Windows 10 (64 bits)
    The order in which the programs need to run: see the details below
    Expected computation time: 0.5 minute if "tdone=1"; 10 minutes for each combination if "tdone=0"
    Random numbers: no random numbers were used.

1. The main program is 'kgo2_master'.
All programs without the prefix 'kgo2' in the filenames are for Chebyshev projection,
which are called by some of the programs with the prefix 'kgo2'.
A description of Chebyshev projection appears in 'chebyshev-description.pdf'.

2. A trial run of the program:
Execute 'kgo2_master' to produce Figure 1 in the paper. This requires that all the data
files have already been generated.

3. The steps for generating the data files: 
% Set lock=0 to run program without lockdown:
%   (1a) Set tdone=dosteady=0. Run estart =2,1,3 for three initial states.
%       ALWAYS RUN FIRST WITH estart=2 TO GET THE BASELINE VALUES OF Y AND
%       rU2 AT t=0.
%   (1b) Set tdone=dosteady=1 to run for steady state with estart =2,1.
%   (1c) Set tdone=dosteady=0. Set the initial (n0S0,n1S0) in Step 3 to
%   the steady state (N0,N1) in (1b) above with estart=1. Run estart=4.
%   (1d) To solve the planner's problem, set tdone=1, dosteady=0, estart =2
%   and doplan=1 (in Step 5).
% Lockdown: Set lock=1 in Step 2, estart =2, tdone=dosteady=0.
%   (DO NOT RUN LOCKDOWN FOR THE STEADY STATE, since this requires recoding
%   to get rid of the variation in Lt in subsequent runs of the program.)
%   (2a) Set tlock=11 to run with an early lockdown.
%   (2b) Set tlock=19>=16 (<thigh) to run with a delayed lockdown.
% Sensitivity analysis: Set estart =2, tdone=dosteady=lock=0. Go to
%   '_fspara' to change the parameters before running the program:
%   (3a) Set l_lower to 0.05>=0.01: reset Tinfa=6.7/365 to meet target on R0.
%   (3b) Set umult to 1.50>1: reset Tinfa=5.8/365 to meet target on R0.
%   (3c) Set apmult to 1.40>1: reset Tinfa=6.7/365 to meet target on R0.
%   (3e) Set improb=1 to run for rho=0, or improb=0.5 to run for a low mu:
%       do not reset Tinfa.
%   (3f) Set phi=0 for a vaccine: reset Tinfa=6.2/365 to meet target on R0.

4. The steps for producing all of the graphs:
The graphs are produced by (1)-(8) in Step 6 in 'kgo2_master'.
% DO THIS AFTER FINISHING RUNNING THE PROGRAM FOR ALL CONTROL VALUES.
% THEN SET tdone=1. The default value of estart is 2.
% (i) For each block (1)-(8) below, eliminate '%' in the block (except
% for 'clf'). PUT '%' IN ALL OTHER BLOCKS. RUN 'kgo2_master'.
% (ii) After producing a graph for a block, put '%' back on the block
%   and reset the parameter and controlling variables to the baseline.
% (iii) If a blcok is intended to producing overlapping graphs with the
% controls (c1,c2,...), set the corresponding controlling parameter to
% (c1,c2,...), run 'kgo2_master', and then re-run for c2, etc.