From 4f2a7075cb0ca8667690558760278591e8f2e987 Mon Sep 17 00:00:00 2001
From: Paul CACHEUX <paulcacheux@gmail.com>
Date: Thu, 14 Jun 2018 22:49:42 +0200
Subject: [PATCH] Add first result figure

---
 Figures/base_infection_first_result.png | Bin 0 -> 43987 bytes
 SRI_matlab/main.m                       |   2 +-
 SRI_matlab/step.m                       |   6 +++---
 3 files changed, 4 insertions(+), 4 deletions(-)
 create mode 100644 Figures/base_infection_first_result.png

diff --git a/Figures/base_infection_first_result.png b/Figures/base_infection_first_result.png
new file mode 100644
index 0000000000000000000000000000000000000000..93c1ff3730b7e23b97fedb72e682f85a1554099d
GIT binary patch
literal 43987
zcmeAS@N?(olHy`uVBq!ia0y~y;JU-Wz~;fh#K6EXc|)Qt0|NtRfk$L90|U1(2s1Lw
znj^u$z`$PO>Fdh=h>cIgK<VYVcZ>`S3@RnA5hcO-X(i=}MX3zFjv+pdPC5!98L0}s
zi6t4~`9;~qItre7$$B;scR@xdfz&$}r6!hS=I1GdWag$a7@FvMmFDRh85mk97#Ufa
zSXddG$URv$mw`cn!PCVtq+-sSyY&TEzSpY%`<~S*;O3FV#3Cc@$=n_!!;+`BcIE<>
zM?2cD9O2AdS+m5pU;C)%GTD#H#m3KkswXr*+7`@or*zJ)AnjX^8Z~8$&+%W3li_w|
zKbJ0@#KigU(S}8RzTGzMI(-QT-_`y9VKIBHfRj+vT9*Ff^7Fp$vwwg7_vc?VTTThQ
zEx!Lg!fGxb0~9E)Rk+UzWixD7TgU=sE2y^!L75G61iYY3h9V~oP+)<$CpbqHkA}u*
z+8Dk$C1K9bcjfzc7d<^y|Ns5}zo({Z$Jcy3S~Hy)l>HePV#0Oi`aR#x)^PCGVf#Oa
z=USJy3E2Jn*uVe7A#TN%v$M_rAItyaw(m8lFlJz2_@y|7p+|aU$Sd*skK)c8>i+Y3
zjz11x7xVK&yM4{`x#f2M|9t+b$IH-Q`1aG&)8hJZG6@H+t`1-BKR+(OVA5*Nvhwoy
zokglG0YO1QF)=o5>N16at-@+P8k(A&#~&Y_r}*Tz+;I=d#mOo8-?mkiFfcIiOnNKb
zQ26Yz{QrtC-_LRitKF%5J~w~g&$NBoObiKYZg0uFyoK@X9Lvw^|Gir8_HB`<sOZHf
zv(0j)UhQ1}|BQ$4jO@F2?>+^C%gcPV_4WU+|NnJAU;aG9qW<J1Jr@n%%$6;GAH=}G
z@X1w!QNo;yi))$B%uk27_5XZmx93sZ^W#x>{=T2j-0eP1thvt0a3I3G_SYBbKFj)l
ze@?`eTw4>_yuWH&ijn2>^mFrUXTQ|f*YB57RaI3Eed3|1umAqtyK{SF%3713K>Ssa
zf92y3Jyr$=g?O&du?z<UY<|Dl>^)uY@15s%(_cz{)ZfC$P*5#ZQQv-1r2FWvxBK7k
z`|aiJeZTVg+*WS!z9LWN)Bo1hxVj!)ZufQNvSrI2`kK1EE9{wKIIVYgxh4YxL!}BR
zhEApwy}Yzkv1MoR^LzFG|JMKBUVr=Y<<A@YYu_xL9(SrdlYyZ}T3jz?hDG6{t$&5p
zd?b$LOxKUE`+m3FnIm|akEVvi%%dEaUzRMZeR^uDPxbeAca?+O^7G&O*Z;cg$R~G3
zb(8#EZv8zIem?taz*EpaF+}U?<;%u}kB&Tg^hoin%-vEOx%noUL0jbSoSkC0CtafM
zo5bUfKMrtn^891BI>FhlIMwUxw{QD?zuPUW?l;FE(P@6Wnc4lN3=9lCzW?fFC$NK(
z%+8`zqnUfYUD%j>T(17lNBz2w-PzaHSXO*U__okE{alNHT-}d{`E{RXuMS^d_j2j<
zXXobLzI=J|-Utl_2A-T#Q#6nM?L7W?>eQ)EpFZ6f@n&JW+^L@TZ{PMFO$w~f&dy#|
zf4on2V?@rSC7$5)@$_TGo^9LAs=1r9uU@(0^7+k;jf>CQ1lDU!UA1P7j?SN--`Ss?
znK^l9UT*H&uYU47dR_$Rh*_7tiKwmpd!$p?xcU0}`1-@@W(QPv?b?;JF~Zf=HLbwa
z)m2CA_4oNNnk)BQzn8r1$GUavu3fuUVG|b@XP11e=b+7*6~$*~nQn~G+4Wj8c$tWp
z*tF@><t>X+8lvSE^enVHe5UjF_xI~puI!Y1JVi4&Y2%H{{`R&@mfK$8WdGms%cAmA
zN@{AV$O8?1eNc!eCMKq)rlxuJ7=n{og3jh3rU~rve;$e7ulam-{+}o6wfkJY&OBf9
zZ0q&7(|i7$E8i#V^y49ay+p0?v+Dc5<2Vw3etJ5SvHn|n{kH?m{3-t}c^DYjrv2r8
zZ|CLhZEbDs%u)UAP35nb%NN(T%GuRev<TET{JMF*&P*yYGE%blOxcSI3qO7Oba}b|
z`FXa#Z=IZc^2x`?$Jdv?%=njLRQYy`V$0{}=i3GT{ha@Q=lSP5nl+VY`lOxP_5EJ;
zZ=QR5DvLuyLz9yqf8~<?A#~kZ`}gY4lFgNS{=eP-cl(>_GkSHmwq(xyxBLCR-);QH
zLBYYrZ%@2?_s+FjZ0q*z`qA6g+?V`i^D5E)@EhURuU`F<`FluwU%{s*Cw*s|rB>zr
zV3z&U(p<U!Yxe(tf3KgDaeDFQ%^c(ObJ6*GO>=K;+4v2VGZT6?2Qe9V1+Ba?RXhA%
z&F8b{Yn~Z{N*s=b{xcq||Nm?KVwXeV`#w#r|FyjS)}G4ZEB}A*|NlF1@wRPdEg$Ot
zo&SGw&Zn>I>up;KzP`E|xRaTIA&<4SwRL0ml`B_T1pfW#|8HUZ#@+r`;FTYX?SCy^
z5f~O`R{H9Sz#(COn~8rFU%$G#`uX|!^D8aW&&{dZ-yqk0G%q4rM=U!71WtUpa`mdI
ziHT0k4uh6}sHj=C6%$fps_)v#Th@PF9ltlu)_l*AGc%2)rKC26d#zZ#Iy)z4&6+iP
z%C;S}xo7viVsHDm^Z);xzrNZsy1cynS@X}M-Qx1LRbOVN&pWC2wMzZoh8{<|!`c2h
z85t3KDhi)HeRg*C`MK8OTiHLu0#Ez*wQV^!Kb<i?pR!-^_N`mTZtWIfU})HQ{(5&~
z;HKo`ed_aTiXJd3wydxFy1L))R|dyJ^Sa08r%s=)|6X2S{{G(HK!?P4cXop7m-$wu
zSvNPO`rG|{a$0}CjMI*pw%Lpf3Xdnx|C6$8Z|Uo6)22-e(Yku?p4``(x7qRBic5l4
zmcG80+7`RF%2ckuU7+mEjg47XSFI@i{_d`J_&OOkMpjnV-R1B36yMz3tp4if>HU9B
zOB$!eSQp($onv2b*A<qsZ~wi!cjx-3Exz1d)mi)d+uM74f8W_zY@B{hX2%1A9mlOV
zY}&MGdVJkaX89inmU>UO`_RaLr=TW0Jbc;TtLy7ZA0BGu7T2>We00S2`<>vhu&}te
zy2FycOn;=`a8FN9|56|K?d|R9k&Jm?UtQf@{$B10ze)A;?nZT)duvnkH*T-}{q5Pa
zX99<s`R!&@I^QciE}Q%+mWzSmz>b8aED2NoT+9D^ZEpEJ%kp=3G@>uK8_j(4@+GID
z&7Tj4XP!F~ue*4v=%bH6y!H2%JfB<s?(XjGD_6dJIz9g133t0s&GLT|LQb(WJP0fM
z_U5L5Q($1=Gw<hX*REB3&HeXb`u`x?;%7d~mM>nc{8=;e%g5X&H`C|uotd;rJn9Fi
zdGOUbts>RL%*<{6*EWG`>tbiWeE&WAYx91WvWeHVEqLC2<a_gW`~7*g)n%TZp3n5y
z%ql-U@u~j#>FLIZH7i%X{PXj(t*x!NUW|u|Q1i?s{ff4S?bOpAU77#;&ht~JPA&R2
zQQ7^_hPnmJ{uVYLD{N-^EW7-2b@+O_UoRHN#>Tq$$-FdAG<LfOt_8TAI1l*#)c^l<
zf6T6um$&{kv-7{ZyW4zs`l~A|&)fh1W4{0ASy0h`{h3hRr^)kOdOL;H)2z6;xZdsg
zeC|zR@G_sue?Oo5&$Id2EdS?0yZxVspkg!eP|L)f!VC;~{B!5dwf%psm0SE?#pB)=
zFP8bvHd5Zbd-wlG`u~)&rca&fD*r3(|L@Q1cE39{&vLWg-?KGZcaDGkx6^u$(V|$b
zw!)W<`V(!79@Njj|KTHF!q-<<o&PTS(<5oTX}cMR;+AdO^ddK<oSLFJb?VfawZd&F
z9_I@hgzMjbuhY+R(Prd%%x?9e<#5^8S67!UTlRhb|G$&fd}o>E-dg0^ZKSOf<1YuP
z3c-O;`Jj<~-<PHOUr(ya*i?L2xqRNL_3Pzd<y>zQ_44v^;+VN#(EZWJA1$21TwH(j
zB^Vf-%dW?kPyG{}nwomi-rAYt=KRF}C+Gh;X?)&Bc&qVjZV^$@-lBtlrCN=O3Ol9k
z?lFdhhVnk{V&|8$srr%;8yh?S-??+=W*Vo*?W?JjDDLZcd3kyIt(76Kp3N_t&)IJ7
zX*6?=MPX8Qc6NDr`LfOa8<Ia9JE6}FYJ5CbwtcWnCwAAD^1A2cn~FIDl^%jhoA)~&
z_Z>Vh%v<+%>-Ar@Rc~%=oUJd)z~IEeoF4u7W5qM&8)s*mr|&D7cQZzzdg-sv&(EjN
zum1Dn;|jz4dwUwQPw;*H_j>=o*9%_#+wQic=i>e&yz=i2cs>~D$M3sy^{Q$5xjFf@
zucKd<{QC3rv#zf0O>d(Vu^egE_4|HheLKB2db?Zi&(F`rMMQWWXY$)5ep`5O|9zwD
zT#6=UW}p_x%aUIY4mLB-F=_~`n9jt&V4%mL%3z-T^3qaJ;WfqEqT)k>xA*A)VNS)C
zd-v`MIL$vRWz^N(?Y<y@fq^6O{=QoFbsOvd|J$5?-i`U+6aRlp&THPB&7D#oUw?B`
z>df^^wr|^(_N)Hy_4?WeQu8Vvb*?D<`RVD|+2;HPpQeV#z1&{!a!MkiPtWkfW7#%F
zxo69KXD<WQSEr_A+;3kQyxh#ptgNJj=XaF(;;mavHzyw1n0(x*u_q(GInm~Aqxu=M
zAK#<D$<^rU>Rt?DpBdl#`s8D9fn4aQ!6?DFb?erCyI(Kz|6a@Ax@C*a&nJ_eeN|cn
z*2V3$I&jgy_EoTXpezH!pF5wPo}M1Fc-^{pE0@pfV*B$+|9{c^;D2BI|C?Oa)6l5+
z^YQrNvhwnB8@c%Xb+x&<xxvA~X=jxCZHnHOa&vKMX=xd$A4}S}Y}qo;$!d2u+=|fI
zmVbX=cDi(&Wb^dJ-aLOQ*yC>}$UQqXRonROr%#_W8aAbzoaFQT$;rv>RWD~77C$@F
zE?2c;-@bqMs^6b2`S<<*zx#&^=6^dN%=h^3w7tR3xw&s&uiw9~{{O!-HrFQ~2RC*f
zHw7>m7%j@*%d<!99Yfti@jAtp56@=jf4lkj%i8VtxXinB7#I$moo!cJbwxu%W5w##
zr|qs>xG>>w(dA{nH{~zq-P`l?!b0bn)7|^!_TGE+v;IeW{fxWw|JHtN{QJV)eplI{
zh0`v&f1C90y739`n4LvWkA44OY+kdcSts9szxdvs%FoZw_usmE_iiA=(%rkQU#(od
zdbVNmu>;|fJ9-Y<9Jc%R=H}&l_vTgJyl}yw`rDfowtM7XFbMN~{ag0(($Y6ajxBU<
ze|Bc(<>&2N>;M0IJ$K%Eu&<Km2QV2h&8z?SbK10N`~TkF|99^FpL6g3y?bBwSbX0y
zpP53Yb$=eqr(1zKbo*^>Z53Nq1~2amS;@d~Vf`EZv+4D-q5}g1gMu!l7bs1<u_@Ji
z{<)u@pU=0iuX}i?wf5I?-`QfOVe4XM{(MyNY|0z|Dax#@tjqo9%O$uspJ{Q{THNp2
zEyin;@L{^dt14#mnmvlI<<_T9;QxFzJbv%Ce*%ZTy}iBv^W67$He6{uE%C&-$Hc_s
zOJ{t}%}uG_-`%}w>-?+nUcn(pyG1>jmzRl(i_f<xT;%`%)z#I;2b(6V`@8iTeucEN
zx|d#lS#zDCC;$HMd*936-I;lF>%#{JoBeISUYY;@=Y0MRorlxyKF_=!SN-<q^ZENf
zpR+ces~x_s=I7Jt&K&Fa{d#qzQ+V>pC!ns7`Td&WQ&Tkc_kOvgafgF}f#>ScZgKti
zeKv`QS|&}JRAFN`BRFVfh`;~&soLRQUS3_@-T4+4JV}n55~mo3EH3Cwl=yS_&D*zs
zU;5WeIn6N13<?aqxMJ$`%COtp^Yioa-o1TW`+RP>n$HY_)KgO!>(*SazIAv;fyevp
z+1wI4&Ye*`bOTa$n2A2y&5%&@``-6`-^1&_hTjxyxpnK7z@h2;zAR-d;uP19+q13s
zc%Q8N?;FRRI804VZ~m?SY+i4C$d#Gl0Eecg=F9Kh|Clv2G$bDUetMnnT%e@0Dc5v~
zwuAO(7A;y7cv^`k=2*hDM(ID7_dRTJuG-JLxWQoBJ#dxjyDgs~?^OBz-*=bKt6Js%
zW5Kp<-!8h#%O;p!sSAHO$Ex(z|6lR{wKw-Wr<|W>yZ=Y`ev6>1(2o7$lizp`%qV_(
zYN~elx+9N2b_%OM+bzbx(2%&RY<ufQGeuC(j0tJ5Fy}sG0D*zwNlO6K5QYiDPMlC-
zbJd9+Id^s#-VZ$x5i;XFQ$yj)?e)L6g9jhN8NrGS+UD5R{>uONEdNt1YmcrI=YaxH
zZ+p4l+*9SE1<@<kt^4+xZSk4+vOPxoOALPfxpz?Gvp4(VGy1YU`)<Get~b%6q@<*z
ztSs-w2E{4oIKKbM;!xD7vXHqa8D}%;{i}pCg0Laa9gQwZ8;hTx>jQ)Dn>e($we1#v
zvRhB$HGG8h*Y%4R3m+V4{4VCI@z-HTwPVEZn0?(%A6;OBr4OF#bTKG<z`1*MoT8(n
zl_qkO>vELqbNKFddf8KRo%Qi-Sa1}~?rP%bbXl}^?c2Ssf_ArL?D7LyBn)9ebl~mM
z`Mgp(T?{)Fr+~wQf#J@tMG8(FiHF;G-?IuHTC?HX&YJ72+bv+h{$SeG^NYZ4&%MB~
z81AqgOI3tAU6i1X-%;(Tu~mcVNjuyb+8Z>OdK}@_T3ow$@#AM#K`(oOcX0`UEE3kR
zFlZ<|wur+ZjzbmdoX5^nTr;w>y;X$XrwYz0YrQ1-a~<~%J(zP6-XQ#^aL9f3;pe(t
z43C>(_HJ*z|L^Vnx{pW2KgS+jyK>jADv@<Y<?rqYI92@isQB&i<DW*1K}^;EJLY0*
zC403GuJwinZNu$D-TM0~UR_yP6MnG5$SRS0$5zoNyP=L^cyLVx7G4->&tTc#d+-0o
z?elK9-FM$wBDL;i!k3?Wn**Os=X0AGZ^a-Hob<B9>b@z%WB=09uPcL>e~N_+obw#~
z`fo?<)(88Jto~p9<@HWU?`b-q1|-Po>v!$?wbkF=^2WvX@^iJnJNN?@^dG7bhq~wI
zl!~wQ?OV@3-QBhD<&CQPeclWbBCsIZu{35+#m7fSyA@k9FE2a#_+!m<<~!H8KE*;c
z8yJ0ie_y`)=rW&~M#az0%rs6H6BE03{kr<(lQq{(CxmESy?b}A*V0APO~D{=_ntjA
zm*0hdOGwBOf+p^U=!2L2?PJ5k#hq5HSP`+es`N@;e!jkr%;nX~=UWyp+q37-nVH5H
zPd`51FK=IG(=x$FZTIfo>~`Phmfx#<KDWFmPN_in*q8qdkW9(Y{H0aDMzQ6>#fu+1
z{cDt7?X3CvsYlXy*Y4fl-`@Vd^L(w|`m3?mpPZa*{GV45)S=2clC*KpRhAwZZD_hV
zuw(ko#S6D>GrJra6SJrK`@1)Z+1J)YZcg)^KkwJCUm4$9y2bVFzFcsgI(6#BeVnha
zuRs50m4H*l->=tSzkYrC^l5Ev?d<Gqei@4imL>1*?ELlX*NY>+e*M~<em;%Q)YNos
z{QiFj8ktQ^O)u{B_4S>uAI}Hs75@75OH@0oMIdZVL|}aU{i|2E=HA|xdwW}Ic*T<w
z6SL31z0Q!z@KeZp`Q@MM7$xdgtzPXv*J|o}EpSvMaIE`&{B4Ae*fl%zUtjM`S^NF_
z&l&nty^cQqC~)Y{&(Hh6uFbz&dOenvwRZpfIX-IK+}zq)T0)&J%QkP?v}yC^&#$ho
z4qq2D(<C$KN@Q4=Snk{jl6&g^|I4uNc(XWqTTbPhjmP_Bt;^o;{hs!{UBGX?oo)5E
zH!r$oW`~A^l$4dZy`BH>%krZ~kJe6ewy6E}<?r|V@=hKqLNl#*-P;-X@!|u4Lky`r
z?`D2mjF}ku5dj|_5OCq%Jvk?fx<3`JuC6w6>sPM4xiR_p&f@1r$%|MVGBY#3bzVFj
zq4RC>{68u!7cN|w;I~{qW=FxpL#+Z%>FLj(o}M1OGGvO!jjh??C0~wo3M;mhzP<+P
zYj-hi*}gq~M}gv8zx@0AY+DS{&&`>zmDSnR_2};R`(}HnfLfIG|NaOZdT#&!r+|~y
zT)V0-85|FfbPB&L*>#`c;Ma3!r+F=9(_lWK49^~-I?>zKEMNZo<8gU?ZEbEv8T-1L
ztE)nd)6cD0v*yj)x3W)Eet&;&U-ZPoXr_*~_U?j*haMhocjibxH-}SECvH#0(^FF~
zZm9hHZ0nXSSFT?D`+(h^L-EShtGdzKbgZu2xf2r@DCqR!)vHrJzyH3k-@jsoMvFlI
z@%XBjO9h;&tE;z|9|2X=0#2TuC%<}0J-E}e;*j(-#vVJk16`*~nUZsB%gaNp-0%1O
ze%C3iE+sAP%rSTF+{?@T_dnbj8XBtD0vUi36%z{!3ww8G=jX-pe=`aTKW^Xu_pQL8
zE>Z1G$IZ>m>YiKQuYA}l?!@u-?c0qy7wx!b+KM=F1o|ljg@jD``>BnsI4|TO?`a0N
z_=mmi;6%y5kbHYz?eBS2uU0;L_AFy%!=gotG7Ou2XPcd!W4U?vZtIyoX=cGM!W|tQ
zt>)SlJvp&0_co}>dGqGZtgEX+SBLHWdM$doUM!b=bnk-x7Xcg#*R5N(e7Sm-|4gH$
z%a^BrdUEpO>5QXI4<(=MY}RBx!431a#P_r_GcMk{SND8wdD=Jq6JOGAy)5~4OSoU&
zUe5H~T<hr`D(`l`-*<O+`TKi&e?L0fU0Pb&&L^9dogKX?<>Xd%GgH&8+qTJB6ezT8
zNIgC6xLmbP@ru=}Uw?UddFE3|D{JfO>gvV15=s|evHrbw`MARQ+7AcWudj_ZFOYPI
zyRi>i5*^!XS^TU;ASycg``g>o@6NB3>^}PFV?|`-&5w_dgG@Mg?b@|7C;or5|95<m
zYd5H8^&-{N^W@`xd%3Ole;)F;3vA21UH0zI&K1jd@3yx3zVG)tZpAgx+t;mLt!*{M
zM{T+9>}?qr7p)FopQiTT|Nl$-)LR!-wY0PzeXLmad1`pv&U<U`ez<ey?fnfr312wj
zrJH%_t1BAL8)|-j`t<1&sME4?<;p#KYJPrts@Re~uX5Rf1qv+&;GA%zQ~3J&`1s)9
z%U7;UnVI<Y)m3nI3JnWe7rA-awr$(it<!7SkaKg>q)C&G9Xlp)=*&!GV`KN_BBG+M
z*XNq$M%C7Wx>|{ei5v?TE>x_pY__fbHf8#Bck$Tmd9fT1cRruj%@@0~NL9;Ptb49)
z^|xEm(hu(RI8I~i*>@vB%$4y4sJOVYGWh2@NE5B?K~?>ozN1Nrx%c^Z9^bY0`}g05
z_y7I<y;#N3Ei&>Zi1n4niDU2Hy_c8!_m??vDCXSU^i=x(54ZI6bfcL%8X6X5Z*ELd
z^)_;M2BqVfe-*B+kKbSU`PrF;Z*OlmH#c{Gd+pjayPs|A>;Jx1Yzfg4?LV%qrB(O;
zd%b=2fs2dX#r0w&wg#`fa;TMiajeR{s@H4d|6YyfRowII)#{5oKR-VozcuUX*K=!E
zuIyyvk2$s>G(3F&x0mYd*1tZS+nF!_k&)|*;uHoRSnKBTxwisd9WncAem*+d4Vg%K
zx%sT)lI`2gr+SHsin^Xw|1!}>ZGFs6qr^ikb?5ii{4C;FxMRnUs`gFG&$4hYx^sZx
z`MnDd&hOpUlPK)fFgtw4dq_|*D9qQ&u2bFk^<Pc6p^5i29mxlvOdZ9yZxL)QzybMt
z{CgHV+Q_{(@0&2QMxvd~GQs8C-kGZzl>OjsnOig8%bzUg+`cBmvVC*<`DeQY6!PH>
zt~njYAM5Js%Kv%5{yCOof;p^_=CIDh%<S227KQnUwgN*zHmuHPK+NGc+`jdw`q`0V
zvzK~LfA)JLBac_&KG0AM0|P@t&9CK?*U$g;b>`1?Ng1pD?^C|7%D~e#6?08ML!thy
zYYSJe*518&&mJ4G?zcBKD(~LBb*rgZx9#6Amzn(|f<ZudYw*gD^73-6sjt3k9AJR?
zBH>PUMh1tXj8%!o<>HTzj=p>MPT<hBwb7g1^FS5FhJuHO8r5y;{`~m(`1s;)GmX<v
zoj!ejj^*X|fenlt@LctM>iQbhmJe5g{Zp%q7+cSsJC|Z439?W`ZN<8EautU@KR-X;
zzTPg&f1Zuy=lj|HOgs+IKxR0$H~#0;cuvJF+qb7zEsB*E(~pyJI<Y3}!!han7J+By
z=gapWUp%k&_qVIpu7!cYwJTTN+}^JL>OQE2kb665Wk_x9-!4&YCl2kfH7nMy?=S1<
z=yOw^J!=-ISh*fwzjqH9Y}!=x_t)2(Tj$zTZrZlZ%xYKK*;&QU&t1K8<;%@XrX8V2
zFTE_;m(6f&{;OB7f|vU}-3?mT!oVQW^!jfhD3AX7I&*RS`Liqj{ocOr?~L_(_Uzfc
z`}d(%?$p#&#g+&iv-*F3>i-<CpJQ2^mYkgI?d{Exc)U+`b=cah%*@Q393DjzaM9W(
z02)u-E9b-^D=Rzy#}RccZEa_cP8TJy?$r3FUtV6`oPPeCjIY|{*4EaUi$z66=gyt0
z*m6en`}_Ox;o;Yx3(c%_`?j%X$E#PbuG{^6lHe|X&+_|QWjlvPyEU-J_4X!rxk{Cm
zf*&6q3LLt*Io+8<JAB<7!(=xOMK8G)fkg`rq}%_#Dd1!@v*%ZJRh3@MjtPFt7caiZ
zBl)uA*Ncmb6<dVW{gzyQY3JGGlL9Kf=K84%b!tr6Q2zd&fRonLE+NjhZ{J>2Eqi-Q
zw0FKiqEnt6v%xlaf>`_DVY@t^;+CyjTLtFT{dy^I=+sngXO2Lf2ak?+uME+;bv*Ia
zm6eJu6*l{7e}4m&yGxd+%u2kp#8a`w`rVGj=g+*azW>{aV`k#tQ%iW$&dv%A47>;)
z_CG(zQdksJ>b`h!YO1!jmX_O-x<=4o8Klz9WBvT>EH`-epyfkvevNXALH4yZ0*4m&
z+etZDl)SjG|KIEVZ30qKQbx;7%+15Y!puxfL&L(3C2iENdB}ZucL})Vb$*YIt}d_Q
zk))0J^}lZ`w)EHkd3>|^`np)_J*ycW2R=m9@o|BHfv&Eu>FLj3uit;~=1oc62?9<(
zt}M3`b-J;;{JnqOr^y0`UR_<i_+N<DR8JwDxII5k$N&4pt-t5Oy?b^&-tL+$*|EDy
zG9@-MZWn^J1q)t#s|fwnOEjtZc{-j?@ytwP^_B~F?#Pr*n=)m})TvvyZ8I}9Z50sH
zi;?(h|NG5m#g?Pp;>$0;Y-Z=T`*tH4G+Jr9>8s0bGa2hLo!!N+udP+C7ERC*hW0uh
zOq=SVVwQJjN6E`et5&V5xo#u{8d1vBkKCkkx$enC_eE~5uC9fJg;)Mv?5`7X`tg9>
zeum_?lm2x@4JYF(pH5wP?ejwZT8)+s>F4Kx28C>Ye|!7-%^M!Yn4LwbD>gGeo(S(p
zPSNPtQ26**gihGjgo|zat?zt%2k9D{O1<0ldfke^u&`;I@+|`UejL?z>&?(?J<-f>
zC*gEsbGm<R?cZ6Z*-0C3e0+S|nPYX>+E<&mIf^Jm?5nBN5wouT_U6HBO+CGPn^L_w
z9?q+Nw{i^##O$l7EXg?3!s)!<KH-n%nQ~}*#DGbt^T@fm*49=a;PkSm;QabF(dff>
zxj)5TU!kR;5wWM@<3i_lx#~9?-~4;^>eY&gd#4=QQ~5b;ecayM+uOc@!KO`{?(QzX
ze)Vdrz_)jIqobl`Nk)c-n%4jOqp>L<Bt*pNgk9CzO`A6P{N8w6E_z2nVr*>e76+%~
z<i|(5#lOF~S^U|YyN4TAQS^lWWf45|c9X`Tub|0+Tjwfl=J_pO|Mvf{>-+yLum81t
z@7}$vtgL~z6lVH-`}R#tTs;5Q7ER+d3l==M8XkYv&hYc51q&YBnz=dKqTs=SHIbVy
zy)4n1s->a9p?Kx`b^g!y8VX?}9t|H=CVHp{buzQ_ZP~h2RjBj(Cdi=2L*E}8D>twH
z|L^NWU8C#MbRvy{f<;9{Qu4oCc>dJaOKo!C%65MFeed`EcH*!udy|o$&#!o9mTC5q
zbxb{}aE~zm)#>_h@!5x;f;Z+qywG{s+@AAq*YU@<x95Xe1dCj|H~sFHvz@hmQbV8`
zw1_$&G1Wr_6c^Uk<~cV6ywY1P@we<2zw>OjagX)cvuEf3zVjSZbR-||+xz*P_5Ls3
z`?u8n{UsRF!zl3>R(nZA|B4hzJbJC*JBLn}L)@8Jrmr_1mpff9>QDzOizd_u3MfwW
z*s^u&UU3(W%{pBS2KQji!3px6E=m(UWGo66tX-=+(IY~H@xe7^m;+b5N-uBgH)}ut
zT({%H$M~P+?C~~p*RSS3KS|XaG{9RE&RCEP3xX#xN)unac=6)Rn>!m4oi*k{MraJ&
z%F-=<JbJnPe|qY^)e^J5*KdQ=<FNd(^V0c69A)Ve-jhzI7|qnt(mGYHctCfx0&GcO
zfw!lMkhFPT%&wA59kJC>8!`^t|Jm{QZe-%y;)R*v5nndAIK6vwW1}*-30XGb@$6Z%
zxVX91{pReb`}->~-ejs*X?gkbpp}W|3%<R%8KJW+`FLOAd4`_UFj!%brV1N?+O}<*
zi_*t?iy$eZX1d$|bMtI(-?}B`G{e6B->IqE?Y8G<8mB)yGqd>DmzOs;r}xK{9$4!P
zt;ib^Vc`|cohZ|>E$7C*+Unxs&u7mVSAKeOrk-(K>TfnkDcz=UYi{D>W4-U+p1bg=
z=dR?jWy_W=Q=eDy=<D@(eK1(RW=+qwOpm2OU%!5p|8=3=i6eS@-rDu+-@kgrr8vjB
zeBHi%`!;Rbv}ezref##EHxpRD|6f&YZLNuaZtmMgX7)vuVQV4;ogy|Sxu(sN_A7mJ
zW8<w`w=!0lW?!3QQ(06Z$1y=19stJ{xeER}`;JF@TbDvl<jEAHTXPfd?kZJm>FVk_
zckbNT+2-@D%hxSkx^(YeTd8X|!Qk%QzsIEW?^M6vtJtz>)24Ut-hKP_t*fi+%rYL(
zf{)d!Py5`vd2?sV$w|7=+ftrcmVi2NyLVT=y0Y?kpR9Dks@mV*bmRBg7+;07Va&Ap
zSx&&m0xs;&2y9vEJ$>JHjYIZ-ECT<|7IT#lPudu9B_kulWBKJymqnX58@E6zW2?D#
zg^ye~5>HRleOa<=)hewnLCnrRw9ejL_SS3pWpTZjj&q<+h|x?NL^r?c+nbv!L$VrH
zPCA(~#p8vsXaP68!gwy!=@K0segCfHtg^l0E<2vvKiTcjmUewzY+&-WHIbsCqId7!
zJ(*$zS{o+hwBvD~wG*V95ADuN-9dEcyZZE?egBK_zW>)RU-rDOeSi8)+2NIbbFHqv
z57}AA;UISt+M8_nUX~8ZTrXa1us*E2c+nyyx&HLHY15{)t#K7G+V}Tue!IY}+qc)-
zUb%Z$RyKBHl51IgG8mX7d<eD$*O2kQx5l4-SO1`qy{#(#fQ1q?mlw=_rQRWs)(RQ}
z)N}Kh;<RI-_jEni-i(jUD{NCwO_?`u-pR@8?(U0K1@ekky?DC+&*}Xs^NXIIntJdI
zPfsu``){8JDo9lLWGpr$9&QsjR2|u}EOudu+dOUMdbyep2TNaH%e}ejDX4HNDJw0V
zYgb$K>&r{;X*xT<O<A^V*|lrew%@DD_6XLE-j;KBS846Js>i+NZlBT`441*nf(gZ)
zE=iA%^?Fa&Tf1si)UTWWcRZA@lL^b=&7RNgyZg|hc2%Lg;^Jp#IAQ(dSDtE<FE8_r
z-jZ?g>gsTmUiFKN^78WR>}-z2V?E#RS+3vn=~Tk9Ul-c#L_l@#`ID2?7cbtsckkae
zxrXRAXlZp|ttOX}&8HK}8za_a9bCj)Uk)mX?D>`M>{YnWI!W017O&#1+=VMvy!d!r
zo?p>tw%OW!`}S?ywr#@(1JH0)KuAbPP|&3tHzb^9SQbBX=C{3a<w{6goSaj{-m21P
zMW?6h_a9A4KR-{lV%4KZk6ye0t<Z>%-=B7NR;RFf-xp?vV`8_>Aw|d?E2G`U(5kE~
zJz`@h%foMmiHBNV%W7$BzrMV=;NhVa8#6C2TeWJ{vSrJ@efzc|{p+i%@&A8aUtH$x
zeY!`|*ew6vp7itcR;2&__V(=Avx+UNR;~K^`8jx;H6SFU<nT(f+*>BKzrKX74%@nY
zyE`+qB_Oz?IKiiDDP;D`{_54MkU0;qv-YS@^mub;=jO#A@OZDnp;&HK!Cy0euU-G;
zH)w?N+qeH~WWE;Vo%(G*`^3fLHQ?b-5m8Z5kuRH{Rh4Ynv7=(j@>A3CwCoM^Tpe2k
z4z=-0gW53Pn*yw~*VO9WFYUG8y5Hvc-}$z4-)ps{ot&g<Iji!)lmiotSd=Y94eDSu
zUXOHH`iWNYf8YO4e=o$47YR$h$+7mkjcc}8*2D=gHf)4-fD?2!Su-<iXMheZHEcZe
zj)&n+60G+9bOyv-1?L`k$CG?~?^SurZ@08Tz1C?uk)Qs2KL30-SA#UHD3MStOP844
zxA)sAd7T@x-fOiz%a+al9LscI0j$D)GVQJ;!;W0Ii}z@BOz=>-b?erfa|d4B*qD5L
zRp@FfYwPs0v!=d}DpH;5_4UgapP5EW7c5X{2?!1rcB=S#HT-j|ctI@8bIOO`@l420
zKR0Klar!yXdWiG$WSu~xtoyW?@?ibY2g`2EWpq1>G&a8eDT8^~1lWkevH31a8w(#F
zTk1XiSklH-tlegR?mp}*PjB6NzI9n_;>E+$+!Ay$va+fkb*jhJ|NROYU;6av)TvWn
zUti}p_;+Jt^46_eFP<>Xz9wT?v}D`1vj6}7Dz=z{2iMiAKg=@Ce)alw@XC-cMfdmB
zZr!?dDj1lWn*RRww)<BwxQ{FBZ*%eL)zIi@>F%S8g4dbXHEdeE5?<5Q*4BoHUx(Hw
zCw|+Ozl*7@txZi$t*or9t^Ip*b9(+eM&@Q5gXnId$PE~dhkIQG-$iLO8N{9WdoBFe
z-^XP-V!xYNquuw<oBeU?`}KRjCGV59E?X0|b=Bt0#n;wEo=n-a;luOpW#Ql<@iqY|
zY3b;wsC9D-?(eGwW$at2si~mJ_UhBqbpQXf|Ns2Zp+j9=T`o!!CF3I_Z{ECF`SsP+
z3Y&XZu3Wiw>(=uleSLlYHIIaAPjM*5?5+Ch9slc+!tw8#)?(eM{-wXZy!`!szkSOD
z50$)%{m1)ccY{_S#+Yp8J@PdKI_M<vAzucgV{9<t%p#7^(9r%hQQIMf=T;4-<d|ol
zf48-@l~~R7P+7$P^835H)o&l97)h!MmA$<+6*9>F=t$>_3$o>R1bsuJqGr{4eK8X2
zR-JrOBk{uJ%f$~5wU$`PwleU^+s%n}<5;+O@#E$8e=|7}4Icm2wVvxIZ(C*3Qt;?V
zXNlF^x#|)>j)Z85stVn@eS2{?c$ft=9{<oTAz=eoz71pp$OE&<9x7$;?%ccx0w3>n
zaU_O;CN*C-Dc;_i9e!noX?EC_ItQbfJ;xv4OkSk_T^~}Z2A<zj`8g^oYSye-c^36H
za`*Su{$3Ng*~`mo(LB$|YOMnAc0BHTlXx>*h$rpO*Hd-(f8P^#Ix)xk#P8l~6)RS+
zj^36NS+Zs4&dS%<*3O+fH#Rmlc)8!+hq+9LVGVbIL+6_~6kB+u%~~!=F8sp$@Ez;I
z%y5ml6W?D-n4zVmwdpha%%{A+c&fkMOn2tk|1Ep}w7WkKF!S?#(mH$g?D0O??5wO`
z@5=W}ul2v+9vK<=Z55ljjQxpC?DdbpgYHh%)xWn+E-5ePR(0@$l>*z>!|K-y*RD<T
zP`MW^l6W+1woyU$i_Omu9r*V&{-5gar&G1V*^0xork;AjAf_Mp=c#`E&sOocjQsrb
za>m)$)+qMM2i<J8{CH`pci`f6=%eX}cZ(T#K#Pn6xy$F@m0VSCT=UQV$!@(9{i&y?
zt&QIPZteDar{2{(?loWJW;J);w<(X0_q+FIOl<|N!?dX|c$JxzHOnmb*4K5WS3W8~
zI>;^`vpw(bG~H-WuhJ(yJpB60%gdMh&!2bMn;CACBA1eiQ0KO7+vdMN2bnbf`5oT7
zw*S#=Kf^pfj$`5W?d2yYsb0KTVG|b?R`%z|$9wngEn2ilz$rUB``I1&zc1VsTlnif
zG|!v|8gpQd$zJjJ_t#U}>r>Y6D{zK43RbSvytPJgrjMIvNBR9)^Ie(>?(jiMhjl7K
zoqc_Mz0&4kp&*cZ<5%yt;FTMavNLXp))i+*d@0z+-7?YqM6UQIQE_qme?J~8wp72{
zxqRQgy0f!Ph1LDe$YB{fxBT*;nZHfo+nbxhb0xZuE^-wKTN?#xl4ne9Ex%WpKH)KY
zWjUybfYkj@c*R6TUCk|WZ*BSY>({MYw+<aTB;e%i-0Z^iWEQj~Qm`7)EeZ|2dg=VX
z<GlYrJU-$+`|z#jPM6nhc=52{dsB@#houRuO|bjkLf_eDxu8MK3-|8r%e}pA)22<k
zcI`TI@%g-JJ*)E1&(0q2m%qO$)!Wt8wIt)-p32qX>(AM3+`HG-YD@Y1dl5Qfy(fK4
zKm3rM78@EWs;c6|;p*zT6bzbo-uIntwsp%Em+Q8mmOf~*V0ZcZefM@A@aX)yo1N#t
z0n{1GI=x$S7lPUdi<O~^9}Y%wq|8RN6WF=8_0`#$KhC}%KYMZiPrdkkGEN$LdUoaS
z?wpxvygmQ^y*qbgWM92~>$^0_)zuX=v-xw5*V3x$YHlvBMXl_MrXM|e^!n~o0mu6n
zykPsiw>uj?7XS6(+~WIokC;VXEH;33GkJnVTsRbSAZsy}$|IIUOg~n!so<S{?c3kk
zpJR_+(23ln;%=NiiE;1Vy|sM28!kIKIy#muS#bXJY<c@8x3^9Xu!Rj0{-9pvM2mt4
z4&a3r75RZJqS3jh4s4LM{aUooTcUei%uXSv9VIV=_Hi>FGhYZP_7Y@XKs(*sN)GEn
zLPJ5z-aZ^){;9_zaUa^yXfRxrF90faN}wg;Z6`-Z$K>SXeXALc*}@Fs>4cS10*7L`
z58dT{bB=-G@f1)3`tW~q`F6!=V$Tj9U*D_mIs468#bYz0&cE?WS)SIB&U}+qvCGNY
z@OIX%qrRV5yfkLA97!o%{_?U{tcBCPsG!vgCtO<cbk3n=CzL~GMQE0Mv|G4nM%xV!
zR`J~R*1TC!!AaY1{Mxkpe7mvrX5)Fav**tIz2{He{B!q;pPhMH=W*a?HE02W>FG0P
z-jwY=z0R0JkRjnmCuCf|p=!;Za|XuJ53QI@Z|p1nHGg0Eig%jl&V9We|Nn$?{}o%-
z1217pdp8@oWSA@#Pzu;w|Lfku`?U$xcZ285wcoLHO*Kn{u@ba@e=R(1T54oy=<2Io
z9M>fo5^6f2R&AX=bLP#rWz%ER1(Xg{{knC1U)k@s+rvXmO-)0;>!_E0y&CR*gunjB
z;jMpffQpysZ8<yN@B97j-QC^W=JU)^4rCH&=a>JtdjGH0?`z*z3vNt4-gp1k+WWg6
zG;zO@0}q``cxkUUK6A!x$2kU_%2z9wzmj9T5wM6wqP+h9_xk$Vx!ZOBecXNj?>qUw
z7v*h#y;wXwrpU857(8O<oRZpcA;V44h2hul{{KJvceyi|U5S1r$FSnExBlKM+CuAg
zzuV<5SNY^1yZjs5{U47>XV2vZwZSDyW>$W`TYkUx``y`Qxw{^9X-D1v^ZL`(l{bq0
z&&YQ0|GmHC+2m4wjw5&O+=(r}`}I!o`P{p^zCP;KzqhOObu+(R#g`WsYkz;sy}PUQ
z^|iI2^=rlFZMS!szrVdb|Ng#O^?z$>e}8*>YwPP57Z>mC|8**S-={mp=VNO=9^Ko&
z>-D<bV)rZG@BMzS`u$$PkN5xoy}!Hc?XMq&#l^+)|9>18w^vb$tNC~|c6Zs=XE(QI
zU+=H`vRJtO(PV$S%C}puf3vH<zpvKc{_mGve_ffLr(E8ceEiis#uaDP?LKzCumAu1
zdTjY!^L-zC?@d13&R_q4S-zuV`TV+Hdp@7DzG?2~?{EM4jIp!{cwWXp!pPWo_uahR
zch8+^ea$mHzu#^1yNeHca(<s%Tl@d%w`IPw>%KR>m9{^>^G9>^yqeGJfA0IXb^YDC
z-)}=heSLjpEsLJ~`8<8vviOpVuKN4`e0n6mJ?CanXz2T|4-dE3zk6QyF6aK9)YH?>
z#MGR>wl;e8Do~$Q>Cut?s6Y2WYg6z4d1hYyzNh0%&AXk?*X{Xq%KZQ0y&sQBKc8EE
zZ?eDL&wt<d*B>=FZ}<C6^}WyYIx0RKWUqg0UibLy<$LA#YpcJ%vpswN?#&tV=GlG8
zj<>&DZh!6Ox68->GO>nj&}2Mz{o-Qx{r~^IcRw=!&&l~Zc6atves1Qs`*FdUf9@g=
zP)cAp=9RB`h}Zm1!BNrhJ0B+gzx-=={mGa8$Nn6gVJ&dICTJo<bLQn`zI(+%>-qX*
zt-o!a|MyLw^}83gmEUfr8-KpBF*$l$PG`sbKTp(GR}=?|sQ=q>Pu@KD)|N+$jW&y`
z+I%>`tgTjazWQGAlO1tq1q4Cq4zw)i%ZrO~6F+=7%>VyQy1jMTn;W&??}`iV?Ra=%
zqVis)hX)#&_h$b8Isbp<+kej*Ox3&^K7!Jp+mUU#xAk^BXadbS#Q(YkT6L&j`}_9&
zdGqGIi&g{AY(7plbh)u5^RnpK?5O|mH-5Uha;ZS#Gv`u%jyWqq%c%E??``G=%}7?C
zW>ZS4tFu!A8F+Q?&ntW~78j}wO-;9^o}LEU9&uu#vi7~t(;nyTdbzCQ$A@<NJGYBJ
zfHuhNwb#2h``+&I{PkTGFPBcQTFqhi;~;<CBjNckYQOG%Uz?xL%WTka&i-$7!$i<_
zq1Eely*i>k|No!o@|8~}p4QvFX5B<k@Ha%Iq%FVP*%7u{^k|BR@WqT{^QzxCH9Y*f
z@B7~Oy7`em)=!DA|NFJ$#^&_%v-bSFvMuLk(Ic;CM{T)k@Bf*)zUJw}cKLlT7Ipvo
zmj8d-_a_V6<+?iB<?DW&z5nUcY5o6a@Bi_2`SWW1zgL~->z{_#rIs#cdhI)f!NuiH
z(P`c4Z*P8rmMosH|Myw2@G!6WoNtdVe);lcRs9{%#Mk-Scg3s!m4Td@ur?(vt!($*
zvfXoM2KpR3vqpf??b1wrON-m)`)_tV=39Q@ZQ1)hpU-vt_+J12`_cCGdp>!+KMLBD
zaJZfSdtCm#JwHD@JPeAm>iR#RDUwIezi4%R-SF0W`?|~Ohi3?=GJL$f|L^VjKhLBq
zCH?;P_IiB%Uv|3>jU5%APO5LuzhBq&-YR!1$SDjwYwiWLY|=TsDTYI^>`m`oCzc0a
zK<V;XJ7`(0mffqP;_-X_{d%qZs8u{}$A?4Q;`2cp!$51nrtkmr^#1R?_oHimT=Lcj
z4S7wEt9sdMey`&B-12=N4som3T>;Hvh4K`r@%_DZ{rs1;{iWOfe)(7`#oTrsT(I5;
zt;IcheYSc2y=}SCE;b*J2!HhjSIcYmp4s#J-R}Gy58Hfao5fz-v+Qg6|Db(VqMzB5
z-#Ns5_FOI4U|e#q@_B5@MOTnhv)AvvR{iOD{lDj{f6M_Dq2hH9#g!gyJT5nT&K#SM
z-r;eTsV;XapU>6b`z7f8$%)GD_x`>sk1sxJn!o4cG5x=n^j(kq{eJ)dm&^Y8bs_qB
z3~UoL7>`|-&fjyfdgCL7%U8qWV~<;a!tB7UXi24-&+AWzFHd@Rz^S9)tomy?frx5#
z>y&~pc0Rwq&%W<{Un}_W+V*{4K~=`DSF4{tdv-4C(5Ef4GGS{0=G}Yttm6l0q4wJ7
z?QykVuZG9hemyE4|L0M+eq6=F*8Tr~@0S;RU-x}?`P*BiZ?|3t4gU*%+}L0H<`B1j
z&F$OwuKj*>=FFLOvAe&$z5V^#_I;(I%nZ_pf3iu$gO=*vsQw6A6Fohq=;Wi%>~a+f
zN>znQAlDm|R4tcOT68tb^~j;pHrwwvS{afbHF4{`_`T=*?gmhIDE*V-^19Clo7q9D
zqe1z$JN{3S3us>O-(~y1mshRw`<5vlTXONy>qQ!2H817=J(TyWxfFfiL#uwx!!GT0
zGiJ}O{_(KA`uzJHkNfQZJ+^<h_aA6m-OHuZ@74W&%Wwau;5~2oox<awDogEu#sRk5
zSH8dCHc$hV3(oC)u17$f3XtpjJF4IB{r)=sf0YYpcdd<UxP`6l-!sPNb8gH5`Q^aM
zRCq?>Ry_3db}L80qTKsCx4z$h=iZjg%kS+r)INT0DPCW2dh*u)v)*sL9%o(k<wdJ_
z+>dG7*E}x!_W$qw|GgcP|4z9Ve*fFH?Q1Xl&&N^r-`c+K>)LxoOF)}%Zg0zd&B`TM
z_<HU3TP5Js+Pr?d$U>ddn{>pTT4bw}w+k?y^D92S&3?zzle?B%SYMB?uf3~$1eAjq
z?Vnybc<|R#jssDGss|!&@2UK}xZiG9`uTa!=hxS5^!V{Ae1FwtU-PThpaUlu3XE!<
zT4r?pY%V>|>ip)U`k&=*8D6_iVR-y*lDA&vlZo!H=79@oi>Xr;oeTwg)6U3tByD<k
zk%7U_6*}D9yh2>b#bx!?u8!8F9M$hS+g|f99Jn9=uCN&nEI0vf6uvF9Hg4r$IIuw&
zBJOZDEG~V<%wC9R7$m9|f&Iy_!g0QU;<BWF33rPb8P=!<f-PsrnJ`t+X|d$d5UA}5
zI(t5?%e{R0a;g78gUY9|w_nLIw$0N3x3U<n1?szWpcK)I9$GQyY>|)Oy#Ada!?8Y?
zx(IkHvZM7ij}r?+!XIcuGJzLSfI6`-G{8Kw^F?gB#_3Hur&|RQe=Xkl+qv{SD}#Z6
z3uKXbLD{bs4nxt;?62o>GbB8Kn)TpISR2Pt7n50r<qBOVceyh&NW?%r`(iyL<f6C)
zzugh;O}o38k%5OD>SvvkXEy1qUX`^qN^oPrT*nRuh7DR!eaB+c($coy&fR`{rlt!6
zLxeU={9;l^z_irJY21p3`u1*&2I)zII<9BpnKva?r{R&!U@!|NvjZHstpbb;3Cd7M
zOE_vfw<IY^ot<GVa5M$vt{y?CyObAUI=V*_D(RdP*&!f!H08`0E<pwcwiZZGHz-b?
znHm`udNpsmac~C%gTV%<)jW-7&Ya2H9*yc|gB>vGb(?g;+u%V8?*Jrh09C_+N(>A=
z3{Z=|MoTK0%nI~5dZ$rC?X?^)L&H^QN(<<Rq}*j7|4o7V?}o)*0Y#;#wbz0=Hawq~
zVFK2%RUOj5NT}827F1d^qi1HMA``;_3ux?b_}0R~zyR~o)=4N{YJmFmDzpy-%ccx4
z2c~Blx}YYVel~^!34)O9P*C;P8Ilb_q7Q_iqCqonzAORRx(1qP8&bpCIC|4WoLCqb
z!ki%)@PMaN%TDG%pKG>K3=Gab5J7`w&{Ve(?ma<8i1#FvLAuQ*L0p`0<f~IlN5<AD
z6hRi`5MYpo8L9+w3y%~e%nt|%D={z_oPb7_vs(uP!+{OZz~X5<qjP%Gs#TyS2*~8N
zCa84|U-_Z!00xG(W~g8TIJ?CvGBGeHJ3u@Tu=%I}xWsGaU|{IsfmrFl`>P$29OQTz
z7_Rj|wKsDJGB8M(Lvlewsfr5&L&6+rpYqxa0Y(OfDySo~9Xc2o81mI1ArYYpvT@EV
zh%pazK*l^c1J&r~3{tvU8KU$BnD^Wd;<yc9W3s&=ydG1KF>B-@ZZs$Z1=0pXh>Zs#
z1e6#U4jh9fMrCi1p<ki#tLzHmHACZ(#};G^j~djNDIjAEen5@605+xonp~KDK^hN3
z?dOpJY3yNxY7{WqtaEzRDzJ7Vh<1luG3ORgRDdcm^W%U1fr=r9YtTe{BjB$ys8%W8
z{O%$H1E}#1s{0tO#io^6%>^Y5K}8n^h6Y1u`nl%%Um)?k?M1MQc0l7h8x;Etb37m+
zGDjHXq5>mGj67JwDagRkU<}oN&O;xP6ExHo&**Ri6$lJ-pvl8R6Rc4W;_L%Pprn(q
z7Mf{6QX8OB27Tb3;JYADkJ?bQn0?WVHPu!O3<)j}#R)NAI0e(vltH0zpb8Rf2Q;iR
zI}GYKfPHxplGz&qAw`p_3j+hgD^+l0Hv~?YIoD5J>DJq_G`$v%qvvMakN=**z;NIP
zB;F2eh~g29S{vq~)52l>?by4tz0Hga3}2uY&-lU*D$RL8t~Drtc*NjEm7-JRwA4s9
zP!elc2nmS>$1OT8NP)m34v9e?F0kGth~9wN;4*Tz2)KL$hvY>Ch&y?&s=LgX**l{P
zlKo{F7!E{1ow(~Xxb}#WR5~NSrrLyop`la-qN#Ku*mF?#gFW{?fq{WR+6@wS*~{vo
zWuhz(1H%F_NE}+I?-fXVsOJoFKf^SLS0Aj=<reHZsiC&$q6nzzw+8C0%aAe<6n72A
z&|v=es|8eU1w$&q<;)BW5=;=2BOnFoRuQE$XWV>_{y7Lvn-X>q*<*-AxKTh!#5%dI
z9vs>{){wAzyaL&^3DELgLVT}4;!8beL&089vCok37vf@x$eA-|-pvaK<(>d;hzZI|
z{lRGnoC0oCK&)+x1Gf#X26Y%1OM|O)P=S-M22wUic!A1gNE_FwBj8QX&%a;!*cccN
zC_#*LUM#Mpv&OmOv0~G19tH-6W-oBi9&m^jQxcgR=;N{&RBbaza6v3%^9EPvqM$H@
z`Rgk*{<uJicv;dr2a(C8{_G442ecqgSmFL&&)H>5)Y^^+F2QH~*TAj0Y)F&rz_FLD
z95Zj`ZCCDV=Xhttz`)R>0dd&t)#6GzcT?dp*24@D&(<+E-hBHlsLp0MHV0fBHkeKT
zS8%X+0oSYZ_!t-(svuz*vGiVS8@P1XC<01k<pvD!8oy!l1#UsHYl=rc3LK0DEs*&I
z$$<r0;P5X;ikgPYkQBQiXuiP0i#<tRptv~@2~jUl_y^?6+i%ar>^LnDYzWD_3=(1B
zfNS`86>JMMALhLlU|?W~V1~#UUV09Re{c+}fwZa)bi5XZXCsD)e27xPE4xn%2!dLu
z;_r-(gg`S4!?6x<KqR<)acXHT<q(84s%~3>8mg}#IangB_MH)E9AMXJ1_p)$&~_rv
zWJuk&6v;yj5MT1_f(6>~_YkYV!De9!j)wyZuUa{{Ukf~{=tz(SwS1i+5o#a+DeNxA
zwk0V^`GQIZ258|QFdNgZ3`m)rvjVGK5ycQ+f1Ns0@Z!soJE;r|3<-B2ae6HdoDxf?
zGcqtV6hdmr2x)K<Isi$m4Od^II?)BvAbPOo%g)mQu$0|k4yo_v9D5Yo2C8?w%O#)j
zxfn)%2OXgVYv0v=*?)RNa{ilbe4vK=xhde_RyYJsFcIuvS0>cx-b;<(5|lctu?$p{
zGi*=;mx>1}Ue-@vcti5@@9cS?9OedTA4wE~-RckvQN&RM@(2T4BiLjKL9pb31(39L
z_!QV8i2{g4Gj@QYkippnoK<;t{VSK$+kX>9$peUz#f}gqEk}3rfD_F%#TIZPGdKda
z)*<)T{nQ&9&C6DU(qn=MIA|MsA@$5tJ<vf{(0sRHpV)gNa30=u+5pr$s7eRPUYiPb
zZ$vUAvt-ITxtuu@1>S=57*a^DX|yYsgola2Ik3G4w7$MG(h1=xJQUmZv)R-il<#;T
zje>$@zwUr@{T6UgzkoCx=6L1HI=O5C<*ul91yBXZ;O7qZo6hDUStl7t?1X`~+A^d=
zJnaH07B9uNxfyeT^IJnPcvQZ@@apc<8;WF|KqdGrb&xmQTEX5pcJ$P49>{1J#1CtP
z!HK<R(>*J2!&6XwjWeh%!NBm^A8fnw((_g=9W!Qzf^C+D2s-B%${qsO_KD!E5nv29
zg=hYmo({j|n*$V&eB8|guFYO&fSNDf;4m?efi#AuT~h>QQ-}k2bRb&S-HUBIle7J{
z;PLB<;A~$BDk~d|!PQs;<0SB0!dhpLX>JW*zhqsDZWj<-8>aoJrIZ6yzcD170GG!J
zi@!#;3kVu-)&i9d2}i(jo}jF=&N*ZI?U|BCQ%q(BMu20m#|Uisn(nEu1s?r~Z3A_H
zxj^Y(f*YI=BaFLW3ru@w1XXPaQGL(zx+1vpJ9`+U_KGh!pwgG0w*vLje8Dln5C+co
z2QpLATuQ9wPJJzqr1gydnw1o&XmS8&NC)|uk<(H`L$85Db`7NAtF#8}MQDyYHocSs
z<m1=iI_rSMMsSKrPzHNbq7b6<iYsV@pi_O`576kU@ccjL|Nrsdd~?pMS*x;D7#J+v
zA<B3-0-1Vt-~V&=e$D5@(fND7_F2EX@gtFeA>k^-P6>xaEE3tEg^ur!hVTD%^=Lb2
z7yY-lx3gz)Ffg180hdNwAudSh6iQgL_roD>!NSwJ+wYWI_ANcn#K5q8I=DFlHIEBq
zp1;jU7neH)hj~l+nHU)UMuNlL;ED^U!L_;%&GIhbVK~r%aqWDvxAs(icDMU9@pj(s
zw@uu7D?m*pSP&>D!{;h^&F_5p^ZES#pFbb>+y6UPzVEZheG6OLT0Q&kpaJ&hcHegv
z|72lca61W#sL6((7tU&DW?)Fj+4K5B`MuH?3yc3uNS{}k=G)ddQPHd6W9gF<6Tk2K
zzE`mDZt3;E59R+ow7&PTFK>JG?cD9IN1mOX{r{u<|BoW~t!-`J{tsqkun@1krU>fa
z9(j5joFWr6Y~fP}1^XC}T?dWP|Nr}a|F>J&*5&Ws*#3MP|L-Zge9edJ^?$GP%iI0A
z<gFiDbW-)J9|yw$4a=8SEvNJQe*QJ(X9p#ZUjm@9qJ6c$yVS#DORuI*kNpPn?x$1Q
zTYFd-4roNiv|q^Ds`TiGRm;j7#r=>Z^4MVt!{gU4E-pTA^ZCq^Cr_;J|2!u=fB$3I
z_bWiDpzZFdQjVi3E=SI+0VnVTn@*vG8ql0radEL=;iJOEpdrZ%3!McE|NZ>;osprz
z_)>KHgqd&4tRMXdY}){ic7{Nvo{g*5?OOHy(fNOG&f9&SdEVv7>-GEV9`~AmlZz|6
znYz^%)Q^TXm!jH1b<j#DPzp|H^#yHJsb*VX2HMNI|M%VZ-{SIbZ(F-vSy4%@>cv9s
z^?N=+qVZ83D+7aDbefyX^2?Dc1d<A6osNJ8`xqD)1Oz`uq~G?J|IRHEac~)nME3n3
z$ISnI(*OSnJZkx+tL3!b?l%V-nXjoZFkG-cB_Ie%`x~;%R2fz$1~V}ncySUu1Xg=Z
zG3N8hD0gPiDFz9_Oka|NnLwS2sF?N(dE1SRyG~}9ES|9m?Cl0qa05d^S^K)8$Yeo5
z<INkB^WQ*<;tj{ZMVs=v?R<*Hx8If>c@)S7DF7nMK}86xOS5$g-$GFTsudI%2iU*`
z0MFj3?vkLwtC(FZ&jc*d29a2LFR-lx66Ot!;C_$<pL%RNxN!?fBn*&7G+TJStW(;y
z+iyEMl8iug&>C>@k)Q!C^Vm8TgX&|39wAV{<uKQ;YKkJH`D%I|l+qZUgZuXzE{a`O
z+!D2x!(9^O2#H9r!%L@wo19NIK-C9BGPv=x;hSIS6wv4w1H*wE-~w(<@6p9PkU)0m
z2NlZ+Q^Vv^J08Cl=$dP|fB&<GAQKZ9z|J{01>&57FmRYHe*`faQt~CtoV?DtU<#-}
zWjF^O4v{FIDJeMDFFmsLLWbKzh}t>e>dK%`2GmeuXy66c4#y68_CwllmEa+Ji9nDC
z5=_3X7l5f#0d=?B+9CEE1c2;6AOKCKl~Wdi+b(k8BXf^|M=mAQ{`MQ4Ipbz=3gmtU
zW^hf+v+L+$a3ez&VtXPuS+JcuwV0<{6g<WUj$H<4@EF3C^MRmBpFzSC<eCL+V9%(7
z(uj+Tp=58G2)K4)kdOsQCp3cw(K<V(x=V(DOMHVxU^iAl#vNZ3Pf7fj0I6=E74w>z
z;GXF1<xY@v!f*^6{1)on?vkMHB-G_P;HGYR`bMXSi1r2GtjOT!G=+hIVNM$)G!DE1
zMNY%iU*VuWA;eJ+biiUA!hJ@d7D)*tS}edRS|Yex^w1%zmZMyW$M4-(3o#NhdU4DZ
z(jwWkn8$dtU>~RgXJB{@woO?Z+^Nx9Eg*GP1LDO4ZlDx&;KhkE9WP{4GeASDB~u(B
zL6~C!O0hiK&%6l*wZWkp4uCaGfP_g@*_6aIGjOeHBLj-11Dv2BGk5`Rp6AWyItcC@
zFeEg9YRiNNU)Cd8Z2)R^9&k8u#&5YWxU;c%Mh~bpFA8etGZ^dvnH11}lq)e&_E1tQ
zs9g-n?hK2;IYwH=`1GbVtFmqew#~c=%{C1eL5Zy4E2O0o3Xa7C6`<fhFzaPJsJ>VM
zYJxH}6oS=VJ$<I*MTr%-XQ%|}w}*fm0^q#z;0+`_w3c!}Vx(a**zcKNRthL3rI{gR
zE(W$4pnMm03Nps9LI6Bidb2|lQVcPi19u^nz`?%?>;i^s6EqkZ8VaxMcgiqX900L(
zjU*_<lM^>Ooj6mrTUId8#|_kB2Uivh33EZ|LBbc3hLyywDMI|hU^We8(#450|7t_3
zFmS(^fx(~*6oNb(?L<M%W$;J>xDIAukf;Ww_cl*RN^xox5Ip+Csr3OQdmm^3#f!lO
zNPI)ORR<CvCH85p7B^!qNCX|208$rV4CSo=h4O}rzrw*o`_KRj1T}>d?&^j?df!kZ
zlR-`A1GgZ7FSv07A2h)k)PWqlruWR7cu3ps?qzUM6aem&ufDR}siR|NB*Z=jTaXvd
zUGdaznNbAE8jwC#L6Ej~ONWLUB=Q;ND1w628`9{7+CK+0RB|997!v3kz_vbig(Ug|
zZ6IIu?D!|@v{<qg67mgSLB+p>Da2EdLHvXl;EAwUsB4%(MYBN-)B+yXKqiI*0?=4B
z=mB}F*#l|;2PhHous}5~29;_7{7{}UsA<2!5o%12G{~5$K!`gHo`H<Xc7qthpbiS}
zHS-`YWbgw8X@MA2V*yy{@<otxp&=Ahq3J-nTT7v0#h?suY%Wy$wN4F21_^7Z2X9C$
zVqr-5vJ#R|!2@0iEMM0Pq@^iCE!a@#!pU$z1Cn*2x+G*kx-_7=5<YYaF*H1cmUayb
zL8HcOfoX0>HgPSS(Zd1tFAr!A#XtaRssU(F+QAl5cDz68)Cy|4NFF^m<M$VE`NhD%
zz@r<;#IV89Z;imMw`DFms|B8YziuNBDtQ<-Tm!kODsqj0*g3-^oT5#|;Ch$AfDI%*
z3mS@$a9a(jSGTSKWxg3fkN{>_#KQ1E2b!rM<8*$i#@2KFK=vT~eju|$h@nAsrlcU$
z;|abULJSFO!0F*`#2SIFlO52QXLt@89)FyQnxE4|tdq~5g=q$j<t7}_4dXa{&hU#C
zOE)CY8@{S|F~Hg_&^QB+#x^Xx4lXg^jy4riWjOHRB)AYgt<^H~CM;|X99%dVKqavf
zG<XkO5L9J2Z~&5%v6Oj`uKw5PH3HLgonC};$p6|8YPLTIw>=mlL<5-^t|)-Z+*SUq
z34-8Kn!!NJg_B`}pwZ^bFI_-!d$H$`CAi9ENJ!}vVn}$Kn&$FlB~sEk02&}aAO)EO
z-MgCy)E9Jt1v+SirQvD>*ae_5xt^E($KLT@g~sp!1_4!u1ClS-I;>cg1#XcfLS5ld
zzleq5g+4g4pk^_IfpV6h`kdsBqo5?$@x=O>Fl;TrjG4W#C|c8|!PsE9H9}D*x&@Tp
zAWdZkh6v$BEDX@fG;}_fV-q-$F>K&-;biDp1DZPASkb4b1a6~l-~g9Xkl8{B!S3f0
zNv^F1#;;`{LFb_B!pXq1`Y6}J^OkJTln?fcfeEBs`W3<<1nMR;Ff<%g_F_nINCkIT
zBMSQ#UhFxf1eH=$@M1^^(3zL~Vz0(&qpp*n!tgnG&YFRtZK4Ka!>u_LeGA2vbmk^&
zlpVXXTpUzFf?VA;S%a~mck<8wH+S#9VDtL-|Hto5Js(fHleNw1P48+RZxx;`${RO$
zv^rV4zLg2O+c}FfHp2XidA;tV^-oXpxgASLNmn`|XYz7i>fFSf(^6cz+1D27I;~*N
z%zWV%b!TqpRhc`H>A(8#AGA&G$-H{Te_dR;TtD;9cdvJO#k<Sd{@!_P{ZGrb>ixp~
zme7D&u&P6d!6WrqF^}PLrw`Bs)3DTW3c~^Mv}H~|BANnz$}~d@r-YY+stg+n9zTAp
ztgM`Ve%{>WP7@SB{$gOzsa(XuV7L0?$B!K^-oNj+oS*#TmTlqYW=JMtnBhHz;ehE+
zXf!bV=M7|HkpBp^j$xm}B9?}w+0PVOU@l+~5>RE(QM5Y_G5tW2kSarn@%&^6SHXD^
zOT(+ev*!(K9^dAJl%os@A3KE@ir!SQ!)jxO6SFiJ75tAKKmPmo@8oli9Pzu=^QJ-F
zTX0xVmEq2gG`r&gJ0&C~J7?#2bHI%kcj06(pYXc2@8PqpAAgiXb1lP@4h=>H<HwI5
z3qFLFXbd6RflLP;&HO2I*oIy42h4k+?o${hT=?<hhYh=8h!zVcB!e<s^_#-*p>?Ju
z$HF2*BO@X6`(lD0d-JWKG0~vvF@<46!PHL=-rU^$S*E4n^R0?Mu!?|Tno%IrfppVi
z0jD!RWe#6|*Ew6Z1ritx#hpS7>t0qN*<T|S$aJ8-r@!CbrM|x20%`%nK8{5!2``HT
zpw3ZI^<tRD#SSfU4=4$%GK6N=fWRO9`k$wN{P^+i?rvqJynVV#kEHLH!L;cx1u_Mw
zuB#MqdRoi_%`6NNC5u=T)~Ct07#JBD+1b5QfaDa09Sb{z8ftz_X%=`_t0y-Z8u$l#
zU8gYIS&$~*5@5-F&=K4?WMG(}zlh~Q+Xqc{#XJAz{@DM2@>4sA4<CFMR%N(3qgLS1
zi4!M&{MZ3YR0~dZ2sIo9Yxo^wqvi-v-H@g1#V~1xCC9^LPA9nErwIfyePG@33^f%T
z5K{GG5X<=f?k*@bKab{tB_;+Xrzs3N2S5SvXy*UV|7Bo-9n-GC_~XvGdA8Dm8=uy1
zmV<`noeR@aLAS4fov^|Wmfsl^I6;jbu;>KHXf%lJ0joqA7*xQ$CXgV*Bv5k%%wq5a
zH88=<1{F|S1I%;~!e#7_6Els|&2n#*ytuG%4R}b2fq~)b0<+hQ2YMf!nyL+2cI5Ks
z@B8}mRp3Rt3=9FhbEYu-ISyKHeCnvL`Q0gRA<L8*7Fcx%HFR_=cJG(7EK+$b3bvWS
zAyiP6!O+%rZ|do39XGb;$6t$vNHn|(oc)zSp_<?BM}o_n`}_Z2i_Sj_x5ReN6bA3R
z2M#!_45@mtusv*jT<ya{t!th`^fFvM?8&)d$z1#Tean}tD^*ojD=*sit@w}q!ms|(
z;(xyHyZ-A#_1b=&AN`QRyXN6Ot8dHK_q{z|#Tp_V$fT3`{cd@@i;PW0!OczQYx7Rs
zP}WzEfF!&Lz2Bqcf3E)*_;lmHYv23o3yT+c{itEqiEVCfp4xru)-B`ob91cA-+jFv
z-`@c`R^Pfjuc*lAaw)`V4(}{A8K<q=e!p&a(Z;&JzixhS2ZQZ-cYlFS1M3E@DQ93{
zaF{Bfdf?L6Bf|cm1$(`{y;JA&^77XII4tkAG$=d!HE7|~xpU`oo<nLShDvp>hNY_`
zbiRH0a%SJ>CzJg_OK<z_|NXe^Z$I}c<ajlP1z9?$s~8HHcNRQ!a;f=nkUf$gs(n?X
z2ICKp+}vE3H8DGlUPE)PgXxE#Y!8%=CVl+!@^a*Qhz}T6wM>RaDFaby2&#?RNpsX#
z+GYFEu%$d~w?^_ah~-sP{hC{TZ)e5FN9yvR<6*wOx>{OV+PxpVt%8AJfzuR*%B{uE
z&mBB?@b#X84-XE8t`5_^B?pOs2Fsd#j72XmEOfrE56f~!&zG>ASUPXsybcYq?)pEE
z<+)p(t`$SnGdz0;l3ccISx3dwso^4!6?zN@c)h1E>@+bo6%_<+PQL~ox@2HrSm!x~
zVds^bHzfu4|GxWv6RgH5@Lj}mV(DD#ayOU!{P!0w1i*@9hG*fkzcM5!9y@kSOB}L^
zqB=YK^&D|EH8nmNiw#LfyIki(ItC0M<Uxi^oH+5OeY#ZZO}&^M8$7aMu2lAFh^cze
z$iDB>)b%DNCNZyDzu&LFZwPg+!!{Sr3EsJ2F!goO&Ohg@-^c7KS-ElJ#-&S_?%QW~
zI~Qt@=pvRA8yCCx=iS_tdhYDCYuENxe;4Cw^-+VCYG)#*ePu{!{PpXXU8IW3ksCK6
z+})2qdGh3uKeTe4;X8$4XNcC+s_JU(>$}U|S`|Jza&oeIzr6jut5-vJ?ks$K?9DZ(
zVJ;m)3d_sh-;3R8QTge~uF}_1Qc_u2S+i!(E`M`lWBU1d9U9Y5FAdT>-3f8sfl2CK
z4KY>~9}-;N+}^%^*|M~Ai??Q9-?nX=+mWZIryEDWLP{}^=}Ay<aInjry5Da<+uhxg
z8SJv><1y*(UC;n|@E8;z!NI|Tg>P<btepSt>sMEo^z`S8@}S=MVG+pm<Pzv;B*EL;
z^Yy2C70Z47`gQ8*ryUXdYAV^lbG!@;40bAB4Kc6Yy;~QvCVqe2^K)~T`^~i~eRbvH
zV)yIo<L}?TEiL%=?OW}zH781xA%56!<!YKOLqtbgTbph5H=D0lf&~pNz~Fi7zkmO9
z%AjFWbHtN#LU{J|b+&bXe*FD@-~Rc(ZMnDo=G(=FhfhBVt$i9K{k$4vvh(ulo=y!{
zN~)^bwQbwCS65d@u|O<)a7@vwAuAw0US81D)D&Jd{piqOEPDL&`FwdfxxCz5St%(g
z35gdkUuueQU9Z|1qbF}yvm-)B?-{hQ-f+`}bAq&czg%xeXJ=<e#s9zGl^4b6#XtXh
zdb+;;Ina6y(5Y;dok9xcWp8hletdNF>FMe2M?hP$@7;_0nRR&2vgOO8qoTwLWg+f5
z5VMHoM6YS~wKFl|-LCV~AfgO;1*bR3G3Xq*a3R3m{rH6o7xwtiumAV6udna(^Yi+j
zV4em!$9wYR$xn+ai;9ZwEeG9UvNL91!NWsJ>!Hy%<Lq;31%Ea+HYq77rAG%EnWz2_
zt2<|AW)`M3RcX=Aos}QOp(@V^s5<1u>?~Tke0ln}xu6Zj7cXY-y?^7zj5-mBvksVo
zJht!4Qhm|xqo9Id#*7&n(@(cLB|beh_4B?zo9F+{xoLDi9U5_)RJ<DQL~P5s`S5VN
z_H|}vW?Ng^^S0mbfLaZ2VP44o$iGr?de>3VWdj={_FUr*Umq8H2G%H@4@y<HHY6Ut
zaN$B5uk^FpU!Ujy-?Mh@+7K;WL@3nkWBl~$%F1B<Js%ETxpGBaKD4Oi;>C-vWk82!
zgPcBLmIh<d$`GyE@4N5s{b&95?b~i+Xbwq;v%Al*^UEb~{ZPHKcXup><=@=en*I6N
z*{QGBt<yWM0@BXFz;IJY)j@CXzhAEfFJ8X<`O9Vh@87?>yHr>I-c|a#s;cVS+uP!T
zQ)fdJOM(=enwo+#Qo75Pt5-$2S}$I_cxtg}_O+069v0>AV%oT&fdtKco9oxb>@-r+
zi`^CS^XD1k^DamB|NFZBs65n_2F{?eWTime&(rZ?b$Pe8Xe#;5u?YNG`QSifheiFr
z8n@R_yOmrxCx8ON$%TiXzr3Ubbn@##cKM!;<@2guoi)ENBO$S(=G4P;*6+_$K-R=E
zY}f+Hh^E=swrt<t-`6+o^wZDJ&##Z#x~k}Z_5I)X&RFZ}>PoJKv}YMQ3WAuPT)(%s
z`uiKu+4uJ@To4dkzxUg%H_JhX+N@l8G6OU}4r*q$OcqkOer;Xs?(+BdPJQhX)%Kfj
zXKQOKEBN*6*O~92$*w~{)nV#HXqmuZ&KJn^<k&Lb*~g9@<7z#WUJ1{H%i<vCI5QOJ
zg0c)K|7lD=J^v#p>n)!Hb=wC`PzD8E90Ix&>->*TU%r%hdwXAkrKjuapd`7ss`Q#L
zG`#HCLB-pXCr|AE|2)5I*RJfn1urixy}CL){yEgXhAW^r`uq2<{{BCozTUsNzkdH=
zXa}Q#NB49U!v+CfUf%2V3j;Ln6dspFN?Cr8BJth3ck_xt#|P){|NHG&ue9=`lj`%o
zyoE-sLbRV(!>f%sH#b3#a~1sf>9oGTpI_QuXb&^uk_%@*?yFa?cJ11wv?xX|eJ|8y
zV2Rdywcqb{R6OoAhgY&XnTuFfm_0uaD%JWr&Ye5AW;@ik-ch;Nm>XIpq@;dbU0-+g
z$`zC6|E^pKSsAh_Mx3kF6}Hp&gCHnSjvP7SQ)94x-aI*2Z*YPxNPJi6YtX)rwzf8Z
zyR6L2$jxcJO3>7NVh$*ytpx#HSTBs>y0TY8jdk&}Gx_zuZ*RX>wK_)ce%b9@(C){_
z{r30n+?n%L59+=^rYG0FeECxQ{cgD4^z?T<y}h!QMN77C2W`BzwY8PZgT^yQAk&j8
zS67FFPVAjI|Lxniy}iBZ=jLcm^#au)H*Vb6k$k-G4I3<p2`S94`};aRTyI<M?Q3`M
z+7{PWRaJd{b~ZUV`S0Jq=k0#yEDSi20qy2Hfch~T1VC4dxRjKZp4AglTC{IpT}H-=
zBHQ1m!Sjm@3=G}=Qy472{e53wfA4qn_PpMXix)4t#zLL<p%BzSS-EoMYa83@Z#gbo
zwr<tE2A<zyU}&i95K@?KWo@neXu190m+JDq{{H;@{O&H^-lvm5>n}h-F~fBV!_F^{
zy7k?TOjLHCSO4$l%a<=_{(gRL?(X03b}K*n_xruR2-nQFQ2RPooCPIAWo6~5-47o=
zJhhXf)oEdj-uXQRPfkqy`uckKyD6ui&YL&y>}>P&jnKk8Vizb`*Z;nK|Jj|Hu__@`
zyq3;<z5RZj^-ndB3m6z0`h`^;^lZP~NQP}4JOAkuXftK?yPfJ^p#@`uCCHT3VQZtd
zW`$nUhUjEi=RSpDXSm*U(4}_X)AdZv%(iXWa^?E<bXbyBaGJufv+T=@i@MR<&g{H-
z>(;LF_xGF@PB;lORuWwQ1;oeSSMP4=kdT&s{o=)pIdfwEctWk%+$p3m-D_!3et!Pb
zs_W}w)&1sxnob`dAK$%emzKGC-klkDtsrU}z^Uck-QC*a>i+X$w&%q<J3IUO`W6>|
z2JM@Bb5mIGY0=IQtyD`vh=zzyAWyOL%k^}CP87Xz<xJ(4vbVRcUAxw?qvmH3YMUo3
zE32#HcK-g_d&_5=<#u(L->XOtG>7Jg8c=@7%*xVI>XWrLdw+S~zI`1p-n@yq;tR2f
zK_1kOtNL=$T^Y3NZ)s5GxyAeT?JL?TBY1kc{`my(Ts{NCfp(CC?%usSwfpPq>+kRE
z{Cs6)FsOl5_2tE{udk~?V-S1y)*kH=otX!1@J!ovYYjUC_r(hrHl&}Qw|riemJ%B~
z2t40&;r{*kN1-{cq2EH2@l)2@TU(F!%kTgH@AvaNpaWB1US9tC)hnY&XhCCe*pqWY
z^4YUzLHinAeEj{dU%M8!w<@$7>We!TbEYs@8e3T0xOlM;wEMX9daSIhY&)N9SI4hk
zzh=#na)oBa0_VaYrYD~6?(REd=KYM+n-1z_Y~Q~9{oUQ>heM&peJ%=OdUEW+g9knn
zFI>I4wes_`Q(s@acriPF-^{x(cmK51WGs5Tbb8#Xm_1*wMNic?1>GzmD=RDb@zv_}
z>moO+oi2nLH~+yBmJ_|pmM!Z5oq4|0qoc!V;fB)J*E&Gw(Ql4{mcs_}flN=%fNtgb
z{hXQKMnL`>v>U}>E)vM}<hJhiJCnXb1GnKPXgV5fJjy5qs9R<?|Ki1qpve34=XLyl
zDM4XjVNkZM`1t74)6?4E6IhPRmS4GgHFYzz+vPBYVdv)B-{1W2=}bQ@ZJsx$`16B<
z&7eFrB@k=}1H<PJKiN+7c6N5IQHP4u@cDT)+^KlG^}1fnjt@8Q?|(Wi`cf}6+5HFs
zRgu5IU`;qoN0{AxC==KC-IX+3h7&igt`7hI=ehm+`}^y&vR=K8|NqOG-*!vd*;#Ac
zq2^z8n8M(F>&urfcgyeFe*S!7qVnIre~%`8{P6IwB=j6=hE>hy%@0iacE<R8h+f&-
zTVHQwuP-evb#+-8A|)ktsze;(lmK2A&J888JBw0FN=%ga<n8w4-Q6Xs9R|9Y@sjG7
zuU}L5LQC8ewt-AKvHyP@x7QLsnpF7gSsSmkQc_ly)@*2Aau8j4wu)iG%&FSp|Ned7
zpP!R6XD8^o6VR~)f*)Tj?oZne9Y$Osw1_2Q*Q-~rZr!@2bse;3ZF~OxiOTL;_Va7M
z-Sq$acy0ZyT@UvE+Iy}2es*oj!Fov1{eb)X)s6PGr8)K2v=`h04>*d)*KCZ@JO5|S
zym@&U89V;TOIAO5dSo%gZ(`@G*Y^Jj|G(l>ME&dP<NJ452B?C@TX!hRSeL!|@bEBb
zU+dH7pnJ}8zC*MyggSsWCxe=d-A8$)&1CHV&2L&CqLpe5?XI1$289AMGjsXf((4Bf
z9C$PP?Afzx*RB-=T{zVP>oSEpO=0l9&MsdgAP8FiH|KY9a&qJuXymQx)nJ^q&Rwoj
zMd{eFV{Gi~+U4tFc7iUJTTpcO(xpi^Axl0O7HEN!gU`l@ocw%wS=qZgi=Tro=yX~b
zpfy#izpAS0)$7-r=YrQDF)%!tJ6R~f^Y-rY_uF!BgIX!!T`f)vL1ho9(y`@*`0s?R
zy6S<T>vut5?pJ6ATfk+h33VI8=s+F=!)PFY63b{x7)=SFbTBBh!?Y#O5@Blrz>cKP
zj2YOfK^Y)mmkn?r{P_Rou75}9*Z;>qo-aQZkfNu;E9kt%CgJ~#57E+rUNJ0I(kET(
z9hf$+ac#Z0XhYaLCdWwyqR*6E#M>kfc`%<0kZiK=T9laaIXll#!!fw2*Z6aEO#R`M
zWlzt%p114!$B4N<?-+l{`{-u*c@Aj2LI5!ds)U&CL=5hCAZA7d5p$r3NmS^@O3;WQ
z2dsz9Fls0`DKRiGAcg6XO}p3`UX}m)@^Z4ezn`BU=unF2?Rl~?GHa?qd;Ay}7#0Y3
z2rXFc>+gTR`u$#R|2{d}TNx(j=2%|7{MZ*d)_7pikA=!5TXS!3+p72N{e62In>)XX
zva+(izq`Bj9JHrX!W76<vc6kff1Yjiwkbt#Zfxw8HosSTJr*>=wJvtIm!^Ba+}ycy
z*H$|~-1*?X?-Yi+umAk~ymX4&;*C{bU;RA%_V)Jede5FcdnE~NaJ~>Z4LX$S?5fb!
zY;0_7?Cj!l|9(8~FFgzGzn4gWY;{{~Xk=7$_2FcHyGXB<t5$8_zP)cD#FmD0ie3!s
zKu7W4+M2z5*)lIp(8vex`}KBrDlh+E%YAhj_tj^!A!pSvG^m!}zR|wDckYU6t`{bY
zSQbnP>XWzscS?KxmUZpz?Z1Ej&i|pCJo`Q;Kg81)?8H}{?f$j&ls|u%xLx$F!UfJ9
zLJMxk?x`reyv+A(+U6-ip;s?1cAu^vANM@~;vffg0ab_a-4!1nG4tD8*kxDp;({~3
zZOD5Yn6#j(LwNN5x?0Vp7c*-9{d~T6@7_}Aj5tGy2*|r1KYm=fa%Iq|7cVmM^WVRD
z^X8Act*!0Xvk*rylyHC)_1pj3QCj%_@Avqx7Y`gbaN>l=ZFOi@(yA-Ijcwh+_3J@r
zP=Kv%sJG4cYIwIb8GI4f6<erD7dGyJw=YNYBsioP7#K#2Bxp#0#-gCi(GnRN5TGG5
zD08$zf(8U~g)~~pLjz^BJcS0tXn6`98)YD-JY6RtFJJ$tQ#~y$?Ufu<;|qaBEDyTU
z=T$y?b91w>uy848uPG>W!7YKUH*em2`SNAjW=j*3D>-J~ny07fuCCSq$ulr4aQ^X9
zxnyhY?{A=6U`qX=qe&%d`Cbj{zFi6S=jGuE0UcI+cQZr-!`x%Dg%+&-`u%%<PtTTh
z?QLzxjvadi8_Rm3vWVq@TB!WHJulzy4anUUkQ;X>5u#D-<+hFd`O?;b-$gEDx^N!w
zmAd`M?=Gk}z`($uGn4!3GuOXCPnWZ==AN&6*W!X=)hp<78wQ3KJ|uMl*8O-oJ^qy(
z)GIGEK!J1c;K5z)5Qi{Wb!jl}NP2&7@BO;pZx0-BxZ4bskR9z|Ktp1*I)?_tXmt+l
zGmLhxL3?PysRF5cJ=&T9hX@11Xln)<5JRCgbHM&qZnFF}Tc#4(YWCc;mka*Bk=-(L
zX|vhpk_&a;_;&CGGBNyFeAldO(L2EfuWwjq7<;q7O@F`OtPWV&x6S1TN)a0t26_Z8
zBJ4YZ!{2jrt;?UBn0WH!$$k6w=|*o`6TN-iJLsCv20sO_hCFi<6B9Etv;R-s<!eK{
zem<L>Usbio6FSlnA!1+BvhLfx>i1i_TIO1p_x1I?ySaIJ>Avdk@1mYV%Rs9}4aQx!
z_SOFO(%fD4w(8~5>A!yemal%Z@$c{N(}VitY@<L&vw;%Af&KF|7<XMibLPyID<PVx
zx3}efZkIO8$%r+XeYWZt%(zyNnj=SCUMKzh^z?JPt*tFHGqb5Lv{hF!CAE~F;mVp-
zt5*B{`t{38^M3vRzr3$n{P~|PW2!s8%klS3EHSfm%KyLfH*I<+`Cz{L6o$LmPftx<
zI%VU=joEco@9ykum#ca)OMdUqXR{wae*E-I`oV_{tgE=M@!mRn&vAC7RfatzJ%A#5
z{$A~+Q>H#WU$>3vt`BISZu|W@?WKQiz6M<)UpV(Y=q82~=z^C8&0ntAG6c*%JzYQk
z`Q=_|^KH4%8Ttjvok9z$*G6q!HRV!{S$cZ9_`Lrg4)Z^CI|&Vh1H2#wrc%Avj)HE?
zsQ+=<-`>~PcWThD-@jiYMgbSGyvWMFwnlR)=)&<!OTCv)nK^T&ZuGV(K~lZPPEFO$
zh8$nO!0==G@AG%=#N0Jz*lJn-@6WBR+5ao*-rU%Db#?gmdT5Js!Q_fyrjo7y{{4IS
z@Zrz;;|rgDczC!A8te%lk;c3}EQ5*?Zkumh^Nx{Wf%jay+FwsjPM$n@^7i|6zpt&0
zW@l$#zw6<hoyG2bGBfAQdGq@9?M<ntca^^W_m1~~X>Wi3{y(2i>+k#V=vc3`_`GPV
zh@#-&;P?Cg{}Yd|3G}M`aFBg|-LIF?c{@`d@_yr+c5S9{y4{Zl&DZW5$k%?kIQxdC
zH+%iXA3uIb=kF<eabe-Td-qn`Mg3WH)8)tF*E)0AbtV+dT^e1{&#rSI{>%KGJ9q9%
zXJ~k5WNTade((3V(i_mB-r}dKZqClm$;rv#cR!yszyIf}bVK&6nKO6(ez)5jl-h4?
zOg{c+=h9}?@>f?@ipNzvTy^;3#f$&`{kyhxZ}ZZx_y7MZ&&qnW>M&^WmmS-tk_|0#
zl}{$7ZML+rskt6ou4@wya?_4w*?J&1eNL~xboO%t>yPLh3x+MXx8>fB694`EJ!ma)
zt);2y*6Q!?uKlcjzxVpD)j!Y8wRU%R2Q8Vj`v2OTq5j#<=ku~-*F<bg`gh^d!AaKj
z|7uoyMC`Bod*Z~2jmgKOp7-1R%6NU~-oZ;%f4^Rjc0B|#z|HMm)$6s_ey(1>@7Asj
zB^xR>Z{B?Q^5wH>#TgkZZstr2nx^~vZgz(F{tPogv0V=@RUTmSe%<E1eFD4Af&0en
z3|lU)-~TTP)DZ$5s{ZF;BQtwcyRorxadC0>wKXeOy*kjy{Mnp`pMQH^Ms4L`>4xlI
zPbT}X`d|Lx!NH%!-t6A*J|34}zbfnZx3{0w?f?HN&d%10-)EElK4GQppAUyaLuD<C
zQvN;peBOTjs;slKOxNn`%Ie&^x3_xr{n;~SL{uDpBz{)ykkj>eEB=UtlXpM(@2ZsB
zz)(NkT$RBg8nmNo%B2jGu=R1a^*d^Rf4g?{&Ye4F&iK6Ee&<>B|DaQCywXwPXVX5n
znWy|`e-Qur(o*lWvDP&|J{$r$Ew8NX+nbx4mrmJ!zpnb<hKm`XMHWXoh4=2=`;Bi}
zOniL2=F@rA?`GQX{Cv*ZTXSvbrjiXR^QvC0Y;0`2#_#v<$4ju|Lf7ptk=t-!((ZrY
zsNKO?|LdM8gF@ofRiVtx%+1Zs-2C-_J|6%1y0f$M+R;x>Pe13cuBu{VV~eUc-T#pH
zfa%sasrpApI-~RVmOccfwClf$*REO>owxJpkB^Vz>wczcKHc~G-R&)z!LQ>vx0Y{s
z0=jEy<Hn88=hy$cvNCw@-n~nwq|dM2Hsw;)f}EELo?oDS)c-%;C+gQ-%u$MY`pf${
zBe%}Z*ZOz=Z!Nc1ES#Io&~7g)d^An0sXHqp+(<~!vx6mSW5}#B!9w9lZQHiQinP4p
zNt?Rn?IMYS$zojrSqlV&yu1`<U%1tKp{>nBQ*+t{R}NvX+|33aZ&&Q>E`D!dUVldZ
zxoq+$o7KkYeV=w#6#n>A={J*C>;Y$o5W|9HpMKW-JY0LT`u*PNL0TeQH~ETw-Tzj&
z^v(YJd3SbveAKPKuGX=-y86JO#NyerW@(9agYN!)@d9++=w5LNs};NM*Zp4mY5%n#
zFVJ+#ym@;4Jv*DVeto&@zkB!Y?xRWJ;o%DdG~O%Bx2@h*)4Q-)D&qWtnJyJNW)TU_
zHxd}<^&B^6IN%l)74_=v+p|#*^!NQxT6*QqojE~O0i8D1*1s=sFS^Dx``VvBe{S4}
z@OpWqQ&@9pTibu`VqV!6tAYmyUS3{){P^+B)0vr>Pu_8CxO(;Kr^2Ui{(QY2fAWH(
z@0$7xK1yradBs*RO6wdD7FK0A&}Lcxug2fs|C`3vZgKr>+qNy6^6&5W`%^A8?|QIo
z*|PRHfvqt&Kvxk3d5P=Cop~o+FYe4<=;rQz{Ma$jvT_rXFPHu8Q`OVc)6dF=#QOXD
zml@dE*%=!fTUyTCD<dHh5fful62Sag-efX|{~0#!D{Z`D4>~Jrm>F^!jvhUF@ZiBU
z_w_}%($CM+wfD-ctgP(Szq7miy|ey`M^+Cd1NP;XIy)Z*Z2`%<TlMnN(kYjwYKMbv
zn5(hNf1z*=)L5`P_x$rs-KYAoyQbJSeQ`V^YqRc?az%Pt+9vzMU%r2j-&vG;eO>I$
z-O=0gVjrL8l?^%i%RY3LO9hiz#DcV(4Ggapy%-vzqGMvtm~{kuE%%#yYje8))7`qd
zx|&OqHtyK9D@s!K{+>$XrCavw`E!B$Q9aWQpS(LeCd#Im_v-I1czCFlTO4%t+24y-
zE;x!#eZFYjx^<sGVOqA>Z~vPGyG}Vz6y==x?Afz5EWh6Eet)m>`P`Hj3O9b$@1E&W
z5n~pyU>e*J)o-_65Ayn1^Y6HP{h6ml<$GtJf4;f-e*OQwix(?f?<;*BRwjJt+O=zM
zUS?;<2>cg5koM`@?fm~gpU>Cd^TFxcrky)~&d%SbDc0TU^l)l;T;<oR;k>dTcfY;8
zy}5b+|9`uS+&0y($;{3+E)S5t`R!ewtTm79l>@)n&#D}Hyy+m*)|6(}i2mO;UsxGB
zir>yUo3nv|f$N5+i(|;Gt#69IojK#9H$8f7z`nHa@9t_YJ$Ufo(xpqMTsn2?)W5&K
zgS~#%*!_O<?*ezjrf=WAZQQsq`M<G|QPHzAGpAfS*v#H9Yd!6+dhf#KSErs9wK^#-
zy>k8f_8mKB%$+l1#tA-9@t^W;!KzczJKXzZH0xJ^MD5DT%FETA*-bb8b2RS!@FHh}
zL+{N5#%qh;ottYdy_b!lq7O6!2fFU}ee~P6Z!cfIthw~(&!7AD?W@+eD1UdSPu5!Q
zDN_RP_4V=fpU;}RySulyw{JFX<6D)Nn7DBJc5~}<c?$Do<?U)JQc_fATCdo4iht%*
z=kf;!8aHm-nDS!5EjBi`uA`5>d?~4}uKxK198k@7tV_%y6t2|PU%q>kjr+rkJDQ9P
z^Tcj%%RS8(6&<Y~w<luX<H%=UUS59M|BQbxLxi)lv$NOAwQJYb9J}Cn_37uIGiT1s
z`M&DU@rxfneq7_>@r-X(Zbik9CnqQ0EUbKS?AWnQN{W(JE4IzGF27e@^tU-p=%M2~
z>k|H?X4Y%(1KZBQ?TlL*v@%5NY0=JO$BwD*W~k5sdF*d;>bnOUC9Eo>X1G+`lCWA~
z=*)fxbO;Nm;mg3#^5F+`T_}SBqI=ANnDKW)bR7f`<9Z}0bvf*&53;Mm@9wVB?6=%u
zVPR2GQPZYR7Z(%z^s^>8IXPDh<b)9KyK}9}&8k2v=nfR8r=~7lz54aS>7rb%=JzTd
zZ$57assT@&@Zg20o_HL5_o?)Pb?e>*>g}%mT~=IN{3Qvr$Us33y#6z%;n%NUJ9q9p
zWz5XXd^D-Bq@<*xV#lXHTeGhRdHw$We*Mc|Ak|rr)eHWuxwp41oic0Itan!T!JC2#
z@9S$$^;#OFIsLS$nc2G9PkP?o-rCyQy}i9VW7bJ!gBAq7JO5Z3WXSScSFVJlr$4`P
zCFI+u6BCu^&Yk=D^JigU;k$S5?k;=#>gCIqZ{EzAJ6CtAm$S3;RCWE>T_xAcXU~}N
z;?0{kFJIo=njIb<etq}*eY>}83Hdx9v<PyZ_&d<aR}497n>TOPT)J=HzTM~Ktx8_p
z>jw2K<m-McT)I@X`sVH1yEkkIs66=e^z`bgs=A*~r|U*<JF*6}>~vl>%!Hf4azS3Z
z%isU|`F#HKXV28s)la{dvnY7*>gwv^;^O6gb5EU{V_B?b{o}{u{>@5>H8p>($JfuD
zJ=?nSQ_8n%?ecXIyTND2%u@#^frQ{+uU4;@l9mqk`uq3q&l<a>AsQmbjvXu78RM0C
zexB{yw{JC<E?c&2O3=!XU!P9v2YWq!{Mgyq`RmuOnoIXqe@}aR=;`U{wnso?zDgA}
zpv3_je#+R_?fH}%9ew-Ip+g6o*|q(ry<a&YXjkRuw6YJEyg@s$KK)#i=Arrb<8gWY
z{eOzC6frU|Jg5YnIK;p(ua=kh|6OLwk{1_h^?&^P{oek^gJ!99pp7N-xWR6IV0r1@
zy|_yD)-!K!Z=XGTwt0#z0|Uc{9MF+J3=HS@J$?H0&4b(d``4~tzrJR5h}PBX*Pj>d
z?2P~gq|M`3@^zp%IiR+^_V>3p4_>d|zs_d=jvY7d-o1PED(Ge`m_qRSm;;;jr=Na%
zd%OP9l0QE_g0>i~TeobAn3xzFJNsttsHiBjqFRvUJ)onB7#L<Bu>b#O^Rd(JA5W^!
zPhofOm%F<$**!fyJwHExV}y^Uub*F7D(DuW1JKiH)h=AS7WT8~`MJ6Mzhyx;1$}*W
z^-UdT>d%^gcXk&4{rmTIEGW7P>`Tkv-MRU~6qNn0oO|^ut88&=jNbJdH+Eb+4LYYc
zeJ1EIE|4q($HR}Hj0a*nAPo^9hS*ym*#x9O5$s72YgFlI2#kinXb2305ct7g*<2rO
S+Oi9LK#QlVpUXO@geCwZGm0Dl

literal 0
HcmV?d00001

diff --git a/SRI_matlab/main.m b/SRI_matlab/main.m
index fb9df2e..ef695d1 100644
--- a/SRI_matlab/main.m
+++ b/SRI_matlab/main.m
@@ -18,7 +18,7 @@ coeffs.QuarantineRate = 0;
 VaccinationVector = build_vacc_rate(0, 0, 0, N);
 % VaccinationVector = build_vacc_rate(0, 0, 50, N);
 
-m = model(S, I, Q, R, 1, N);
+m = model(S, I, Q, R, 0, N);
 
 t = 1:N;
 
diff --git a/SRI_matlab/step.m b/SRI_matlab/step.m
index c874baf..2353e56 100644
--- a/SRI_matlab/step.m
+++ b/SRI_matlab/step.m
@@ -12,9 +12,9 @@ for i=6:length(X)
 end
 
 if length(X) > 4
-    infection = S * (1 - 25 * options.InfectionRate / N) ^ I;
-    % infecter = sum(X(5:length(X)) + I);
-    % infection = options.InfectionRate * S * infecter/N;
+    % infection = S * (1 - 25 * options.InfectionRate / N) ^ I;
+    infecter = sum(X(5:length(X)) + I);
+    infection = options.InfectionRate * S * infecter/N;
     new_waiting(1) = infection;
     dS = -infection;
     dI = X(length(X));
-- 
GitLab