<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body>
<p>Hi Scilabers</p>
<p>On 17. April I sent an email, but I think it wasn't released to
the mailing list, I only find it online here:</p>
<p><a
href="http://mailinglists.scilab.org/Scilab-users-Covid19-model-td4040626.html">http://mailinglists.scilab.org/Scilab-users-Covid19-model-td4040626.html</a></p>
<p>Here's the email again (without attachment, see above link if
you'd like to view the attachment)</p>
<p style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"><span>A friend (Lars
Risbo) published some MATLAB code in LinkedIn for simulating the
infection with a company lockdown after some time. I figured I'd
try to convert it to Scilab.</span></p>
<span style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"></span>
<p style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"><span>The code uses<span> </span><span
class="bold highlight search-highlight" style="font-weight:
bold; background: rgb(255, 255, 153);">ODE45</span><span> </span>and
I'm not sure that I understand how to convert this MATLAB code
to Scilab. I hope you can explain what I need to do to make the
code work. Some of the original MATLAB code is found in the
comments.</span><br>
<span></span></p>
<span style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"></span>
<pre style="color: rgb(0, 0, 0); font-size: 13.44px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: rgb(255, 255, 255); text-decoration-style: initial; text-decoration-color: initial; font-family: Monospaced;"><span style="color: rgb(92, 92, 92);"></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// covid19risbo.sce</span></span><span>
</span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>//</span></span><span>
</span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// SEIRsim1</span></span><span>
</span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// Susceptible-Exposed-Infectious-Recovered (SEIR)</span></span><span>
</span><span style="color: rgb(176, 24, 19);"><span>function</span></span><span> </span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>dydt</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0); text-decoration: underline;"><span><span class="bold highlight search-highlight" style="font-weight: bold; background: rgb(255, 255, 153);">odefun</span></span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>t</span></span><span style="color: rgb(0, 0, 0);"><span>, </span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>y</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span>
</span><span style="color: rgb(160, 32, 240);"><span>if</span></span><span> </span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>t</span></span><span style="color: rgb(92, 92, 92);"><span><</span></span><span style="color: rgb(0, 0, 0);"><span>days</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>R</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0);"><span>R0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(160, 32, 240);"><span>else</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>R</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(188, 143, 143);"><span>0.6</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// Change reproductive rate after lockdown</span></span><span>
</span><span style="color: rgb(160, 32, 240);"><span>end</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>A</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(74, 85, 219);"><span>[</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(0, 0, 0);"><span>R</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>y</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(0, 0, 0);"><span>Npop</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(255, 170, 0);"><span>...</span></span><span>
</span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(0, 0, 0);"><span>gam</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(0, 0, 0);"><span>R</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>y</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(0, 0, 0);"><span>Npop</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(255, 170, 0);"><span>...</span></span><span>
</span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>gam</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(255, 170, 0);"><span>...</span></span><span>
</span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(0, 0, 0);"><span>Fhosp</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(255, 170, 0);"><span>...</span></span><span>
</span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(0, 0, 0);"><span>Fhosp</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(0, 0, 0);"><span>Thosp</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(255, 170, 0);"><span>...</span></span><span>
</span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(0, 0, 0);"><span>R</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>y</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(0, 0, 0);"><span>Npop</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(74, 85, 219);"><span>]</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>dydt</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0);"><span>A</span></span><span style="color: rgb(92, 92, 92);"><span>*</span></span><span style="color: rgb(131, 67, 16); font-weight: bold;"><span>y</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(160, 32, 240);"><span>end</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>Begin</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0);"><span>datenum</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>2020</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>02</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>20</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// begin date</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>Dlock</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0);"><span>datenum</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>2020</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>03</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>12</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// date of lockdown</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>days</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(0, 0, 0);"><span>Dlock</span></span><span style="color: rgb(92, 92, 92);"><span>-</span></span><span style="color: rgb(0, 0, 0);"><span>Begin</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>R0</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(188, 143, 143);"><span>2.6</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// inital R value</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>gam</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(188, 143, 143);"><span>3</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// gamma, the inverse of average latent time</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>delta</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>1</span></span><span style="color: rgb(92, 92, 92);"><span>/</span></span><span style="color: rgb(188, 143, 143);"><span>5</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// inv time constant which infectious people either recover or enter hospital</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>Fhosp</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(188, 143, 143);"><span>0.16</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// fraction of recovering people going to hospital</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>Thosp</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(188, 143, 143);"><span>14</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// average time of hospitalisation</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>Npop</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(188, 143, 143);"><span>6e6</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// total initial population of sensitive</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>y0</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(74, 85, 219);"><span>[</span></span><span style="color: rgb(0, 0, 0);"><span>Npop</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(188, 143, 143);"><span>50</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(188, 143, 143);"><span>50</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(74, 85, 219);"><span>]</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// [S E I R Hosp ] intial cond.</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>t1</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(255, 170, 0);"><span>:</span></span><span style="color: rgb(188, 143, 143);"><span>60</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>y1</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(50, 185, 185);"><span><span class="bold highlight search-highlight" style="font-weight: bold; background: rgb(255, 255, 153);">ode</span></span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(0, 0, 0);"><span>y0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>0</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>60</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span> </span><span style="color: rgb(0, 0, 0); text-decoration: underline;"><span><span class="bold highlight search-highlight" style="font-weight: bold; background: rgb(255, 255, 153);">odefun</span></span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// [t1,y1] = <span class="bold highlight search-highlight" style="font-weight: bold; background: rgb(255, 255, 153);">ode45</span>(@<span class="bold highlight search-highlight" style="font-weight: bold; background: rgb(255, 255, 153);">odefun</span>,[0 60],y0); // run 1st scernario</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>idx</span></span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span style="color: rgb(74, 85, 219);"><span>[</span></span><span style="color: rgb(188, 143, 143);"><span>3</span></span><span style="color: rgb(255, 170, 0);"><span>:</span></span><span style="color: rgb(188, 143, 143);"><span>6</span></span><span style="color: rgb(74, 85, 219);"><span>]</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>scf</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>a</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(0, 0, 0);"><span>gca</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>plot</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(0, 0, 0);"><span>Begin</span></span><span style="color: rgb(92, 92, 92);"><span>+</span></span><span style="color: rgb(0, 0, 0);"><span>t1</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(0, 0, 0);"><span>y1</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(255, 170, 0);"><span>:</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(0, 0, 0);"><span>idx</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>-</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span> </span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// semilogy(Begin+t1,y1(:,idx),</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>'</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>-</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>'</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>,</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>'</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>LineWidth</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>'</span></span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>,3)</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>a</span></span><span style="color: rgb(92, 92, 92);"><span>.</span></span><span style="color: rgb(170, 170, 170);"><span>log_flags</span></span><span> </span><span style="color: rgb(92, 92, 92);"><span>=</span></span><span> </span><span style="color: rgb(188, 143, 143);"><span>"</span></span><span style="color: rgb(188, 143, 143);"><span>nln</span></span><span style="color: rgb(188, 143, 143);"><span>"</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>xlabel</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Date</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>ylabel</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Number of cases</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span>
</span><span style="color: rgb(50, 185, 185);"><span>xgrid</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// ax=gca;</span></span><span>
</span><span style="color: rgb(100, 174, 100); font-style: italic;"><span>// ax.YLim=[10 max(max(y2(:,idx)))];</span></span><span>
</span><span style="color: rgb(50, 185, 185);"><span>xtitle</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(74, 85, 219);"><span>{</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Danish Corona lock down on 20.03.12, vs 14 days later</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>R goes from 2.6 to 0.6 at lockdown</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(74, 85, 219);"><span>}</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span>
</span><span style="color: rgb(0, 0, 0);"><span>legend</span></span><span style="color: rgb(74, 85, 219);"><span>(</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Infected</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Recovered</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Hospitalised</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(0, 0, 0);"><span>,</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(188, 143, 143);"><span>Total Cases</span></span><span style="color: rgb(188, 143, 143);"><span>'</span></span><span style="color: rgb(74, 85, 219);"><span>)</span></span><span style="color: rgb(0, 0, 0);"><span>;</span></span><span style="color: rgb(0, 0, 0);"></span></pre>
<span style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"></span>
<p style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"><span>P.S. Attached a
simple graph of what the output should look like (file
covid19risbo.jpg < 50 kb), at least partially - because I
deleted an alternative case with a later lockdown date = dashed
lines.</span><br>
<span></span></p>
<span style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"></span>
<p style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"><span>Best regards,</span></p>
<p style="color: rgb(0, 0, 0); font-family: Verdana, Geneva,
Helvetica, Arial, sans-serif; font-size: 13.44px; font-style:
normal; font-variant-ligatures: normal; font-variant-caps: normal;
font-weight: 400; letter-spacing: normal; orphans: 2; text-align:
start; text-indent: 0px; text-transform: none; white-space:
normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width:
0px; background-color: rgb(255, 255, 255); text-decoration-style:
initial; text-decoration-color: initial;"><span>Claus<br>
</span></p>
</body>
</html>