From 9f20dee3ab9dff7c948770c39ee29bededa659e9 Mon Sep 17 00:00:00 2001
From: David Lenfesty <lenfesty@ualberta.ca>
Date: Mon, 4 Nov 2019 16:44:56 -0700
Subject: [PATCH] Added plot of diode data

---
 lab_3/data/diode.dat   |  11 ++
 lab_3/data/diode.gplot |  13 ++
 lab_3/data/diode.png   | Bin 0 -> 19838 bytes
 lab_3/data/fit.log     | 288 +++++++++++++++++++++++++++++++++++++++++
 4 files changed, 312 insertions(+)
 create mode 100644 lab_3/data/diode.dat
 create mode 100644 lab_3/data/diode.gplot
 create mode 100644 lab_3/data/diode.png
 create mode 100644 lab_3/data/fit.log

diff --git a/lab_3/data/diode.dat b/lab_3/data/diode.dat
new file mode 100644
index 0000000..f5bb925
--- /dev/null
+++ b/lab_3/data/diode.dat
@@ -0,0 +1,11 @@
+# Temperature (C) Diode Value
+68  364
+67  367
+58  382
+52  390
+49  396
+43  405
+42  408
+37  415
+33  422
+25  438
diff --git a/lab_3/data/diode.gplot b/lab_3/data/diode.gplot
new file mode 100644
index 0000000..6e49d3d
--- /dev/null
+++ b/lab_3/data/diode.gplot
@@ -0,0 +1,13 @@
+set key off
+set xlabel 'Temperature (°C)'
+set ylabel 'ADC Value'
+
+f(x) = a*x +b
+
+fit f(x) 'diode.dat' using 1:2 via a, b
+
+set term pngcairo
+set output 'diode.png'
+
+p f(x) lt rgb "black" dt 2, \
+    'diode.dat' using 1:2 lc rgb "black"
diff --git a/lab_3/data/diode.png b/lab_3/data/diode.png
new file mode 100644
index 0000000000000000000000000000000000000000..6e77109c67e4bf932679c6a4b5514652e639efc9
GIT binary patch
literal 19838
zcmb7s2RxR2|MrPUWK^QE(j-E|C^MxcS!F~BrDSh1uF9w+MI|H@Ar(?2vn0tLWo2dW
zt?c(WbwBrb{?GG&-`D;5+<gkyb)M(%_dCAF_c)I4bzSxNAtr`(3<N<i9XWhJogk<w
z1VPn8PmO<hm!gz{Kj@5<4jmvC$p0d<VuJ``J#pl~z7x(5`x{;4PfYv}A0GO=SMtkK
z{>#^n9$mlf@n-Q?0$zGFE7&%(es?~7G|M_)J)wf}l>Ca0mo`UsZA$r--tvXm5Jk<m
zgiq%h<N7w*gX}Bl9{dcb3YX%rkvHjl7~&%9X>-3-dd$K_Qc6jPeHAW3_4z*grR1N7
z_zm&Lx=0EYL7b;h(-A}@6~h{WpwRvQf8&^tp`4r?yWi?{>*{~TP+wjzDOvmZ^R>KI
zXIIzz_wEU--%1ejIXhlbR3c=<Lqg(TylDIxv*FI?yu6JvGBU!#!aH^(KYjWsKR+fl
zRg=#JS2TT9hxYsT@0pg3?^NXGPADlQ3@PQz3}!rj{MfkQ)^>wzw~EJ#re%*-Gk7N_
zC*L^AL=eLt9-Pc=NVV7{Az_egl+U<b#-aCaNJvORf}q$1OUuE|&ikRE+nS%_u0{3*
z4H66YA|i4;7v}qF5>vFFbDYAL*~N>3G6kijgCD#ZU%Yr>F0Om)I(adE0<_=j>ejNe
zr#cSQ7Y1$3>8_lfo)*w}dbn)I?%iyD`3h00j#L*fTyU_r&or^c)3q>4*5=8%FF?!1
z%PYc4pK8_I6cA&S{I(#xQy)L{`GR)G<$NOnwNsDHFW(c~u9lkn<=%aMUBix}N_V$g
zFHFaXtWnDq6BCo>5#i&DN=eb+XH(zR(9qzmz~`IiGTwdXW=v>k=)r>rnV6U=bCRN>
z@M>$<tl^&y!|jbTZJb}8?l=(o<cV>BUun3MsiEO|<M`Ovt}C>RqobpP%{d|=LmeFx
zGebEGBL%##(h}@Fz5mO_B+XXtG-9d!^l8h+jr(0ozo*-Uhle+NF3cSX+H%3n?2GnM
z`jN@W`JUK_SFc_@3=dz%ck<*(H(AfdcMdmZ$@3uM73Jw;Ey?BSxey&4oh`Kb`T=HU
zW(NlcH8nLZ=B<18zQb2s+~?}*6O)?5S`YLDw0=K^B-@zHU)Z(ND_5Rn6v}yB{Ep?V
z=j+$6g-^e0|M6ob4-e0}b$g9IUMnsx)^Fgz#y$9D-1M{jb$|bGS@(B-evA`N%ID6V
z)7Ni$svOo``jxR`I_%EB9`V1u_Bv@E6`q{n!jDEHSbA<Qu3KDtc3zt7SM+1g`q^~J
z_|wfODXk$rJ-zN*!$x-n|Nh(0^)6gjr)|o1&v<x4J!F3QT6wKR{fV)$RrE`?Y}ry+
zSZHZ!X}=XPu;s&^7OGDjwsO<m8x;<lm~5)YMYE|W|D>d(q|D9DB`;l05U<X#IImY5
z2z&VO82uiX(e}X2O26D>8>4W3vf_-X3^E+IZ{F-(-`8QAE+8O4OG|6-i<h2w{M0v6
z;ERKRG(qt5Fz=%wh&VPjQ9OJm)m;34w>%a8sa(F1aSg{AZ0POPS(L8Qa9s_J{Z3dS
zRjFS6+?Ow3Vp-zj<L#H?Cng?$llY+S<kZpGNw;+AQd-*hRzLD0#HSjU-lS$<rGVY{
z?Af#DFJAN<UqTSYSE%nZN-tw12$Q^^B)sDZhCi^RrOV99s`<@DEgY@JeRu`hC`v>^
zR#sMY^vCz_qmq(t4+;T#%&MoiMQ1-eC1>;FW5!T+#<`ME7P(o~{rmUtc#oy!vtGe&
z!7i-l-1f=;;lqccqqgsWsFdDZo}_6HcY^3%JxSu&U++4<>iPl6D?o+0tJH%zB000R
zO_?rl-n?17xeoH*m^j6cMaa5;xhr5=`Y=igm(IWTP4;qU`Ii<MwjEz8W0W|ypJH3J
zYS;RBywjdx@lEM#;?7u{KmWO;WMXvGb;^L2Aa0WP{e5xaT4fD!-}kQkfu!2EvKQuO
zVpSq4!?8dQ=xK#SaBk%0f3@``=IE=dD>$(egdrDOQs(*d=eguOX7gC&?46tz#;a5^
z-o5K4AL{ZCAE_tD#*CA-qoSe$_uy-+vjpf?^Ytr`VPbN!w5TX({+j%HdS+(vojZS6
zOFnw^=<M0EH<;(R7&(NCOH1c*xXT{NC@sDhtMWS<d#CB~p5o$T78VvgA@}axGZ(*<
z=fi}3%EBuv?L53bxTvJ0Bz&*KHogO!R=PMl(=WO5_N^qd?<#KYfgeBA_4K5TKGKI~
ze_Tj<_3BGuVMTd)q!c|th+cnigoGpibq?L-QSYxV_4V?SwtDmB%Ta~TpZq!F3!|eq
zDbPoh#cmmD&Pjanf_ue^ioVVVlM_?~<wAQ6Tk+l{Sr!fAc6MW{1vNkGTeN?^gO>v!
zSGUJ&#NvSKojx6Z#Lmu6LLv<?F?9})-g(*by-ULX4&gjn=PfJ-y1V0EzPxqo*0-uE
z-E6l#VmYs0SAYAa>8dOxm4T=7IyLoGoq2i}>YcOAu{xT;{{H2RjJg?)1MP)D5oxzF
zGG@x7_5=JyPQSCCpLMx-@uHcByL%RvU_8{v#|OC1&o=q)w5E!)bW@Rgm~6Ur>xbs~
zNr}tykrWHukGV8YV!^@N&FmrSxg#SZ)RzrsCI@zYt4p&o(A4}^@#N6<?c1|5GdFQ4
z_~rlb2v7?s=vb|hn}S*<>pqQZ7!lu35cMO@nO)&D_g5X{y7DqUzIDWT*|KHWKY@+k
zCMO-ER6n;`?-pBNxL|5}Xs)HJOTuIJ*LK}hy?6H4S&d9gs>{pQ&oOeYTc<5_|MKO_
z)Hxp+tE#KpN1SoSm4bz{+@@^x^d7c{Q4^Fe#vgfx8#A1APo2_ysdwQQFRl)lg&PeP
zG4OD77Z0;AHny9Yv@omu_AT}G>l3ctfq{X^#s%B*ml>Lxc7;h;VVNGgQxVk+8yUzd
z(dOCH+nYV{O?U0uwTU_Nd-<YKy<F$&>+5G{-2?>%%_<ujX0SzcbabjsPM%yHQTB6W
zL~S^&TVntIJbzA!#q+{;BTD3BhqdJV{CqW2wihc43kxHoA1Y2yUte>&-B$30)YQoz
z`F^FPrCzLe4!5_osA_3(tX{1sbM49B%L#cXZ!3E!M{)VexVyUta7%#%)PMhOZfPld
z@!NA2*<Xe~K3=a@-@g3?o{YCI`IVNIW|#hbD0}8*%_KHj<egJP)H~f+Bt=BZN=uJx
zX$23iCx~7p59RmQ>~hPW?Du7EweTT`XAHE(dP`Y&mb0)p+1WK@U6Kg;SxrSm@~&}T
zuGzH4gFbY$fY&o~tW?U6efy1T*CK&Al=?dc8vOfc3Io*HM9DJ;h_aF7dUq87&drvj
z$FD>J_i3OXFzD#MMHDA!pFbb`p1z8em6EoWn&1~;OM0~v2Lt6dEltjIUb<0-AX>cR
zCr%#TX1zian+{m!?&c;tH?9%mJhO_R*uGcPF5iXKd;E9<6O;d~Tg)Vy<}#{IpUw@A
zqP~<1YN}OBgCn9U@GmlW^4ry=ET?pJdyWHVW%g4gJkRrGX(R9ZrZ=^wjx^5@@4_$N
zu&{<lGMB7htR;x5xQVm5z1FLDPl)HO-F0r$rcLGHw1oUri^XBvC}!GY^ee6(0F1G-
zvy+&tD$8$}t{MdJq<G}H;{u3nw6v&)@-%Eop{@WezzDC5>j5WLf)K+SYOGtg4i}V~
zmX>4PN=@Ez^_^1&-d=tHrdJ9I3U2Q1a#QVqE$(-1z$a$m57g<_j*h}YLT=;T6}+Ca
z`$$YgmCK+k-e_GC3F49J|0<Vx21drlyLRo`9K<24SNrB7hr12R3m3zh`5)cgYgk!x
zyT+FB)z#GjQD&$65|xycdP-@D$WYbXVGiNoOQY?Do0ab4o`22E#6Eor7W?(<SBIgd
zthXna2*ccXJes!*(>h=nIyoNh$BrF4aA3()S;KOiJ`HvC?pt)Y@l2*8!m7EACZ?t{
zGc#6=@5J`*-MeQ`T}g>Q>xPxrJ*GN>C7PVC4%#o4_aa&b?R8(XPYmg2xqdAxQ>jGN
zG-Lcpb@lt{-;}{jPdGR%Fi>4pRo;mobt>l1JFmXTQd?FQkZ<Je<8!pr#Kc5MH?{w(
zmciGmDozQ@@9!M?mQY_FKeU8;pS5+;ojWUh^G3$Ug#-mP?9We-uw}Wxh9@_)Mk^Q9
zuC%xqXnX%ry8Fx|xO$9I@L)q)S#z_TfWU=MH`gjFD_5)IxVE<5dGv^fgCqFC18ohB
zHk1}!U8O^A4t(r>2?+@#1XAW#?K+~NQ5C25gqeHqZs*~a+qYQ@XHofkzTV=NI!Hf-
z$Nt|7%e~`Z{BLdAK79B9eCjTxCf0L|USv2oKxNgCTJrt7TmFsJjhz9!a{8{7`T5?g
z>L*VsDnHum@CC)@>(_gsyCo!Wh?SI;?_JV+BLV)GnVAVt++|w2xQh$VS8})HjE;@{
z9BTHQn|5IL+kC0y;X^KkYvJMCVmiM&iL5K(;f>ym95^*~=BW=Wsq@Rrhe5exrDzGt
z=$S7()~2SW+S;}CZ!T_Cc_ehH1U-dvn0Q@7L(be}L*O2~UuEB#b?a{A8%3$s_Brav
zpbpH<@!)CSb(Yo-@L!*^edk#gNm*I_^{nR2*{3JXP`4P_h0b5NknT41%EyOJOhbtM
z=cSt;-!d&D2zvaw%%FQ=;?=eP>@8k3)HgK=voY-9y7jcvfUlvx9>>S3DKq5Goik_8
zZk~F3=Z@O`{a43D%+H@UOxC8Oqhpw17#6y7I65XqktC2lMj!dQGn|G{AfD@`R}RbJ
znFj_1p*ImR$UgRB<I3x^b8~jKwqz?}WhF9aqpRCsX;mU9Ah1@-QRp~P)&1q3VAWIQ
zyoTje#3SDowC2XfK_fZ#4i3$(eJ^g_yjf%58O~Z+QL%!R_4cDj9-YN&4;?D7<o!80
z=`q~0Pg9d!Oy}gu{H`&;9NV0j{EeU>jROa~5)v*dC{W)TqamdIqRNcgttG2zB;pG1
z32J%FP1Kw@b7pXGFwo)o^XE^Gg(#4C^uP2a{6}4$G@X+~zk!}fpDY5BXl;%SH2GA&
zplm(<@K^u)Vu!wO15H_B)9I<H;?mO6s3}=l;YA&v3JMx(YscTMW=o>JcHnW)Ruy6W
z%m-u<nTv9yeZGc-*xSlk`k#P$k)!`YAQ~BAgY4PvijD6yK;7wRmyf@799Y+Mj-Kc}
z_99W?gQE8Mmy!}DW@Z5HB!_|9<=CeLv8#21km^3L9~T#wiHV5_x6x%JHnJ(daH2|H
zTFg2yFn~gA9}C7zHAr)xk?avZ`|L0OJKl<lR%BiJ`Bu+TM`ta`*xLcb$4t(gS;5ZU
z+1~DjnnQV}-j{G2MJaD?er9OjzI{)hK3$lu$-zQ(-=ZOU={e3^-n&lD!|nL-(&suU
z!@_3>RT;A^#%nT$pKkowo{6D14+#!VyHK^|)P5?;-VZ#YaU57y3Wbo3P@yj@EoG*s
zx3;#9lwu=t<<<L6%lA|5s{5LtH8(q(YF_usuB)`Y-9$n{;@v9pw}X<}scDO?{|>1K
z(iA;sM|>||exGzE@BA`4qL;=ua^C<Db*MIZqpYm#mMuPe*|!hOPK`}YP6i8~e!Q8J
z$Su<N{Z^Wej;T30at`ewd)epu)w*$a<>lp5V_gLwKX!L_m#<z;oZy_K?90;ocZ4ms
z;MUTr0rd3t_MR=eUvXyz(Gn9kA-sN)auOZQ)u<@GMgc-z42zXcc4icF$t7r2{;ED#
z8OV|cdP5?i;Z8Quxb>$`pB51j$t+<d9{jQV&y+{S^H79A-5+8zM&8AO<T7HDo;Nl9
zIpnzzMmAQB@r!Ruh0d*?WDkF}g^{;H!@_Kg))3v^)N~H^tslIhXXs@*-@(R=%tkx5
zNzW^$)Yi`K>eZ`2d2U|b#MszFPF#c`zn&MZXf78w_t@xYZ(pB`wDhpk8iM~3Xp#}j
z8hSqRncus&8n;M^`e$)Dp{12*kkeFAai`=QG?y=5zF6|EpkMO-{d*`GT+9VkRjXL2
ze6sh?P!rq#Sx<D-)b0cZ8kw8-cX$8lKXhQj9bex=uCA_V$&HMRl$HJc{FE##UOe7_
zd(0F)x7r%L1S%KFc~Db(rU2aNRb}itZ>FXFOo^Bt%ow_H;|2;zs@v2ky8LjNOPc&_
zpFVxswoQzrLdKf^C=@K~HVJn9r&h!zW{Fsm`;6nm&jgWfq)OG*)rA(vvpL&ch`rBS
zR{Jka=_ShHn>P;fG*{Mv4lO4SvTgod>VNbf2IDC>wIuG_EoWlNH%eCJ@1IC{{%F`a
zlcR*#cCM!~23=smE$*?A5r{x58f^&q!9Pnu45h8Djh6cT<UqsBh)%Nh_n{pte=eb2
z-Y#dg>-C#A9b;4Wrc^|F8>)c#mD$udfgFQIr@w1Vf3E}QMELw~-UQ2v?%KV3capXc
z4GoRUX1a;Un3y*gYlHzx6f>j$rFR&dXiH5^jfsh=v^#qA=%-uUMdwmF*ue{1VnqL&
zoHQ|QV`ViI;V~^BC|Fci7eZELfO^df-K~jl0Z}6Uzfl}r!I2|Js-v}qWL-u#@wNBY
zrQq+EehzZWy1k8$?+LTYKJtHs^X@dPcGP}fvfpSw4l-X(*Q=Qe!q~`DvX6eMp?TgB
zWu~_acsV>g4BZKhmvAsAFRzucaaEtDYHqpQ)YMdVb~fG`H0UYV{O92M)k#K?l&I#)
zswy1xU?H8At5?UM-{s_#a-X*6;NUR3@TVkozJ#8jktbn4*9|X_?K0L0UfSB-U0YfC
zn`Fv2Cfx;}U}bG-X_0q2NEq&Z(?t7`M>I1NExcL5&)%;I9v&XL0K;@UQ(N0tVPTw&
z6@(#Cp0iX>RkfwBPs;7r@V={bkmCGUUtOqriWaAE<0{I^J6w*Pw&}d?<7&>%&cBr?
z*RQeEXr><Ru{)=zsHm;|fMgF4rtB8x=SGX^B9FaFNpb40UBS!CyK2={UsEdi#OV{c
zsD;_V!CEM#_V%x%DY437nyRYn%gl+R8>OTgK>POWahaQ*0KL(BWwLDV8p5PcgNlQV
z&CA={;nLQ~>(`Iq=(?#>@rPWjc^MNOt$OHDy51lTE||h0@KnRsvNFq*4jqbxF4a&x
z^K+=}Oc8O{(uqq)Tl>JlgW1!)3FLb#QWHLXRXj;4DPyJ)?%ygaH-uE7U3D7Lm-Z$i
z521GlRbR@?H9Rr#A|c_46f+^J7V+{WglF7uCyK-D>}>YzXpz{(Zzku?HKBDA5wV*6
zHNtRMU@NFo`?IX9Ip{N7yU(-g$1jP@?Dkx#3H1_O`QSmH!<#D0%dKr}t_K99rKAKM
z+Tu<2X%zEu4$&CsTv1W`3=9mkw6u6Guq`KBnpZnwC?}Vb@|a1GfyQ{xwr<oEIV&|Y
zv)){9Dpd|_!O+l~&YTiNd5p>eBP)Y<WstwWf0|`uWo0E=ykM1W#L3C-iYFnvOrhlR
z{#-%CrO2+f#>Pt96Tl;L>h<{w?2M(dD=6CnJq%ihTk~q3>$G)saj>%=%qqE!jzRva
zNA@l2U&E~iZX@|}3sYs6m1!cMy~|1*u2PXxK6><V)45+Q?`hsuQ4vhjo|Zc6S#R>l
zyaRd}(g>_XE2Y>w;_eNZ0K5O0jV=}CsL}V6<+yI$G)rYArJ0EugQe0zA0jFH-cRCe
zddy8fc<|ubHF|WJ@)6~3tH)nosEXbywRi7eV@6GF?Jqpq&)Akn12Ypv_wEH3<r!K|
zr)uZE%@FEK&k>Y9Ik_Crf*23x&6EvCnzXt*JH;J(Pr!31R69?os+hSiXGdi}6P<Wv
zix*a~v!lb|l5gNg=nv0gV!jTr7!Z-PZKYHRbZnxcRS*`FpyS8I@uQql9LL#2ee-%I
z#-M!BgdU1g&CS*B=;(lMg4WcRg;!i${HYW(W%_X}<;mm6u5+`g77cW}0~Detmka9~
z8qzEp9Q(e#(1<-=5l$=b-pcw~FOqV33f~!t6C&0{Q97s>7rnxbH;KHSQ~LVz<2@QZ
zA)rk(>%hw8Cnnz1dMmVMx{TS{+CDlRohfzacS{Ok%-wx$9^MZ;f`O3{z5x!Q)%Uj}
z0iLsYA3sJH9Hk?sdb+z&?37egVxh!k-pyv8sZBN}bw>0{r%yLRpivK?m$dILNBKG)
zA-lcn_UIKELsj*jo*tYFrTaqfLAu-h=Ve;(E<21ru_{P%9XB@KC}ubEBVX&K-kn22
z#O*U!M+vK@*|9P?)8YpsXkx0kmlD%Y5^2rN&4t()mhn9~6i`}H5_7I8!^x^V(j!}K
z4WaLooIC=zq~_&mzXdv~yC)11uaNknum315&B@TPWht+ls)mNKp`mqmIiHi0Q+g*2
zu@*Aa5vL_YcWgFK{XBLO$|<ah%`2~O*tn7Gi-09>GctIkoy$LdJkTW2v7`&iue^9Y
z&&L(0Rx4Nj`0?X|QS!U(EB*4zN>2zcl577|Y97d&VBSJPLRYR{1=V|IBtjgdx_Z<V
z#~dVT>#+y&&=yY0H=Ha7p+ZN_CFR%;6v8(|?K+F+W`0WS-u-za1F`(tl`B_tQqG~)
zkQ{!!k_fG|$>-bPd1p;bPH1W#Jah;@3EKOqy4s#JsY)MWiEba-@-K&Mn~cVK)+I}q
za;{wqsYgYn42ngz2(79hi~>CGI`mPnX+W%cz@=sS`k0nuBLd`66f<y?&?)WOxwFY*
z#$wm56jY|mO!kYKZn0$<JL~lWutl)4Hh`n+*HPO6agdbUbr>Rjz<;2iZWj`Qm`X~_
z<g2gvsido`>+CG)utIaeV}HQEi@1k7w{G6NLID<pn9IVvJLJ`C*REv?-8*BLM<c3g
zVX?h&1#vfQHuQ(LRV(OeGxqAQkv6%7*)hKTKF_t2jfjPHwW$_i>FHA`Ib+rF8UQJX
zkk!%Y>FHr%x@Kl(2KTAsc6DSavC4N$Ozh)YnE#mr1qAksEu?<bt`eXetO_6vwCc#m
z?Th?sVS1vktf}$DU+n5H$8NIn@(mcAq7>IU8q13prCdg%0|Lqf6EjD_`NSn8JVy(+
zqHBk2q<}3*eAE4V6F%DyW6u?(yNB$)jglX@?tR)H^D5c(L)(acZxKYeS7Z}S(x=)Y
zp5f1THhL>WB`4cxYil1q{CT{ms^c#Dir?*cRzZwW2+Q2<tP|op25IOffP=2Qyb*ON
zrd8ZaDEmm-cFAD0uSK7swy;oO(<T!*6t?aiI?BpS0R>8j4|n(X-|bI@!o}#V!0h`Z
zKK=rf-GR<F6lQk6vxbK9^mnq@h<76s6HLpN;s5^EuUAx5L;<TpdplUi7v`E5wreNJ
zH;HheowHO@D#FfTT6hu@vl%K9l!(!eA}>}n;+2)^Jpp&`-n|KJbT6IoB>}>|t+)4-
zo}LgZJ)AR5(yQPi#PY9kgw)m6n#_s-RnnSV1vhSN!|~ke<$URq9IOLDLBCTu&z{{V
z>zeuXOEu_KFm-Oq^5x4R12kV6lF2;2h9bP&ea%j&M?ZhYS`Kgax^`g0(xtwOVlDro
z!vLH=ed2_PnOS|G`r2Wk^XQ8GQ%(2n^16Tj8-yfK{SSG0ServP*W3vS5$MsXQ-T(`
z=I~-~OBU9O^rdc{#UVJKrd=Qz55eh~?+}!8KL;D1NKsRI`(h}bzgm4!Zei8RWE{pb
zDJbZ|t{g>gIL<F97$$B}4@%h0Oda{dwJzn6l#{T-8<PMf4iWv?_QF*i1BP5z=(^R=
zTpJojdVA-$=J}#B<DANGt!LiX{MII~q-4NF3GGsNSeR5f`QEFx9D`x>W>K6i;Pj?y
zj}E4`G&jFHA}!$VK35RH%R)sLCIZ2AF$-Y7iHM9yO?5<{lkdj?0uJRWGSpR|<xHN>
zNli^zs|#n(c7U(YNv<G7ecrqg!}&X*2?T<cz46d3dwYB7jJPLH+ORQG&XokJpuHzf
zep8}?s=YATVCBalGU&(?m!0haY==|S+S)pL6Nd$Fk@L1MAu+rDt<4jt2SGtq4|W*S
z^Ep9hw(U3yh)_`Y09!8M_Cwpt36WIOtp*Fz!QS4q=wi$iJd`V9o5MgpSFOtH8sp^T
zT*e1Y+`-Py@8(Th(cRySXna3<JRA}b7qU8mcd0^O-Gv?!5B-3ZIECtv9hH)5jdWjM
zU(v=JKcXm9(}(!kre|hM&z>DXPX)yeEDqcs#AD~+sU@oF5cEo9U;cz?0{#{nx(>&}
z%xtek{q9T)4tWdM;d(HDcJ8!=#f9Q7<urI;hB`B@iSbuYRV?bFEJVesSe4DY>dYrW
zn4AXcW1c+OSmsRR>h0LE0|iZgQJY;}&(Wed@yc{NyC11WfkOk6x(3pkd1PGe(dXlI
z9S`4&J_RsEst@dbF-dP}7=0tD^ywVhT|f?SEK$)<ZZ3GrS~NcQvw4*QGmrGI{>S*u
zl?ulnNej$Ki0UL7@Dh)@a7N+A0l|%uC8eZvM|V+fz#$k{?<ZwFbMs7CHvvkG-@ikQ
z@W58SypSQ!2kE-p3MmNOi<6B_)RvgUE0$(v>vo#Zhqn@R%oZZXzdC{qI!nTSX2w|M
z%sB=mGSCtCNZsYg%t{3utGG#bRKbghSw0YV#-9aKv(qnma?szkNu80A@v^)qHnRM6
zI4`lWu>j8Inwm5?Y>*s^zba9srKWy)TdH&DP!}pa4fisXgYkytMMXt86}hTe%|7ET
z@2`Rv%2-v%cTu6hW|fo|6!^M&Y$v1>g6)dj)AJ9zBpO4`!WrPSQ278#pV*#!5*t@6
z_D_xad<N4Ix@7H#ke^}tTb}`I0!NF9jjal&rg+Q=v>=e+j_k+%`}a{rtDn3=Xd@2C
zWA|=*c&1Rj(D_DIF|t|1W163vK~3Is=7U#v_xX(*pFpYpT3!8Fgo}vd*}9e781qoI
zAQ*g9J{no-D)+jvqQYi=b_(a~;a-PlFJ27X)6&;TGJFsAib@<EarMbfYvwhSgJ>`$
z-KWP<umB+I{boZS2RG0B%vqSj$w3+b^$V}U$HxcN5vYOs0&e1b1csuEK`fowmDGmY
z*}r};DqKTDrLJBH+TL*owg!7WN41Vf7e(v2_71Q(G@{E^Z2`U5A|au*BN>nBgsuoJ
z`n*<fRURovbC3a;P!Dmw$Vv<kB;SjM5Uc+BRS~R8(fRWJ2IMoEJ!S^(Y&?YjODKHi
zw8~6|-wvg^x2LCYR}@z(mncyM@wEgXFi`)-`4b_!#De3@1iy&L5PIu99<$D<0k`?~
zujAzGxO0pOlI8O82CJMbcw0#4$T$y2CMWM`w2^Pz(RV}}+p;<9l9PkO1a=4@eyQ{-
zB5oHaXBb*ODCJ<v0M4t}PNRRgF@~>zQ@42Knkp*;2iHe_d*DES>MWU4k{3V6eFd%*
zCHUiA>e9l(5DBYCAQB)HXb(XtK`h4^fD!0pAr0XH<8@v>dlsm_jg=m(4{$2Ubt^X3
zR7WSosErhg);wjihQJ+G(gaMsOtVzn@;elpFCDg?Qym;;JrjL3fbuz%b<)Tb7oY5k
z@N~4XL73?2{Tl~1u&@M-kB`57`}X6<k3b)|8TIw2rn8xd&q{!HsLb8RsYF>n-aMhB
zbA-O`t<6R;oa+nb=HF{;4=bbC2jv^NgN#(0sH?B!^D67o2m&vu`S$JBwQJoDf;Z=;
ziS^q71{*hCot}dI1-Sxu3g`Z{548hd&Qn7e2ol(o)%dEce7q&`bVlcoA3IGxvzpn$
zf5HLY3w1=yrMFrDI9tf1*4f#K!+}S2BY=q*ZTZW`LMrQQvNm4XIgF2A3gQ{sE)5M0
z7niALClhds_eBhjmsnX@)j+t@*U-{Ri;ni1PXEJKj}H$&4aMN%#RMZG8j7Ax`gc;`
z)_JvwkMBJMfX$mXZ{2Fz@#S7qQxka(V`;j#o?D&tpKSJQG?-Yz_GQVStSCJfEG%NS
zrDkV0m6xkd|NP4cVP7;tXlC<?mJaLm?*X#Y{7~C>?yMxOu~)7Wy+3{!;bz`FxqPGL
z53=_8Id#+jG~g3;(*_}yApn*bA3ufO22yAIP9AjmF6mk~CY|+oXqV^i#)d*Ga4jGp
z;QIBbTZ*d`-f^#?ab63@)1`~j;iW@s4kxv}osEs{W7pUouG0J4b>pUrf<>B>jSm?a
z{Q&bcI(rsE*hE^+TuFKP?dh_z&!6KU3rzp~Q(0S=O+zey-QRDHkW_nnJ2sw1!Q;)`
z%Of+9X<KP7-tg&DC$2Pb53k3sGdM#wHY0F9fnpl^`lTJV@rr+FZM-7CEtRkm&GYAp
z+W~)5RV*RH?b}C%D3Lz92E$v@lz#=ELQbw3F%8!(VI`?4DcC1UM~`k35iv+P$0XML
zpN#rrw#4l;xK!zH-yYs_T2Bwfdpz?V!UANkCb)q;x9H55yxY)((P?c_4vkJs^h906
z34xA=<G5E!YL$XzQ|1&5M!bqedE`DQnj{=ys{udd-8DC@+#@MT;+3MJjpRcp?dVHk
zRAWQdaptE)Qk_sibO3=ia1R>ATgrT~Gkylu&}^FB#*~eXjk&qQZQ4H3($USsAPkII
zy?QlH4J<3Tr$Jkk>(QRVycsXxm33RQX3cr<IK<*wTh#zW!TvkD`e<?qXzv#e9QX+5
zaZ|t$(Je70=KWjRmiIoL>6{Er473QQb%#&QiZ8T_LQY;EJ*eoP@gqY*bt7<n^pOVP
z3LYcx|2@O0tf0UYwaD3zrkojyfSCz=Xm4+ClrfyT#eeEOTbpzS$5T-2rM<1~v-o&$
z`8bg`*(^)b_^9Gi4h8Tqky1ZfuKimUZCJx(!9-*6cbA@un7;_Ni&F4ZiZ+fG_D-%R
z*+Id1hEj;k9qQm}U%2eQL@r+J!^UsG?yjo~m-BoyC>i;sdU~g!9xW45W(TF&z+>tx
zE*Li&+KTYgZ=+(-NR*HHD>dp&MW(6_5vwobRl5^_G0++3#v&H_`}?D%82M#z4e|z9
zNNmN^>9+ZWh4GIcU%{JwCvQL`x45_;J_lIcu5-mCO+j0F^w^7`y!RG|iE=A&R0J5{
zCgLilKDg)FlJ7#)S@VL}{*{MVPKX-5)Xxe<XwHd2l^JD~LY3EQ?FhM&jcxqbh%tnx
zDy#!RSA17M14;_YFJy#=_wFe=aT)F=pOSp#*RRTl57&Z{!B2aXb{x9KxpPlVcH&5g
zi)*fEr5quzmj<q=3mJYKyN~@uw`9ril$*j~M@L6_c?!UTx{7LJ+=PG9TNc<VwiCT4
zp<5xB)Binv*Kg;3@Id)4>1+GR#gl&m-|o3V4os*m8#Y`S6dKlZW?f29wmFZsFPqr!
z=r{>?;zwJXv5`@naYdx%?jB|kb>WXj$zbOwwgdhBT&%3}+%9mfpG8FQ-cm^5BP||&
zdcXX=mnO~i^_!6Ukelh>^~@$+?_a0CN4}zZdJ}&s6%FOb_52@RUS8o+PCR>TKUP%*
z4N8_XewADLFeqr3Rg*Nd0ssw89VGe8OzEE0(4HC_$*m;|8UhCGbg4G%e`e48S%`8i
zA@+;t-rZ!I5d|0&lHdCDA<fI@&rMS5IK|CHckFNhTR{Q{R|S(E?bk9s90af;JV(Tw
z&=EpXhyDf%EGl}N3_Aecmr(E1)8q2ZyMF5yO0a;iFwYL7C5;A*M5O7et;eIKsLPd$
z4<pg$bsHHX^&mhC4j{H>b@bs)>%`3N#>Q?H3%Psu0DbG&RCPP#9DHx%h7GQb48}RK
zkPEzKL$3(<MrN$tV;c=S^ZtEyvAU+FmypO36IGNUR6vjchfxcPJsyrtVdvlwvjt5B
zsJOa@#%f+ptW5y08U{Qv4mi{anlE4iEGE!A>t3D37Dk|G5bq0vCN{#u6>R|$h^h#V
zC}o5n1*jy%#UZ|SavY&h=Vg#_r_rw=Y{T&3Iklsc6YpMo@YxBp!KD#$Ru&e@B@xmt
zJ0r?el$27ElGf}n`iR|vUK`1>xR^gYq?`g2P9b(i?+^NEpL>@D6mWK1(GletXgUi%
zeZuBN`6&-q9DuvT!@18x*F=+Qs3A>6P_Pd@8kn@czJ4HM6sSBRFy~>&hlaWW_{c~s
zcwXSuV#Jx7jjP*iTyTKOKYaPp3zCAQL^335kJBspiGYGjO?wRYoay#Y!eePG2#9s`
zUd9L86|0?Q)%*=&g{Czkk7%V3C8g&>tGos%fcAVdZvYw`KR>kTXHx2Rre2w$P*FmJ
z_1-4IB#H2x*HKXk7!>L~J~2cy=u>hI?)vv^_pnf*ZQDX6Z5~IvFXKa)AmQc9r_4JM
zX27i+^~{WpHtDZTeubnkCAkTG+=mYeuM^FTk5V}~Iv!Q|_U&6|SC_7ij#>~Xm5{Ko
z$Mkr-i8nxtX)P$#Z9XqnoC6SRle1^l{BXIy+_J$49=}4}mQH`<zGRExmmS*wYfSAn
zs4)V&C81&$CwjjEK9CwXeB?;oJBN!WPOy>TuRrUa8p{ua7=^cO+bJe?FG~%g13>S@
zf_pBagy15zFLcvj-eUVTJGc2G!zMp{WCy}Gx2|7@NO;l83R$pINHjv>S;{nhktO>-
z=MmJfZv4S`pTdQo@;grCP%e<@+OcxwN+V0l(zcOJWdhVL8X0~`Fzs}pD1yqNE1E?E
z+fzy{!iXXkx!Z?@_fcr*JM+4oNIlK|vU%5*tEARvoJ1PhM|;yy*wFwDw6`O;_0`wz
zPyA+v`|<ZGqPhZDbK$}T`&df00zfQ|{<1tP0)%XA)ga4G?A_cbv+{~Jk;(>PsjUS@
z7J@g}7hKjSV`gs|))6NvF(~jFNh^ew9)yKeg)8>_QMM@0%q%TCttIVaRcVnj`=gKn
zqW^gb91-r5YXb`Uy1S8SWxaw&p~?$q8e)&0a2i>&gvcogRZieF`p1#6F*~R&Hb0g$
z+Eg!Ioqv+{b5I&M$&*NV+MB95nYr^aTkc2vx6Yb#<v$A+;Ry(PDs^QY=@k)tUhrqh
zQ-${feBgH(N!@zA7c1Bds&rpZ4^lzD+oWeIFGOfQoh&DavT4M7QB^M?!3aURBHY`<
znKKaf0~O(V0=l}9(POYoC)VjdO^AQXDB}rCqwhc>ipgph7eD!y-K(Go;%Ep5&!E|<
z3|CxDzFefr!JBJrAg$%*$|Ip+l|AJ%D0!5HG&bDfQh!S{DC-&{oqe*)<m!0(?PqI?
z_$KIgOkgg54;=t#MV}=k^xTR(_CP%sD+q$|?}=t1uo2llGCTUtgcjR?i_6@^<YTL~
zw2VyAmoMGF@FdN)LoS7cr4V&{aQzo@n8P#vf2}WQ>X{iB0-w*~rpR12pp7gHqr=0j
z^h>a{7v~V_F@N&*t@FogE5_YmRKFp0bu1{JJd%fZi#jT~eMkQ4=8OZBV99+kO5V0@
zTikyPIBTw%LBRQL4HRDkmyUZ7lq$JR3_m9dta0T%uB*%Gn>RJnyuh<+Rh7jZZ6Rc*
z@KyZf%j3>pQ&XpNJmyd$106acD&f>aGTT3?Np2&GH2o6lWbNnceEGh1h(`msP_8i4
zaWX*@?Ik8-X6qb&xy^&fDD2<A+3u(kauF6WHGfu^drnwp=lie+%w8wz?MR@IMV*=3
z%)-JQq;IUN3~HZT#&LD^z(<doFbuKT4wVl{u9%wNqb~00tSxmHGqi=+f$u$u*xg^e
z3fd-W=|>wjq%b%m`2czjuKMQn>+#Xiz--pg%IfOzwg9>5sVVp!g~t|y3B4n4rO6}R
z__yK%E{Yuih#Z-6WA;sfqhcxeP}1fHWP{3rh-G|#OuN{@=J_3L&ZMCs<ar_998gq5
z1b1ceHp6BA3M>y;js*n<q6vZ95o3zxgRBB&n3*6GXJIAxFR5a)b0YnUauQP1p?Ts&
zdARoqY)Yi8D}z7|VD5u@R_+qqwhi3!oqbQG<<f`nbKyL#Pv?L4FQxml)GiF%gkQg2
zh-F4b0`{0qf9;NUZRG(1KU-n1!x}giy8DQM!2;YQJaxP@<YaqC$NSI6aDI_QYC~t&
zcthvE<NtS`B}Ak%`XT7Ga&zNV2&ngzR!U8Y(Q2$aj|YxGNQOh7rk!01d={uVBkz~;
z9(yRM_UekZ5RZiA!;lai(8&~DB=4D6SU?$2BD!qTkxqxXnPthgXi&y`P*wf5D&vc;
zN)k|>mXNK@npQl1`$MYyC>Dh12b6fE8Xn79?K(gL*G5EKmvf32H#Ny_-Re6ibYpRl
z?4V&voq_wr*I)yWo}r-yw0aY7|BX=Zo1$B{`10+j{=)qH*7o-OHu0EOl9wkU%C>qT
zqqTI&5=lu(e(@HHTHlQI;?}QBcbx#&ZAOdjIQTu};luvU&Mh@1fVD+b=l=(#R%mKK
z$ABC0Vhs!r&qDS6_3PJz2Wxt@Uk6r3OW_;xLug0I$p>hGu$Oy3lMWTqWo*`Lr;aZM
zf{#Yxi7w0!Eg(`%igSOb(iPK3UZ`C80j-R(a)n<$&d_xClK|L{@E_Px)bzYIbC8G{
zV+V|jjEsaTa{PhlKB}<&R0!*0>#s#S#<VUvQX`HmUqp)t!$pI1b?+fuSfoAKY(z4E
z3X=mzzyCr86c8$w_Lzn5=#J19r2Qn<3A0tJg0MseS+4g6h1|)qMVyRGNC>-qdl>DM
ztBVWR){)L#llQ#It5M~K4nuazM~^a$8R*1e%I61}k-xW_k|^}?sw6UWDC^JYdhdCY
z-3VGYfqX$mHS%iI(EwmfDJUqYCegt2cfwgT<?8S6y*!<9KQQpatQPyTi6g2H$*!~x
zqUb+MEb;#^VHoqiZR9@>S*@hIzv_(#+g)2ELKbJj$jPf$w3;tpAE=`i)FLP8PA(kL
z{*OhaFZKATQ=vgYt9<h+YHFP9?SZp_LM#6fJxh81{8;h`%9Z}jEsa$3>N+}w=Taoh
zYYmNzMj^&RZ?R3sNaA<+=<0%8lGcV`acEkQ1A*9>pP!G76wQ8}lm%H(wKzBmTCE>H
zc>;;#hMylSt9{KJz{=a+Q3X_kO<7>J>5wUY{c=_~3`+;v`{crW^TNyzKaoeLry#N}
zV`4&J3(|$xZLrRpoDzc=MDz3W3(8u}MR!NAgW*5jJ*)lS1;JITR>9|u8=C~TYSffn
zoQ7Qc71I=d4${zeV!pw7NAb3n|5-*JhI--#F>=8TQMo!!4X6tU1zM)Js|3O`M49Mm
zHa0cYQd6s_uYcAJ4~r;2WVFR)tVAOI1W5GkR9D32G15N=e?F7v!!2puf_hSPfBQx+
zbi_fY_CmGBvFmA0vFkdnaM;1&O@B$no<VH2O>~+4{KwMi{a#r$XG6zDv}>R)rJ$&&
z4_X2A<gH#%r{QjM@3A$3rcTa>1RzX^l*#AMhf%ANwb7hRC-1QbW|e#E_(q~x#@HAR
zCY)^FJoICXVqlRESa+J1autSR6&0d}A?JIs!ViGI+S6QBUY=`|4F8ht0XNPjt!M&e
zbFF!*zkkU<Lt6ZiUbe(>X2FTE$bDDiAnnCe{!Mt(GXmA)c)u9!yZZ>Pdr4Ut2>IhJ
zm_J#@2Qe(@`R)7n&mp#W=FC2Fb`lZh%kAl0i(`|`1qFxpRgjyIOT+RL_(@Vy(u3U>
z0<*K{(WtLnwF=U#h<+x#VGL0C<s&zSk%n2|&cOg4(!@b8+~1I9g_ZlSV=;FJ6S3;$
zyPmQQXS>g!Z6xz6ot@{9n8ZLq+f^3X2WX;9yrDNkyO<r$W5E+e_tuQm4itsWO1Ol)
zyu7N)N=WzU;oaOa{r#D^n^6uTRZu{&56rv;$W#Vqf=G;~{P&4svlF>s`PdSdCcbJ-
zjE|$dfK=_hEWcjIb!znNw}o}$=4<>fAUY2F%GA>GB_iEu#Rmq$gM+tbUHKpKgS}$R
zXXfJ9t@qH(6&NLh0mYc^v~Z4*;e2o*ui+xkQ>MiJ++&>jNB@ZiccmF_i&{AglriFW
z|3*K(skrpdWicIjNzU)q*1omxfg|RJrp5vy8&(qB05;n-lkhUR1vi|L>N2!-yEkw4
zx(!WLJS+=FEzp@XW~7~O7h#sN6Hej1e<&Ig!y+&aP&HO;JSe(L)O$~!pv_{Q<M)i-
zaxVSrJTw{)#52o^@@iY}$xv>tYob}o{9He)zOhke+qR=BDv593{(=)-z2O{l;>}#P
z&gDaRSXZxH*%ahCKh=e-4)QaK%FM&VUwXuRm!Htr7nX4WMN#;SRK2pYJ~=ctilM~H
z6)&9r5>Oz)c(Edft=yF3X`rhs1&Jz~(&)fgT-<5_b5YE!8b&!C3lG=Il^a6JYFpj)
z4atWqjSH^393q|x{l4Ps*T$qX`va8LN@l*I*=6CJfnHb8FaKFwT;!Iej?Q1ZSs~7~
z#cTye(vwrejhYyox&FXKWhv)hD$F1L1p;(tl|bK;lI`e;7cWKVeOznf?{LYD)q?Lz
zO#6prMUs|0@d=jR%3|yrJT6XIzP1U{Gg@4pfJ5|{OY2R&dpS=~Drw1nAFWWtM_tFe
z$)wKBwfjrf*0!E}LOr+zJciZC)O4#%P*DeL+6Q^9LfT0!Sm0=Zde1lLr1AE_!HK8d
z;K^G}jdoPe0t|as`>%w01+f`QhyBOd#4g&WqGrxxivuD#R?V{PeqSwC?DEw*Xvpxy
zN9!QZIcAX_3nFA!CWeK>cVDQ)L>;ClW45yt4q#}Pll)wOsHig+E=&PLm7e*L7rnSQ
z{L8wUGn`bWUHvX|>x=DOD)a(oKEw{>*J167i<x5nk`$V7S6G~f5fKY`fwG5WpA|Nq
ztzOvJ7(TuK4CX7l$|7LhKJBcjsZj~pspzVVU`4VvqFg=hjH6@JSS=|jsk{##zP&t+
z2u-_b@J?&zvrEs_4}U4m%Wb{Kjq{|7x^nt-N!tj@975I@rs2U{=eKW@pr5P<nSY<O
z3)o1gkC>CXTfDZVb)d_}=IkW06!1(?6B>a)E`6$6-4;!mI(ObbI*GZ&d-oh5jUyus
zef~R~!B9E#+9aeQloHebs1st`I~Mv}E@wo!J45g<DK3T&8DD4_dtz~;IqjbrcC6aY
zLgSmpSx0ktW|&Ngdfi5>_{h!zjA7;GMm~A6zxXq5;w@*ruh!$_WHl8P5lKmH*Unt#
z#Cq;4hz*TROk_c8!rQ)h{`^rIs6&vqw_=)<CmUt}VLb})^It3LvXguWh*TZ>BPm_;
z=X}iXEimZQ)cRs^dJ~X5)Zfp>;Ej%zwQ_OT2JuJwF`u-W+(yy%X8r4yF60>xl%c$4
zNrZXIX99C4aL-?r=?->)wFJKT@%pN*kEC5H2k7##1tzu`Vw`GdC?Y5buEd^8moC)_
zg5F3EnS}jFt}}<rinrDX`okU){3GVm)K*scj@E<SKM$3)5oH}_<+}Mp&?1o{j*4pS
zC<<P-cx`1y)j3)1o@1gG&X+b_B6_&>d%2>ud%3nPmOlE8zFLkgraRi3n=s&utU%#}
z=3R1MZfik^ghHXfJA`5+xp%K@#tGL^L_zWRq4231<}Os!)d|Sx=<0?=M2OU#{XNKY
zx&8R$vuE36BZk`Ao^;<;H(7eFdbs|dnWL-k8Xu8RZg`N%Dv^9<(P<>0BVw!Y8T(_C
z9YV55@1gNqygGoANX+BMKalT0aAROHHaC}9M5#EbD7%9+B@I019FTR+$jC5xijf6}
zgN;>HW3cx_q+LegZ!<73Frz)d&=05sAhs&-8!F$24@Q;IM=-1fj5~7d7%F_0^T=bk
zXaCgjB(4IpFp^r~Tz7>mI!%^vpUHs|U*0qA<PH_B^3ihSq5I>YBx_3}J%XVT_?hQ%
z0HD&6!|*sDnu+?#YxewxPys~5M^LjZo3k4Sj&%?$Rtu%a1)o^f+|u~=Lif;i%t6>$
zTZ@W{mRDDEiJ6^1_lZ22{{<}bOlO!C4hr-wut{N2(WD0dUT{`$oZa|w-1qP>|Mu-M
zg4+}w@zks&%Ji!~zHS|Os!RTD@lh{=c)rmH$N)0=5*P(XAQ0Qp1C*~)37~Hjav4Vu
zpa2=T)izy5m^&AsDq{#r(5d>wj-&>~iw@5JLnDBn{r~)M63YBVOUrXI%{U=F%(W_R
zvI{{fvTC!AWTvnzTuN)F3xqN<1sf(lfS+mnmWf23ynHSk+Q*m;0m@(l;b;Y(dJWwd
zHdd>0SFW4T>_vuhD{^*!Zam1u_Up5giP_U;2+fPdCngRJ4j$Kd9sJ6oAr(W^2<j<X
zzQ%5QQ=<>5Ds+c29bYs)*#bT)83FOjN7{A|5+X36;>4Wb$pHat^e<UJcflFNdLLI&
zS!ng;g`?0_5>|sLWfTB(&)aog9k??^rGFQF32_(6e4m`09Fx&Tr)Kw#ds@+5qN+GY
znBt|87|G`_sI9G49Z(|Yr}Pn0S)7G{Sc&U-%pi%WMueddro`^5sh;gyFoaLBO;{Lx
zix(?Gcz|3^PCsfopFF<vgbO4})D>}2QS|fB{LWz@%A`VWek1_V5oA(r)AgXXoIR_^
zy{|FHa{<$qKFinsXl=E@$be;EY)p(Y3d`ak)Q^-h*-KFt$m`^ceBv}0C)ot6y1M!+
zlfpn~rG34<w|+M`kdSx0nFkAe^5hAxoX0Wl!<g^d>a~QH7DAk|s~_U5P^f~Fm(|6;
zyMS<~>p1D)nv_N#34*$!rK95lR))baWzeip3@YGd>@XmBWaj8MCrpxg=w~|HoH^5q
z&x6>|7Qsvmpvhm)0bvDq_Y}b%C~m&X*V+!ZyhnUy7HQ|eK+IC=;4(3ifc!B|QHJx#
z52S9<O@NOhoO$TrL3|*|ZhSV%C^>PDgd$`o{YwLHNT(j*8MKWU$#pPD-OR>-L-P~m
zSS55fJ^{y{)fO*@z_pOx+j4~9as80&w+RVlx=w84Ykv|T+uPqii!+P%aT#Bu>!8b{
z#si7*(U|8|=4UGn6;s+-P*`{Yy$l2&I=_6^g}I-b{a{O=RNc9A=cXw7Bz#f?^7uj3
z;>4=!RECtBYj%7<_G-E=FV7eV1OK2JI!hNGUS}DYp8Cu+I6aYD;83Bj!5}Qof<g&o
zEf&6ZfL~=_gq)}Ts}~rOfDW&|Y2}2CE)>U|ftV{;ser4DFN7}vu}~>-xKwHg+RfoI
z8^~O4bvzP$QYR&<#S8)y7o@8%^k9asNMBiqT_sGsq^KyN<|}dD3p0o`iCI}P8#m^a
zm04j0svnZ+-!O@Q`uu2w0E$C)`A(C@;g&61U>ll?hlLtnyr@KfN0ve6e0e0UQfUhr
z&jZ<pY&UsU7g1*Hmda&{8!90B4z{$2y8W^SH1!r=gruL7W1FTtzWLs3a6+uQQucfi
zr>vVZLQ3*;;)F{s-+9XA8hnhx%IlC_l*|s@uY^sGk`)51T(i47$CC;;dYD%;M%V~(
zfq0WrXmT+DbO=ALQh+;zgHYL*deh%OCQXt4`puit)-Ab^NH8;p{iUv@C2{eaP)|u#
zoMl*d5&hvy2pkx^Y^bcMfz8*AedT&OsN@{z*(9GNVFC4rMg{e4|9)fAJF0of9oH<t
z#sL5CF$VIKHBX*gp@7AjPJ2Oy`!H<-lQnV+sOQ<l#Ds)~iV8lN*yQAj+}swIzMS*h
zjEw=`yngvc<jk@%J~sf=th%PANR~;@z<{J$+2-pg!tifIhtXt%n#1iZxVhGNE=Rw!
zhDXjrZ*&Z}V+u#~!iBm33z*2T(ht{k{cQFW$+~ju*73BgR-9>)rB=RQ$O)fuQs$v=
zBP`?A#nWH8%TRA(Q_kQmyqb@nnOWFT{e4WBk`KMW=^Nd8Z1pc9g|}UKt=h}m^5ZsJ
zPo<)S&CU(pEnE0|SZ^D~4>j-xW{WhJIi?pi=4}tokp6DJmwP!Kr^1~13H?N;bUm|x
z?{c;_HajB9c7+tzFMSW=5)a{-Bo}keo(7`9hg&pP^X}2jK4;eWoY9%YNY%u|$;Pw#
zr%$`T*q`t<w5-cCB_x!#(J~E5r?2~a%5EvL9@xEZdvMW+^TWq+nh7TpzP@If^;mF|
z%}GyRMt7CnucN!0KjauYgLhGp(xa};M}A0_n(UlwS?;C4clCXF-vs}S`tKH^YcjL5
zOV6bg6cqF-a=j9A+8YkKfj0k5iOJ5j{s&vFo$a)Uy;V|-%Uc$rO0H<K#WWAas>kCj
zq}27>rlTfkT_}?-8yXtA*c@{qwC(LlHP$6~{Q<`}1Bz<Wd-m`hsT^?R@KYS<i&&?m
z#&#rd<#jLIMxA-G_OTNGXEl9h14hZ8O?EQ-UR_TgTofal2^JzPz3E70*O*vX$s3{W
zpZm>k^=8Sh$0sCAstv5$$BI{|>e;={pV@aKyC1?l&&>qyeHk2lDOBWH$4L>kUcEV2
zu4Q&Mv*P^Wz-XuSx+fBAbnbHp6eN7?5%D{)fhQnVPk7<5A)$N4GnW_^jIO#6`t*%X
zmBHBf_@*01FqiJc5vfZVMTuDkv>9K`CT1F%d+bjoe4YO#vR2+n@(T3Baw<oXFzn5x
ihRqa7BeeJzXBxM5rD@JJtDr6r#1X~g2VN=|Uim+PXTApj

literal 0
HcmV?d00001

diff --git a/lab_3/data/fit.log b/lab_3/data/fit.log
new file mode 100644
index 0000000..3e477d6
--- /dev/null
+++ b/lab_3/data/fit.log
@@ -0,0 +1,288 @@
+
+
+*******************************************************************************
+Mon Nov  4 15:58:15 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 15:58:43 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 15:59:43 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 16:00:07 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 16:00:31 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 16:01:04 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 16:01:25 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000 
+
+
+*******************************************************************************
+Mon Nov  4 16:33:45 2019
+
+
+FIT:    data read from 'diode.dat' using 1:2
+        format = x:z
+        #datapoints = 10
+        residuals are weighted equally (unit weight)
+
+function used for fitting: f(x)
+	f(x) = a*x +b
+fitted parameters initialized with current variable values
+
+iter      chisq       delta/lim  lambda   a             b            
+   0 1.2400030000e+06   0.00e+00  3.48e+01    1.000000e+00   1.000000e+00
+   6 9.6489173663e+00  -2.10e-08  3.48e-05   -1.668582e+00   4.777908e+02
+
+After 6 iterations the fit converged.
+final sum of squares of residuals : 9.64892
+rel. change during last iteration : -2.10241e-13
+
+degrees of freedom    (FIT_NDF)                        : 8
+rms of residuals      (FIT_STDFIT) = sqrt(WSSR/ndf)    : 1.09823
+variance of residuals (reduced chisquare) = WSSR/ndf   : 1.20611
+
+Final set of parameters            Asymptotic Standard Error
+=======================            ==========================
+a               = -1.66858         +/- 0.02581      (1.547%)
+b               = 477.791          +/- 1.272        (0.2662%)
+
+correlation matrix of the fit parameters:
+                a      b      
+a               1.000 
+b              -0.962  1.000