<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
</head>
<body>
<div class="moz-cite-prefix">Hi Stéphane </div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Thank you for showing me how to work
ODE.</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">P.S. about Corona modeling, the model
by Risbo is nice in that you can implement a lockdown of some
sort. It's a simple one-step lockdown, where you have to tune the
R-value ... in reality each geographic area is different, related
to how much of a lockdown is executed, how dense the population
is, etc.<br>
</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Cheers,</div>
<div class="moz-cite-prefix">Claus</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">On 4/30/2020 4:51 PM, Stéphane Mottelet
wrote:<br>
</div>
<blockquote type="cite"
cite="mid:d0c33eac-46c2-2b54-793d-035bb4e22bcd@utc.fr">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<p>Hi Clauss,</p>
<p>Just modify the following lines :</p>
<pre style="font-family:Monospaced;font-style:normal;font-size:16.0;"><span style="color:rgb(0,0,0);">t1</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(188,143,143);">0</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">0.1</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">60</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">y1</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(50,185,185);">ode</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">y0</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(188,143,143);">0</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">t1</span><span style="color:rgb(0,0,0);">,</span> <span style="color:rgb(0,0,0);">odefun</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">// [t1,y1] = ode45(@odefun,[0 60],y0); // run 1st scernario</span>
<span style="color:rgb(0,0,0);">idx</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(188,143,143);">3</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(188,143,143);">6</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">scf</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(0,0,0);">a</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(174,92,176);text-decoration:underline;">gca</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">Begin</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(0,0,0);">t1</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">y1</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">idx</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);">-</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(100,174,100);font-style:italic;">// semilogy(Begin+t1,y1(:,idx),</span><span style="color:rgb(100,174,100);font-style:italic;">'</span><span style="color:rgb(100,174,100);font-style:italic;">-</span><span style="color:rgb(100,174,100);font-style:italic;">'</span><span style="color:rgb(100,174,100);font-style:italic;">,</span><span style="color:rgb(100,174,100);font-style:italic;">'</span><span style="color:rgb(100,174,100);font-style:italic;">LineWidth</span><span style="color:rgb(100,174,100);font-style:italic;">'</span><span style="color:rgb(100,174,100);font-style:italic;">,3)
</span>And you will be fine. ode() needs a vector of time values and the output is oriented transposed w.r.t Matlab output.
</pre>
<div class="moz-cite-prefix">S.</div>
<div class="moz-cite-prefix"><br>
</div>
<div class="moz-cite-prefix">Le 30/04/2020 à 16:40, Claus Futtrup
a écrit :<br>
</div>
<blockquote type="cite"
cite="mid:8a5ed43a-39a2-342c-de18-afa310eff192@gmail.com">
<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(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></pre>
</blockquote>
<pre class="moz-signature" cols="72">--
Stéphane Mottelet
Ingénieur de recherche
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
<a class="moz-txt-link-freetext" href="http://www.utc.fr/~mottelet" moz-do-not-send="true">http://www.utc.fr/~mottelet</a>
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
users mailing list
<a class="moz-txt-link-abbreviated" href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>
<a class="moz-txt-link-freetext" href="http://lists.scilab.org/mailman/listinfo/users">http://lists.scilab.org/mailman/listinfo/users</a>
</pre>
</blockquote>
<p><br>
</p>
</body>
</html>