From 5660b379b6a6ec54e049d8debdda62df21a87967 Mon Sep 17 00:00:00 2001
From: caes <xaw5719@wmich.edu>
Date: Mon, 15 Aug 2016 05:51:55 -0400
Subject: [PATCH] writing final report

---
 press/img/PSD_1367Å_{σ∊CM}.pdf | Bin 0 -> 20375 bytes
 press/report/report.tex        | 192 ++++++++++++++++++++-------------
 2 files changed, 119 insertions(+), 73 deletions(-)
 create mode 100644 press/img/PSD_1367Å_{σ∊CM}.pdf

diff --git a/press/img/PSD_1367Å_{σ∊CM}.pdf b/press/img/PSD_1367Å_{σ∊CM}.pdf
new file mode 100644
index 0000000000000000000000000000000000000000..3cfec8cb0e7af1473e92287ce039656a47ba6560
GIT binary patch
literal 20375
zcmagG19WA}*03Ghw%M_5+v?c1ZFJNfn;j<|+qP}nHvZ(Cd+xp8`@Q4)_u9MGnzLrj
znzgD@bC0p}EK+%4Q5ptX7AVrXv+@@xCIWf_TLTLyZf*iPX=58xCo=-(Pn9AR0RaJ>
zsJWGsvBT%tO5e#?*x1n4$QX*37s}De!C2oK$}K%XW5foR4Q6dZ%`Ans?C=dBoHmlz
z4+u1M{~U%c(sF=bkF72+cW@&qS3K&Rcwftn!O_j#sW(Pg;d;bC`Q_z=?e2ZNXrLlu
zininft<Wbli<geAVc6{VIeJCJvrrpF#Jk$rmX+emK_6Gf@9LamlyDDK-b15RdhX(%
z$2|^-7`gSE3`%YH{We=J3Z~A*H58bvVY13d2~8B=1M3M`m?1}Ek*lbNrBS3I+d>KS
zHSCC71X6wkv>#%S3v^d+3ZDir*8DWQmX_uW;)KEOx6cRj3|TpNIQiYYPBVWPulNVf
znyOmF^pL>`WIv;dRA|C=V&p|mZ6E4&_I&30{M&owOY3X3Ci*D$*5e|cB{Qxqgb4yi
zVv4&41;>$rCB->&yqy3=fwX<H3{|NTdq7uH*~>f#?#$Es%h?;9=*ss61c(uZGGVA)
zKX-lVN%X8JsV!!-2A)KSKeOiS5F=G(!r^;gRgmpM!}AL>>ZkYWTAT43hJ@4u)ku_y
z!21R2dbalSvj7Ul+Gk7Z@lASByOBKpyok}stbJvsW>qsfWWNMEYwJs314#y=^9JOR
zA_oHOOZB)>BC$aU?Hfqzq2K~Z2BS#^<&j23lMczc^nZaY#eiUj2Vjz)kPL!-+{1PG
zl!?*nF|8c%q5!aXK$B)hY^>mY?`kk7)m-S&{;feX{eYc`GTAwe9)x4ZY*H$MVl0Pd
zOdVcOcHlg&j&kZly%0m83YSZjMnt8gY%C>=3gy8>cRoiHUPTVgm?)gVfu>Cpdw5AQ
z<u9Ue*(qp*ii@TY&!Z7vMRBZ$h`T%fX<1Dk!I&t5jHhR!P!Dlr%nm2hokWHKo*+2X
z*wUV&&q|7sbav&1iFQR*O)l1}{u*AlS!{*$j~|j_Yt9mqMS?`6-N8B#0fotx-)$~+
zVX>z}sgWX6o2CN7x#~9tIrqE~o!k~lsBH8(ddF)eE(%ro``l%uO*E~K80dJ4YhW!D
zXj5qS%E`~icUwj=onCUCa@N7@+LDokrHTZ87}MrsPq%wjf1sSk>(G$>xpg!RId56@
zyagmL3HRm)xh%t6OG1*crXKxDFEgHRb|%Kq+kVLdgTzHRf3u9UzJ$beg`fEFa&Plg
z>}cn0eSNuncu#9?>Ug<F9a75LALm@$u%hm-4Ln^{=ffyVh1=(^YHB2&(u5K6G^99s
zKDNBKaD``hj$fI$?4yQ(^r^y+Vz@7)Nr1w2+hkrd1;Ez&b>D1O>5`G0BZJb-7=r?)
zpd9oxmHQA<&H|x%<fp9IQL3-Q@&?uqflk?5)5;!3e%w^=q`$r({b1IkSHQEwJ!81U
z!ocJb3xl0(athuX&)fzX|AR^wy$aNYdh5_!3%e-=%ejv+Te3k}m^A9gw~}pHDBa&@
znp65$Q?Nq@(6$UrXfRYW`u#-I$P)Sp#(lr+r#rR!tbgj2PrC^c^|zqz&-pw|XOHIY
zt7R6M>Sy-4aF!`wvp((vEKd)ZZ}ac_RfvP7{9p>vC|rmYT))OZq&ARR3WX`xtrLjl
z0B?E8gpF_r#^HQ}fcrQZ<Do+BDOAEXqU`2r5s>8wY`*Z(o4I5xD&SfWIylp|Xv$6?
zJX39RaNrl}4jtK|KK;pqjRtBsha2VtvgGSIv7CDmx^ji@r!Qe<?!L`^9?7GL#!cyp
zZnKTSv=<<E0AmwMc@=ITB*#I+kRG_^JrN?3KZ%?C+4v(MvyRwfgxXMkiT^i5(<Tm1
z;4&=zh#5+X4+tH;vikHllQsCx%kL{o?XySuZ;clfRZFQ0Zby7Nwv7a@vN=!sPxwzy
z0Xa~{Hb(#2aX+uV_Uf+<`s@6;MKdt5GkqC+iH!f=%KzHDMckakl$<^{ZUS!ZFU9AJ
zPM`i$AfOWxAYdS%Gt~bqKezMGioZ)bQCl0Q|FSlrp{JoIVEHE<d^s?oVW43kVEsq+
z+3xfGum0Ko-9G`Hvb&uz0iC?QsqtU2$m=^8f4XM)ryZS)v5~pHpsgE$=BJ9DfQgl!
zfRTZfK>IIiq0cZr%^eB8lGERrNX)_3+3vqW``iCt+JDBT^x3O}zKx^ZS6_zi|EMJh
zoE)70isCDNp-(GeV;6HnV?{B+FOz?bT+!Ik*4e?(*pcAtE-L=lsMtT{|B)a${vFo;
zmLR{n`8$&`5HPT_v;Q@uzZHxG42+Bn%>Sif`pi)*|3^+eYf(|wTA61F2XhSt2ySTs
zb9H^_&gjknfj&Pk_Xlxxbq(F5?P5u4n{d5fKkzy5ZF)GX*PQ$QXYK`~rh;Mr)Um!K
zG?zd@(HMc%kseTzeK8P29V5jZ5K7gBteo~6oC_2jw*$;KerUB$t_esvfgBUiYS8xv
z5H~2;?3###K%R>ZC6NsO{%lvlQH)z&!Leynf>R4XfJz|sMK_H__a%3Lve8kk4_PO>
zeb8hsuH;UoFiv(4PE8H<05?G2HK3w76=w6p!@~~~bvHLR0J6-FXn<ZnT2-`rfZ-1W
zH`6u0dbns8eoq8^{@aZDx;jQs^+gZP@67;_W21g1zGh|^D?ru&NUHU!(a~@=02T_V
zAIk1@@xfbffxwwS1dz4hoE;6f2c3Kc(%oH4K#!($qgrKF02WdRgE2`zhF*D$-o?VE
z`qHQHAPhVf@(>!>kew@Eeq*7#J@SD?Xm9|_Y@|IHJ<~bb7+!4JT)b@%?W^9K_;Lg2
z?H>>3+1T#wLdB!KsPT4Abs_56Se*cYxOIqaeZXk~2mrvDsG8uhq5(`o0^HKq89o07
zW04cv%<AGZH}quEz35=}`ZMZ(f$nCDZpcOeuzL!1d-ra-$%n`d4c$39i~>Rr_+5Re
z=i|@2Eb!Pz@LrU_xGEliwn_hk1StK-)%%I;#~P#w+<kqfZy^=N3Lv0#Zold{obPj7
zivQC^c34;-cxRe_a2$4L=kOSip1zN5RA>V3<_o@{-^YD`sUl<jbJq1cA*HqLIl$(V
z8v4b|178L}^6LAZ4Zt1#c>&;SA#Yp&@Zd1YqfL6pY4N~l*X#M{1N4Y*-I;Ii2jj&D
zQ^$J_>2v7xhUl9vs$}Ex@;gV;>_;{7=9*mpy`QTT;4Xl#+wlD2lP^4DeNFY<D<ABO
z?on!|@7df(TsrgMV;i=<C-jQ19Eab~ut@hXuvvi-C}Z6#{>nR5wv8>AQ+-`E_&1J5
z^*8%y01Wlyl#dGm7AYzvY|re7k(Qc|NWinP_nRu|?fZhydwq`vpb97JS5;F}Mb3TG
zLa4#P4zNz`_Yehg`{+mqkc@GL%YPvR5M557+xfB{{LuqoRR%zYF1GV9=Uw<=w8#M<
zop=EaRQr*(4MNWWl<@~=^u3@FP+EuIG_(T}4<Ih^ZG-Uwl*xbfBf5hQA9-*Z^YLfo
zgIF3@Bj<-qWgD`p&BvggFWtE}+gtF4l%spY<VHf9Ei7482g1S&6Y0C)>V^1`?JaHj
zGn<E+>I?K$>w?L)m6Zjs7vKAsfbWf!@6|{3gDxN~K-|^FVVm<z$0Yv9pwT`2BQ4x3
zxYt7MBe<7s=RN+V_VBsw;5DuB$7jvnCHy0=fAojn#jxyWLBsh0{LbL|d4=2f4y&<|
z!1EFGu|=z-W_s|@_IzM}Wz*(s=ve9+0ql>*uk<)Dfv!hc_<~;{6>Et8lxf^mlp0@N
z)t#eD1=&1^UVOSWg}p$|TI(8lc@kf}zPUQ#bTKH`jBob>!w|m@k!_wK&}zt1LXJP4
z%}#B4RhGhK0+8oZt4@R5P-$>Q5@)3=H2QN;VPtq`ANsO$__0e)g8z2r+MBiGW|B-E
zyTfLyAw6De#l{1zxiTEyt9pkduEC$Fy(Djn%|(u2FjLKv1qWa&?VhwuSqrN1`ej-~
zEa1WlJZJ!f2sOVX8{WCj<`hrW=$A{p1E;eD|Bp;uXMRXcUgaV<quYGCZ6PTe&Pu6v
z_SS0K36lwtvj)I!?Bh9eXnk1gD?zp+s)%@lYx*IVF|IyCS+&|YA56hm!}AT%WeYj_
zHsM$VuVll?tM}tqqw+$T{E%QN&EYr1nrSzBVWhxFrj&(EFWSqhOm#1Y`}v!*bGYt`
z??7uQ;NkvoE&8c(1iM2IsP!~O3uSIPt!#0Qza?SMdKm*H%op{ScN9oSGIp+M)3b{J
zn`V{O0x$@sl#V%VVl&s^cCc}k(iGqCM3GXPA2TZ%vo2>7HsbxmKhRE!?SRdZT|)9C
zEfKg#EK7dxH7OQ;IPp@c;%R7l-O%{`a4&b|Ts~#>=hQjZTs3v*bUtisE3o%<!@mR9
zquU9)R{O0*6|?~KHo&f&|MsK4$#0R1vBMi9{*T~=ka{N*AIKtWPcwzmd1$$yR0TQV
zite8HHaPFb0W`<bQG3@Z1^{*n1>aAmT%Z#R>yty{LaQ^Nd~MxTK8V-p3UP!9ZvW%5
z7RThe2vCIz3q=w2#N>7U1p|&qYMYi|=G%Tz?GN>B*FYSGL+Zkja`shUFn6(i1<@1^
z@iOd1SU+~ZYhp^SFj;T50@>(-H#OjV3jx=2W#JauEKMfU8{TKv<F%f93R!sf{SO;E
zWT>gLy1+H4pHCqvv2L}{lRTKk?aWLeMib>_>uE*0!Fbo}v|8*{F8Mnxgg`=U!V0!o
zL=Gdbld1HZ8O$qv;#FoR_Gsp*zQqgb7UHL=l`HN=jQVH5H%>nTqqLQm8D<NYXI6wr
z!vz*ZCP)J7w+icoM*H3=hE%EeFHWnZ%oZA_q>12J_2?Q6@FvsCesGoOmyP_Gl89*7
z(X>wqj&9T0b@~OVB+r4ly)7=p(`3&n%nx1s_7+n@F8>jPrrQ)Q3>rhoK|Y<e(T#%K
zVEkJh!(nkP<wS>_%+?h{(f{F^^l|8Rd{MrzctS+nzNZk<OIaGB`obeK^AcN)DqZrh
z&^zhTQKn_<Cv$J|6bu{nFO-Utq(}?%r_FV$Oi|>43Eev0L^@gg=M>r_c+(Fw-?|^8
z@6v|AkV<IHiv7scLo#;U#lt+<sIX{Ly~lB2iMIf6FWu5u5-}9yuU32dfE<@;d8QJ!
zVVpdM4|I68<IDwz&5wBP?z=`C?8C#Y>_v7-^PY><Z+}D((OZUBQrXWM>^-$3q<ejm
zb=fl4BsDW9FdnF;1G@B~v-aM0@#)BXwt~+*J=gJkF9gpGm`-QSi7=rETv@#kUXL&+
zz9pE0Cf>|wIr2#+q*$KlZ68?CoY<fm#jP;Y%;H4}gUis)-S^Y?&CrG53%LbT>YJ4}
z`B)Y<@-Woa;p-0<=efcKUZ|C(5`=X~nK?bEmZw(Umm@h`B6qS`e<*S|D_?QW<$8`X
zHm05!s(=;-`WvMIqX`zM{|MT|7%7LF>hxR-M;Px$$YEGaUK1nLh(TqzaB_agAqMm{
zRDK?f0JD59TYE`WWsGDy`;K9ZJ8U@ZVlhY-Pc&0+TPq&8IK+ni{mM6p*}cg;@qwgg
zT|XrPW`V<~<NZ`>-?#S$fSeA^4#M2>3_&9wthC^68v;~2ee8%79(<Xf<8*kR0fb4@
z2@$X1F_u!`2U<};6W$miQF+tu#spHiSC*XuPE^ED0q54X-eeN0z;DMSP8t^&V0bM9
z-knYoGw#k`TbsFFl<MB9%b}OBg5<h9I>r#8o))g%3T;lguHR344Ch+W9=mYom^m4c
zGE0d_znlJGS&O9_09zSpNzbRCi|TwQtx37G9v(JYGsuEWi9OJCTyh5NxQCNXLF6_V
z6BBieBNgR@@-@&*M(5)ULgAI)UVvt*YN7Y(#&2k}T~LPNYgZ&vgFTqh4m(73B<eB7
zXt6jRTo{*K9G$epLE2563(QL-`RQt36;X}oRffz7z(qXgx}OPPvVI*Y^lmz)8&Ntk
z%kg9EJLBZg1~d~N3MeVn)(;KN^Caroc7R0e`Iu?k{>PFGJj#WP5KZ+3M@N4xjUvYZ
zfiM%yQC2&9C*ut3k;`VQhZ5|K%7{{dn-N|v9($A=U}C*I`l-IHwuzF09>}Obo@2S^
zmqflTqi%?ZdF^(bZX@;+B)JwHZWHszAdpkQ8?SLc3FRwN5149=XH^;3ksMN77vmOr
zZxZK6Jd8Xq_+1yeKhSaz`BCVPO^DM-L?_(c7(ON%><b{h`XtrnFMPXVL6E);kcm^f
z>KNB}*#o4uRuOH<y#WiAeF<0><K`~>gYu5PlrfL-)jJCNof0ky$|Yz6=Cs8&&wdI(
z&r$pfR2H5|q%UQ%GM>FaK~ZG6X}nxz;{vM-xGIMcjkR+F=8>4~%Fw#Ma_cn+Zcnrl
zGi7gKdbt=;-oAgwAYOJR=I4FQiLqqwhmJT_;EbYNKz43#xuyy?)VgI$;OG(JoR!dq
z?WJtdwRg|Vn&xm4A8RnDS+GIhMK<~4fPgvNf7rMR?0ud50>MOX3U9L;s?IyPv~A~0
z(_o%7T-fvLIy5W_r(y1hf!J3{Lr##s(zW_D{t@@t-$!*XO_8!)>=gKjTDaX995j9A
zy-BL{^|h~t5-E0C9FSD<h>8X{a;|x{Fy&9vkW9N0W+xlaUHmOxdNMkus8z@;*KGE4
zVqmU8o;XEF{QPA240P?Gx0D<?L~8KyFi{Y6aLy3CCqqr&Ybpv4BG?Ej)tRT=jP{yH
zxoyW|HQDcAH{-{b;0OD{4T#Cu;+q*+-U^p;EoH~CRG$REroGY4rX@cAtVqyAknGfC
zwj4d@Eh^s^w_X2IdPp4pvJZPMWcLQSp~wg4rek8uc6Doq+n6C--YlmbW;#Wawgep5
zsh11foV2KEt$0}c%g_)nl9-a9zMoQ3djPsl0&V5vzw+yjfn*OjaSNo)js)Cp_~9RK
z#htZru_C?L&V}~&yfl*l_TJrzm;l%G585&7=GAMKt5?s_gVf*U1YW7ObmQ3-&K;+r
z86Bi;iPhyDP6SJtSDc%>0{e(prC>JX6;g6^PL7?oQaw4+>Wjbis*(alRL-|%S#Je<
zQ_`pdJDVd-6$Y@3Zey;RC{ND9Odlg^9BaKv2i{UmFzP25ox9?yT24G^^E#O27UoM#
zoSY@Xy|2D#bWrE;%#2zh{>08$KZ2JU>530BY1E%>ghqc8NDT}}ay1P&Axc8xCv!u-
zL40exTLG!FPI6X`3<8O37dlB0;kqf3&0=OJK{E10ubb>Sh^*Uy(TWMTguEP$>hyr_
z1s7!Fgk0l5&w{BWk@=&+X<OiD5<|NRcsjaYlzu@}Z3J2n!xJszL@Mgubqc@nW|!cI
z<0j>nRn1WN+WR)>nG>K%hV15f15@R(=$$j}OK=#jA3z7)5a=p}THy@{bL(>_(L@5g
z+Tz<AjBJrZt6{Ak6-zs~2fY<|mh!gxZLvb-E=^RWC!v&zE;%-tI7lQyhWcliT;wU7
z`@y=xwDPc^3IF?+fQ1dzX|Mk9cdKyJ@w(er`6XjldvZ@V)aqN9N(VH?9mX7UMds3p
zNFprUlg6Pr(GLl)eW_aWTe5k}H%CqkRV*eCU*v==gVscDc<>P;S|`q5d-`rpWnei$
zxI4mHZ5i$``A+CcKeQZp`d(NSHeZZ?O!f(1IC2=JPB{Y!6XRuuD1Kak(OPTV>;6nk
zK1Wg5{8h!%E!VxvadX2<9R{0bsWG9)U1MxPSnA03JSSwX)3*<zCTv<u<FSE8hK0XH
z?0VED0rO4xuw41@`4zjcQ!aB&oG*r}YFUG}wZ!e6RjX0n-V#4Sgl&hcciyyVw>au<
zwzJR}<2}b_nsr1YL_DQbNAzG(0Qy&=${a2e@@)g^MG0~ywUOv;jG>LbmpZ{z)HfMs
z`M9hkL{^~=1sAwRyQMOkkpls9=qDqWm~9m`j6t?=B>`{?Rx_8pOnxUtWsZrQ@<6LO
zZcMKeM8@?r5_IAmlEQ?!s!?yn&@lvE5}J_^^p46jrOhG_;WTw#x$tO%Ewobc>Hx)d
zu|2(C5fuBB^MN+;j=uuBM;R3bBYW<|M*xSjrD@iQW{o0xFzd|Qcus>h5mFx&NsBdB
z)jcE-Sg?HW1Eueg+}ocS1**4jzKwQcW~Uq2!}5SE`UE&7b{Sz)=jNeq^q!)IF~Yo5
zHQKq{m=hZ~YXBi6qkk`=ZjiyxOo-k_)v_#b=RllmkVQGEK4wfTAxhyTRr*t9MSFLQ
zGuF0M!rpw`q)10oPFN#z6EBwWbErT`G;e%x{oUO>wu>*SV)hpm`;e(Go^!WKvQsfp
z-N$K5k$vFNeq#xZgY;F8inmuqJW}7axFAw^(_Z@cF>{9kz%SLpI7g8&9&j7SI@c~(
zN&IAa`HyGG^E)ey^yMW+w8+#sYljz+2DB=*6(a(~7CowKM9OSu>B(S=Z}&?SUD5*G
zk!$s7&ixdXo|t%6ujP`=`>8JRWg@^QA<q4)8VK90--qVse`GQ*0kon8F(hucLZUOC
zI0k4VLP9Z5>PuhLk^#b6V6P6^x}0)Qk7@;7{s_~W7d@9m<$gJEA&1{>P?w1Bqqi8i
zH|294>EG-ik^UZbq+_Q{D3F(L&fI%GDOaS9kS)T4qI#yC?G=HYcbKY88-OWuZu8;1
zuZ<;O{j<5g@P`tk`rGM@lk2`aTui8<q+fvbpM`|}edp1v7A0KF0n}-6nX+ELGiFN%
zOmEOSnV%dbCnm%F+xiPV>kJ+*U3PUs-!RA|Z%un$L{-9QUt<}j&EYWu47Zz~c(2P2
z#wb_Qh62J}(M3Xr{czG0%Moejs@{&*;rW=L(JyuHEq1*N^OVV{@G)B$u*zS-IYxR*
z!@zK>e%!E2zildpwF$!6I|WE$RuyL!_1mk|&yu*0I*YLM#txrnzLISAR4QTXj>r7Y
zRcEaT#%q7f8G}|#iP8RfgZ^`zBR%@et(Qd>5j{tI`qu?x-R6k-#h}LEyqc~r@npO)
zWI<(Erc($%-uHZGZ9)jXaPR66#gdAUB>1rmKI~)5?Hj}oQO2hdvb3{VgiKaU2a95^
z7D6V>eLIf&6b=hneBT%xf!J@y7*(fB&@JxCC^tbU>aNFde3K&YNOm)AN^+2jvJMG|
z9OdyMFR3jL>e(VxOS03Lj!pHJGfWEG@Waj|%hyKo!+Q#F-G$Yh0n0W&-Vt+D2(T(Y
zf)n$_Q3%TNZ7aUDQ%-w@;?5737$M`dmb#a1=Zey%i9DKwEYTewQEGk+)lneMOjOlg
zN@%;8?(Be*_(Hl12p#|0KVl)@|E0MLTwSKKJ6A4!!y<v1lSi5$2W<L1copc?-U!oQ
z*R+25hbNMw;2wEOuKGZ8#k4iGo^3Jp4<S$G{+J&eU_xkS#Jy2uGqUKWDMPbEuZ)AI
zJHu<uZ2|F1l^u7jl*R|wxo(mW_dghi+gWlkuFAdTsI+AxFa<$UQxA^Le?f_fo^8yM
zLvDNdc{1_tpIQr+_^BfETtXXIa|tYo962Y&(O-`9)dGl<CR{21L5Jh5jVD0^?^6Ov
zXC*)mG}UMTKZo)NKtHbc&cL==dwE=GV0^AR9HD4GU4*n!Wetv=HWq(0YSUW?J-LNo
z!RZJxYXhG{L~Q||!OIU~Lm^`+^6VaRyYvYrZm)rMQ#vo!3tsQkkTIeJs;r5Gcm|v*
z?}A8SBiX{74z!|o{H+87h<V(^&5W>!E*g)nAgc`Ob!>B2p%H{2+)>u^`0aj|Gv(JP
zsml%S%BEDkNk>%+9O`OmG!3~t?eNJ{VjnWUfI`(|;PdGTR<0tNQiP{2M>S3JlzyTG
zB4)UjKYfSAz8fzgN9RhUMTf&XL2Ap1n*jdHQWY)rv(b#XVviFLnPIryh~$=oBv(uI
zhLuVc>)@s!=U}#r%2GK98~B-z&L0A;sN_zTg+!0wd>z?JN@u$%kDK6RFmx&J8GnUV
z=o<y=)rdb<Tid0q3lglvFavlaa^G^X3W))-FhBIha>z+QUh9cZ9+9xw+@;&oU9pBT
z{k9ri{ji1`q3$U5)X_N~>0mD5<vQZnMtpi~Fz9pTatO7Mm<%lFbM>Q{ENd9`)T+G}
z_<j(}KT^1{a$wkE;!*TTR};c?bPxRTAY8T!mUMLu7AN_}>Y39CNoL8zTb4z+x;~YG
z0#go%R*AOCNvH`QLBIB16X>k9$G%2RQmq%hc;cAdA5JW$X{cAJFNiczHDo_BetKJv
zH<u~RVS!zAnC63f?E!p^3$#LXvAd5EVZvGYblt>YeIezq4+!lE$w@hiqsVn8oybs;
za@U6j(9flp=)Qaf7axRkyke%@@Ml`}W&h&k&-?kKI$$~V8c7g{7K9-l0`GBTYISow
z<Bp|2(**X$+l^rvzk-5wALNwyZ@AskJ|Rjpw82NU0Dzy^$W?A7d0nWYHf^;9+*Glb
zp%PcuT}SuiEtM+8;d3h~#WJ&+#Fta=wje;JypCo@twB)Wq{`lQE(^3aD}MOW_|5)S
zX#s$^h(hs<ghMt%CXMb{CqOVTa#nbX3kiRQ={kNf!q!q>luGw`ZjnLIUNjwWQpWEd
zpC<ARReC`h>R9Rb<#8x}w>^sGgvtHb2Q}BebW<+11V~!SM;3cxx(XqPv!m=%sp72|
z`19e2KAer9!b$%e$`Cl@@ux%qj`Fk!Tfds$r+jkAWrjD!qtOt|F0m?Rq^>t3h`Ya;
zjV5r!>(LUh->~s8C1l}<d{@^nc%TQ$N%Gw?w8|L<+fx)8s$1AkSinkIi0NF^9IGE%
zccNy;`?rC<0p*qAUaFg_c(<8LbgtxtGqz%d{51Tli!J$mQcWS8n${qTB32+5x>qFg
z8`OISZwM_dt1N>3wif4(4+ycM#1av<krRZjm};IRrJ~Jfk3ATET?Yw>UH%|*>U=fa
zhT%X3$;z3pVs}r-Da5=G%aJCoCdB-xL#YSt!|^IdeR0)xl<J}Q8kx0@?#(<fHl<o|
ztZDZ2yB*3$W4<K{2_ZFAGKJI_gSKL<#z${Fp!56@AIa{uO9P;DAnHw=k<$5LRlLWw
z9hEj605Y;gB+BQ{!9sE$$WJGvp*xt0M?_=WTNYa$->`BZr=y>uD(qbm9~4;y)vk_(
zflHzvr(6;OG+L!!-UupVLX2cYx>lx>^fj{ee;LiG%!X-(%vx$B*Z8|xJO=lT++_j~
z)FiOUJSYl>>ZBj(K!|F$U~w+)sBSI<(#6}s&f$$6!0Vzib`5X}w3;p{IEGz$QcvNp
za~JFsKX0<xwUUnnv^%2g3QzQ`5#ZP(msytWWCr!UH4TcNxZfG!uvh0!=tv*3+nj3u
z1U-XDChxWLc(-Y#Gm|O^VQ<ai6sL75iq?^T_JYmH94qnNLo*gw6Nr@src~RuY(CN(
zVHMf5&1g7~$*oDCy}2f-3}z#)%dQJipS#99XX3)IVi{I#qI|oYekUIoTx0hM<f#Lc
z&9Q!y5#p8YP)GZdsi>c&5}5>V!9gaPXGsEPrBiPt7#jlqn>r+i8k$4%n4NBh2Q)?(
zKqf3b7_I<i4Q$xhn$lzCd->+t@4?!+2*l-a|F$-8>aY$v*AQm^maxO2H#cj}pW(ls
zkj+V`UMK3%La!<jD@Bmgjo4A^S~0QOwuyJ8#Z$2M3d7SWkw>bPu|QprR-ooo)MW3p
z3Cw|{?Ri>t2U(krlRPY(V58b7Qyb4X_~CbW5E%6qrTn(btl_Togn7n$8xn%U8n__2
zg?-6YJ=d>F_S48OtVL+37nnE-00?=l_GsllAVLF1yktKBcOpu7^H2BuD~WY7`rquz
zw@pE2EbiK}r_iC-+A&mx1ngDP<&Q;NBE*X52uiM{d!=?1JzigPwP#Y6Z{!nz8IR30
z@08=IGG}ygmP;O4%DiNi<eBc4d5`BnG6knfof?qP`rg>U+xFAK1`bhJF}#L|JP;^5
zY>4<L96eUSC?EcucgHbf%&}h&FkU|Aa&((>8248MqX28($&SaB;GZR2QAc#n?3+fk
z$Czt|JLc%PlKu3%*488zQ;10D@(g~A_VMj>gTB5d>YR2Ib?Y(jF>fI{eec6bFSm@%
zLGx$@LKA0kCc)$=9IGf1ATpi2K}&Q7;6Y$_5R%-17{4Tx4YdG+@y98Bc}>~e0c3|h
z%0_E5`{km;$A}#Krv5wrduw2gmn|)oRLB#&Xe5x1h8K1rC3s4`#Km^4bm~nJt-J<Z
zddm~l$2~}k0@4bWPXl2YKa(*!b&qXcRBeOzofmj#a?%pkJFgn^8U2q}?%cK_*fgM*
z20O9z=bnm#k#QDYmXXwN2awZ8P78<j4Rzt~lQ+|orY{C4_(Vl?XR@Wobp>=lq3iH@
zvaxd8>DX0SUcA?-%>>RMnEBa@?a3@DJ^-PG70+Ymn!U2w_*8#H;p%1-a~QHVHP@MX
z^kf0zoD|owhIO}Ne@aT-FtU|AdU8XPu~ydGgkS`!*G`J1E30u67q*fSnw}Q|=Y**j
z;m2d6sxp9&Zs`ncLviUhFwIa=o##<17P_b*#uzwA)gXizlmkpGpMycB%*{2?)-9VZ
z_JElr%72V0&yUR}2W3&o$WdulvBfllz{usqX^K<e3$!LF*&j8m3-`)D-!9%GBvcJf
z&3n(YYK`XbKqpr^cQktHkaMb*muHF!Ycy-Q+;G5llt3_DeH&U1D6o{Mb)TLH$95UM
zW-Z;FdC;;cMa|!9%x>ij%(RxVwrurZ59^WH)zi_zz@(N1lQC0uO35ehB4=UrG&zPY
zSq`FbX%Pc}apJ;!(Q||v*5jrY2?dl4+M%R4ygu+`4Ek9^HSwb(9wmU77Ays7dDFsK
zT{2h!t~&AI0XC<$BRU@0m*_pCc8X3o2$bIpQCT+Qjw9fb(DA1=u;kucnODQrW6YFm
zeMDcd`D8Ye|EY<uD1$w<x>XMSq@@uJ7cn8w0%{TNYsetXL?^&=-AP|L*|3)C^B~@)
zH$5GZ{Jm97t+__jki)c_M9yJg69T_`pwhZR-<^K|C5~=uOwX{fZ+JjdacmQ^c7t`!
zV92Z!s&Z{I)N?ssafUxBD}`?Hc7{toIkNUSq<ViAcL~`|s&_NkhY25LeNkzvqdHh&
z`tK&oTMCDP2gE{C>?*2pKF`@13JTnV$lE<g&XT4$5Y-SaY4{RPLb~4F`I28(VpxTq
z+EBmoVnQWiOBZu$Z%PHI_J~fjLy<x0A!#psYP$PhjpF=TJk7QYm5b>9Knfjf!{glc
zL?ZZ`3S@0(XiRazf&q&B01#Yhz(2k|6%A@zD!SLaY>J8*jKKoiaPA&%Pg%6m(zJv7
z)z?y{pO}kK^03R7MAnF~d20I1SBAZ9bnFe)C3|&7`%rB(UCcq#LH0|NH`hCF?}G45
zw9>OU=w?iQW8v;Qe{o)y@psimV1;(-D|eVOEJqt|cu6&hGAyd;gbT;^P8q!4YyvC`
zS`=o2w)6KLM+KC~))i5LuJskS^*5kr$Us=AyNzTRSvO!dE#ovoy@U(s3=Q8e&8>>L
zo;4>2#w~K&_)=6OjquVW6~tK2ghIz?+D^-qcqNHTa=yJQE>vwTQPmN!NCZMPJG0Iq
zE%*N}gALCKWiO*`bHsOT%W!_O;gPs+cA~+IhUFpwGHz$&0}-^ZiHN8uJxSPz0zRfN
zGd9B&Sjf>_tWs(;A@9D}_=QuVvy<W`jwxuRk)~^I)({Z!`9r(Y<?_SRO*zhvz3-m&
zadv12`;)SdhS@2=4i%@&Ww=U2^N0fXd56<(vgCjwSqqBhTil2O_>_1B-0xCCq!Nn?
zTE7fCvUwGmmy>=MbYD@;qTn(RMf&TEfhw78pBePC60u9BnUpN8xO>Jhlr-3@Nkn1K
zh6dgED;F@1e#g?i&?Fu=>2#%MbV?ks(mfgdrS8^QTs^@Sp^01jm8u;1*wj2X2WeEM
zpP7YYqWOkKt;=yQ`SN&~{i4{+zPTA<#qiYjpxO`$a$L?u3jJHKC1gqZ;)9e=znT>2
z3UwYDHu^?WiSH35WL=K9)3Lo3ET#`=angP>+XU;`v-Cue8wZ`p6{cwJLb+Yf4qNY=
zu8*JGo65c85hnBj(~ViQQQ(-(p2sRR*zP}?A#(4kF(aCIk$%(Uj!ima1y_4nq^&g6
z9gF~DCxbL@j#f?_(88{fMx2FMU2gfXQ<C61=t&14lUObDJz|?2Gl?QCMQeF7R1=^U
zrMm|3T5b-f!#v9-ES}trPyr3GH`)zo<WEEpFkF>eHg!mh!<Xr}?4&GU3PVTRp7l+U
z^l`8U*Jfbw0zwv58@Cd3h8%>*I*V1RLbG$zHsEEcFZI%}0@#oY*M!PVJiX>>q>>j)
zm_1i#%?&S&xs1ddVg)z*Zoi6Jp1NWyEar0yck%MW<6+1!Aw?-fCz6symh}9UcG4Pb
zp27N_px=o+Y9j#s!3-lM|A9L{sa&0q1s;7{)UpFH_SyM~aa&LqtQ@s}AiSMjnBJFe
zn<?DSloz~sv%1CxK`(Byx!<;sbLssC0UnQ?tNrCHK|I4u&`FBt`v<)G699UyP{GM)
zZW$OXA)cMuRrI5nT^>`DnOyQ^dnEpMEt6Vks6WoJD&FLrOvcjaYr@5)T)YTLs(jJH
zg}3KV0?8scH<gP*iF^hrp$&bNeAt}*PcBTdaY0-j9Fx24<LI~r-uxex2d~6xQf$#L
z)jU0&{TF_eL5k(b3iO()HwnZkhZ57|ni#tM#A!|%l7@ExFG7^Jk2CE*f^(TMUe*_O
zXd_q4kV=YI`OP9ayXsqa6)P{^hNyN_*gauQ_koq+p_B%3EL<NRVM=~<y11}^mnyPs
zMMCOHEM}skHfGk}K)gVaLZ@L?sB+3NnM+ZMl8;j<qbplp2=5g7(a0Yij;0*ZLo&cF
zRwK#RQOQ~>)0afOugHueii#dx%4<ef8ClFtJ=#B>&kJ}e-9^I{shhEtU~e9t;3V$Z
zwHND=ucMPS^a8mvh&2-IM8Ip4CGC|ALL%x<(vpXp`YeE|^4xJ1PpP-89MLOW1}qyi
z+N-I*S&M8wVq2>vZN;w4lghFf9!h=Iv^0*yv}MU8x2Ic%7iSuJKe9iR%(c3B^;`XH
zVTbZtv5sBsNpnV7q)ZnN6RO|7gi868uF}x@Y$iM8Zavt2caW#Hl`;rBjQm7Y!&NAe
z;!Y7a2%)&kqA>5AmpuAvvupkEQkk|{;RMT$*ONZl?XcS)UYUYrA(Ic2sS!WUylIoG
z;@@#HF0p^~5@rVfTsk=ZAyh|K*TERcU_chm6JF_bGnZOIdu%HIrNs?5VTi7I>IrBB
z{hMK|Sm=QXHebHeg3{o2bycIFa#Lmcn^;rRwO$qh8x~2;c7J;fTs!=p#6glUNso#+
zFEAdWJDZBw^;G&V7BON*<151U6Y4jaibdz9sCNd5Ek(zpq%((oMv&(hC*^hcceQ-^
zX2)Kia~lz$?bUbWmZPc3lJ$-_dWf#;o(f1AlA4%<7+jHDh3!LcS<+-|->9zo+tqcO
zZ=QbuwyalylB*x8CG59GelK?0obWLkWd<BNjKOe4sVWViJW<W-Mr@3OpVB*wr=AvG
zZiSKuDmZtsof(h|tnk%Z5o6EsVRAwq44)%;$|Z5b@Kqu{yGgZ@p6ghYB-jNO`XB$U
zX1HT5O(BvJw?g6RbmHZJ)!dhhh?AnFx};~`I!$)IJlTBSaXRse*5xrY*6_D^&g24*
zb}uIBhp1_JCIBB{`Hj%PsOZe@yD{z8t7gv|-F!SZe|>VlGSr@j3ed|LVu^z7P`=RK
zTlKC-UpKg=H&JYZMa~Y183+>?jp^zHP9VMGKiCJGEL2eQds^-E6!8Fl<_sl99nLg9
z+rA$LQQgDJfi=K+p_pB_O(tN;QCDz0Z{gC~EvGqZ11z1EWsClXh!@gzfL2ZsKsbub
zn1=6W?>pX<#b&cmg{YSe5hI@@_YWrv`oITdEcBi;>gdD<*qD5-SZG&Cx}60^PABrh
z;7*u1enF#$%|)X!C9h!x%*i&(@!&2hX@W)%keg((+Bf>=sI@|L*+WlR*JK8WV&mu{
z1Sm4Q<VfFfx3U7zl&O<n411Xx2pZI2JyCDdXRV^?euG|F4V1h^XXVm8VW6vT^5_fi
zkdQ}MtsLa1e%#n{VT8`JjdkEbdRVrXE@woofPirvnvG8>fcMF~)KfZz2*tmat>a7k
zmt2bzGb@NBf}<Q^p7yms!9*)~Pw>LgTT@fV1v@U`T<z9aw^Cnj*tJdhI9uP-o(1!k
zS$#}YxJtpB(k<hW8vTmK0Q+3PP13M7qOstd4fV}(*5U8Dk7pSiY)WMNJgRGRE}OB^
z-({<$<vbF1M25>vBt>lC#0uCr9(&nOITm7^#NF;1&o;4FtHWGAGAHwQi`<MC+Ztg7
zS_zF%&=0{W2{zDJknxiPqxCPhvwf{iJaMMSz4=+A7$%H>Y_p{{M0`vPzJJU;cyJx}
z?n@Tc92bs;ZKR$OWBM=TN}r4^!aLF2#b9*hpGb__mDChZlhO}SzCut)^azN|qs2{K
z(`Q8~#7ByQxhxdk&^+3=#B{4Wt9Bk%F2O10vROThkQd><rgHZ7*PykgDgFri9qDii
zUCOpa-OR-rxI1{e)W57_%L@cYQpyf)P`pP$q6%2X5g_(Qp8JRscO=ZvT1clY_dAos
z?7P7{4XmgEIE@dtqovpikFcvoO4G*ng}-~8(pc0(YN}S0Ie(7q!Uf9_tN5n8dQ1r9
zfcAqFD+S@0!u-#lvaYK8qB5;?(DIt{-i!{D!--lyFMAs4dOB%#c{_7T=JmZS#`|B5
zS-HhDdH8%ulo!iKaa~^|U`T?`c-t(+x1iZ}rznp>^E`8r&5TV=8XCbL#9)Uf;z35D
z!v%h?AUBXF8;k?9?vf=Rj}87r1&jp}p3;w_h_@GBAyzp$5RX|j<w7`}ouXuq7C|i7
zpU-&7CN#b_d9;|R3(GFLzDK)8L~)oh03(0E7wgM2_O#}8Z5cMWTjtokJx9V3LIN!D
zUQeH`^|*dX?bdN3%J2tf`SpGA7Hdt@DR}dnDMzhjS99prYYJW#B&(*_*uwYIfk>Zz
zHSD?COHR4HWm$Bv4Vx)a-Q2df^dEx8I*{u~j@$AE$HJSaFyrLHc;(|<D;Zs#uWLJF
z(3-OZcU}upn{RJ_8Xuez8I6IMmNE1XrUe0Ol5^}*?}AC^t%jJN7&!@NOY<!uEKy4@
zS8J*+2@YAAU3H3LNLC~~`O!-~i^(HN08I!wC)-EgCfTV}?m(dhp?5fnh4J$XO!toq
z8^=1c_4{3?BokD4S$%KW*k#pcU)`g7lkfrnLMv2Mc@CrWYeDJYFqNU~{DWAXcpJh>
zY~&A$nNxHNzo&4PlNE*mM|V(@(}molz2~h)pn)*yV#3tB-oB(Gppl-5pPSPyhu7~P
zSK|FyY&*AdDR1uIw)D`y9~u(y2Wd5<-9P=N&pY~le-gvp?ERbmh>Y^ud@oUM*0!P%
z-!>i^I)mbUr3>jwLG2X0(<*ovQw?b=PuONW@bVI`2xO%hNdf{N*ecjD@j}B?i0nW*
zN=ZQ{{L+}Q<AhbtTm&i-ylP){7Wj!z&rh&mJOy;*j7W<&wx{Qr8i#>CMa`E;g=9tG
zC8gG_rI%y_msYW2MYO2N=!t&0GEl77(drxb!o4aoL+8Ux-9BnnkL;(<nysyZ$T(f|
z%+G%waNT|tZC}ob<-M9YLDHy&r(Ggiucn;8MLQ~slyFsz0CS~8AYHt2F$Ret$)2M2
zIZfY6#n?Vyuaz=x-T?hM)#vG)G&MMeY_jUTJzIFPmvoUy?D2BN+pUDglASm>@1{(f
z+-u9p%D>ay=eOBTwnS=3NY)>?%xa^PXG4ir7K{mIrBF0Ze7j8qtQ~}U=?Hr8#}#+j
z2`K`$*_(>Cu+<Jb-RJSU*ptWYiv2Pi3*brNb|wdG460FM&kN862`;)EtdY)MS<v0J
zu%x)+hX7l*(NK$7%TkQ2aAKU6cncP=wr8m|8z~Vk@$|c!iN^-~PfZ4pNT|liu2Hbu
zce-H~ch2!R=3By`ku`n=)mI>=f>_*q@s0j(sR}kH4KRe`KH~_*Lw#k(BCvS}*Kr{~
z0`fe@{=`O(MsJRV$Tm6N*nfVzu_y8R$)ey6%N^mtbelmviZfJFaeVAeOcpg1eKe2n
zZ+N->;28eoB)967kEMgMEQqHpm0IClPd(#`I=aAn*%C5OA;zR7jOE_ZfuY5AN=&3s
zL(X8qDe69K6Xl?nyK`z7LS*JQnl(Bc0o(K3<4H<lpQgx5kZq69T)=?jqhqLxk?&__
zWo1OyM~j%KTiz*m&#4)dP_1V*Wmw!LlN^{dLd=r2uIiD2of2>92S0!*0JEm5gh!<Y
zW@{4TG1lh^1D!+cqJ(w0JmvJl9zJ4|S}r(AQGV*|4j60nVUmi+!+{(-I%fhb1HNA?
z=^I{0yk%o05Hufa!pC8swbE#nMkfSrruz9s;JO$ja&$w`$<&!peW1=0AoEoyg!7^b
zZ&ZWUn|>+W#v_d}3BMQp=P&Zq1jpU-h3ANh7a=yNBRlyM+BR#sf{aT|*vL4P*J@uM
z>4)ee8BB_jS5OV1=GTNL&&3enL{O;Yx?c3Si*E3zchUp*lfbB48w1ctijp#@$_W&4
zwvmF?jc-0v4h(#YvG)Q>GV4ikUauv<a|X0eJh{>%IQrWufh;O&+6$w>$yUZ2(ERg_
zQfWJ*;rsyu$9eiYx6q5WVC%X<Y;3E>Nd1&7Q5(@xCWJDZL9jd#-3M;Ddgd5ErDLDJ
zd)YN-<cK(JR0oN46WUPG?|WchaxvtM6;Sf&O`m2rD!RFu1g0_v_>(UHV#@-e)4(5S
zUV>aklqVo0K3+zGxXU&MFC}y`rVuNY>WYfFs%^*5Qy}Qmc|xbAIshxKL5QhP4nC!)
zKi+5sk8buEe>egQ!AC^UJCWJv$8+7r9G~v>OZaI;cUkK)WDIb8=1;MDXmn>b)j&59
z-(svkUGARzsPOqT_WalOM!!4v3t2)TS={G{1s(2eY3Q$$1o_Avu3OP<V}%FyRMvYH
zx6mwK;98Tz>~579z8wk5;1Lv??`Ztg?!L06_3;75&qWni=UqC?mGwMyO7$D_YFxvh
zxND818vO2jF1%}1fQ7GXBwQ9bZyib91M_R5m}Xk~1gO7@W6BRXOUl1N(>=R;vu^Yg
z*XnIM<b8egvOVg<0P~luzwKF}e$vn>Dip5rj-|-tLX;xDGmIZE=>sthbx&98=oD0R
zYU>y7H4YJIQ^4!Sv+MboL4QL`bU893BD&Rx80mLow-@CPl}&xxxnJJPIMso<z76Mg
zFi^15WaSfi1}iCvx&yrYu8(M>F`d&%_`DBXuC8$~2~qarz@|%i0DGU1x}>2Yvratx
zhJR7Q5(mKad8z_aJrP?HLofsfMetBc1bT)8hjG5qCxxN3R`5a75XRZ;bkf_n@luD5
z3KaXT`C~UegS^CXsJ@12Et6$YFmR@E2S6Z97BGi`_JCV`l$M+*536}eeH~TAbA7`i
zX~?8`pezhYl6roqEWS&vYPyU7=iA^6qdt0Z{cHBAT~qQIn|7M#dXw}R!n=N*oXMhC
zzK(%eFkaBghzgqP(&GsL%{O*e#Q7pHDdYeGbSDVtVxTAtl6u9e2V0GAUB^Eo`amoZ
zyz;#4Kou2^PK5)z8-zPT`wxFn-nFxKSs<5*!sVkvCn3PlwOkvA>AlTra!)?j+E-8o
zwUo(4qN`QwYsHU}#<RP3FU>z1rzVwzd6F~?0-XIyM|fg~VH0$Fr^4)4HO}l6VtF#E
z(R14Y4ThZe^7~QiSpdncAdLOJHy}D|@rfn~kY^bu!I&(}5PZs%bpqWbErFaIFN@?i
z4Ni&bweu4unu&6-cJpvUJ`yAcSv#gPKt^eZ$ozcc)ltmT#~$?i2m0ejkY6#)<Ra#-
zp6RI-@WmO76Cels$CW9Y+v5lH@~W@L|ILei@k;;Z#hBPQ{sFcazR;`x!8S4chtc}q
z;HFPp>fiLQcI=2v4*`P6gBR$C6XKIUph`8h$`M%v?`$o;UJrZStX>3r#H%-^O}wP`
zB2Le!9lg@&{lp9XpK7<CHC&V51}GoJ51x>4hy0_J5L9U1cIQ{rd|H$`!iY%m%MTJ{
zU+Bkkt&Oi@&?YQfhsWZwz&#KYhqq}w!j;&E#}>~0ihGPTd=7>66Y>u^kIbd_be{L>
z_+qFXu3au>o@;d)D9`A&Z<>FV-ME61vUzgSwf$ZUCH_%Q&9xM5NKWzSYbi-+jIrj5
zFQWN^@;%WU;m>r+%&$g3wo$0EWJ|j-E&tM9=O3J;6M?9XGRcL#{l>CF+lbtZnw2tT
zF{K=8hPq=fte57CVu36*{x$lhfi)Ps5hA!F&eZ8NT-fZ`E;cY0SkYard+rgbtP4RM
zL+_Uir-t~5H-o&uUSQGwY5HstZX;J%Hk==(r=@PX!8L{cbSExUUe>F+4?cED7)n&T
zWTQNt4Sb-@4@Zo2=?Pi7X(=sN|M&d<HD~{tU?w)Ue^2j!fWUu%yf45{*x1p~!Q9Tt
z*5NM>?Mo%AZ~aN@sY$4)NGMaunp+z<J1W{*%i3~C(<mC7I{$5FrEls;@He~m)kg4d
zdXI*Im7Rcwk^Pf(qi13u;9z9_7sMywq;F+zC}3l1WlTW-mj?ky!!HQ%)8sGlD~>Pu
ztMOk7A$>b>V{=op|ClH_8C$Cou>alsOA!A@_=@0jeHs&d;d%d!>wVFR|19ahA;B-x
zzx*pX8#w)iC4TAu?ov?S(fF(W|99}_4vtPjX8H~UOpJf|lGgv1%)t1!r!O!6xwD^~
z-(PWkv4A2rpWM8;jp=6}<~9O0j^_U;{&J*dZscU<_=!O>vV8qD{_h796Vt!yS?E83
z%>T#de{WAuz{dC$5DUu}x%oejf3#m`*1y)$e=NRw{=_r?D!;VM%wO>7SNRpi-_qCB
zKL;ZlI|1A0=YPid<$;O)Ya{xP_1D#3z~|@cv{kV&|9ZT@=PLTN`R@s3_+Oq-!0`Wm
zM8SU(n2f9(tp5jJIi{>(r?h~{``~jRKJQbfk8MbkxCd%DBF<lBxVfBcU?+$p9g&f%
z7g&i)SH<8N7NFQ&nT-1_adMtS8hxfuqVPvx3Yub~Q_+!-beuFk7M>Y7)2{;SX)=)V
z7pY%I(x%G`_wUcn6Ay1M9Mk$phNGQGFwz=r_H7sZyQk!lL*9`RVk^(fp6~IGBc%Q8
z<Q4w+=QP$r9k8z{^M!-F>s19T{R_&wF>eI+HJo9U?J)vXuh!j9xI8CpR0{P$m)70;
z=j(HS;#<-#O@k6L=f}KTRfI8+F@wCl-z=Y5t)|{G>=2KXq)n{9$2R<6n$5_RQ5}2z
zHlv-uD6G6s?PUIX?^F=Ef!CpZUi)$;(6P4BztrYD*yJiDb^Ismtxi{sn}I>#CT!!x
z69sqA5JuO)#6VSb^gX9O)zBq$4XZ{S7tYf@MhnAn5qdh_px&buFJk)+MB)a>p?&+#
z1FyydCq%CTQ^Q0^Q@;~Bf5CR-;n^V^3i{N9VGeh+hG5|M@=!fm?%;l@;M-18zwhrz
z-3B#f<wYi{)983rWogr#lAnw{zOC0SC&OYvRKq9OX}95xahUTFT*W_Q?a71h%#R7p
z#y!UpX{w68V;T?+k$Z7w@fv-jN7uZ@@hP0xb9^BSCJS{CR2p^|;qF_{)zkN1AzRVg
zsh8GN{lQ{1CzfBt{F~rcp<@e~NOuiMhd?T3%y3=`lZi-o+$wE`PI%)R6Y(YSc7MUX
z`|<VOajOXp!GSscM#T!Px4IhN;)bt>M26GL;9J^8#Ks0T-)~O`otx`pcNv)nqo-l?
zrVK4;-jzk)XD8pD`c8BQ(?aFKC?`qL`si`}AZ=pzP}QH+vh|5)MuD~Mef0%T5!}S2
zy{&YL?H84Z#mTR`@#2J>A}MbH*TRD?5w9x{@YO2jFWQm$dnEn4GSZ^Jf3T(z1Fb-Q
z+C^7^r$#sT?^ilnQPiOs{(p^}dpwl+9>;grDq-B3uoD@H8HO>>T%V!KvBbJxqg-Ol
zxQt}R7{bXospJyL6tzoo+gL1OsYn^7-4lC~h>94Z%OR~=$DV3C&oiUkwpZu*^ZR?g
zzxjPHzt?=<&mSLCdx7J<mMxT1Nr2_8hxyR()r-&!dp&riujn5W4z|&U6Q8oGI<p!N
zc6*oME?jbnwtigVP)Wb#IY&~BuiWjMq1kBfOfEJ>6lV|WRtuYRkl&Y6w~ZYBij{!t
zPX^653VnVNx4NN9EGkjVQ^{$+77HF!+_CVW_wU%wI!y39c5iI6Uv*$L$<ZRIz|7#^
z-Zo{zEB<#yW4mp?%2R7TVRrIixk<Np@4z~hJ;8Sl-mKH{M{e{k3ZK0?JY`-u$LYE7
zB&nmCF?dl-BLcU3I%8~J)}0|954u56T@&>uX8hM)_U{PiHc6m)MbepAD@U{FEayfi
z{gU)RH$R5~hiBD!8=K|VY{l*0RYc@GiwK(f?}Y)4z_uO9GlKGR+PGYbap|eJ^PHp;
zd0l5FREw@1sJ+$uI_KShoQtK=k1hvudyn~>J-yf$kT;ha8DuyGI=-Dn#-%*hx#8S9
zqVQ6kUM%45>ezt62j0}mXI`qFb^K%gv^*!k6v+{D>(#VM=0kKIJ8M%UL*L<}@yBTa
z+Uv#<h}wPgms5+_@y;Dp9ZKS=13;_Vj{X-yP*5$g^DfXOW;WoJ&;Hgy5U6$hXB)Ww
z2G0EFP8@1)wuIkn;<6UFbPv1c@O;#^w$jlDA>dD41}^OQr)8V>-W1!*W3v6k>F8|X
zNHi|H>W#MirWncY6im%GKfLj?(n+fK*Xidzs`J!YH}*n56P)UlsN3&Cyy>wj1zGLI
z>gko4^Hv5bdq>wMj*nf(&Rq-Acs)H>JDb#a-zWX5TxoD^>Y;7!03qYhmm>|9*g2be
zyNQ;{+F`4y`)sY2t%DR|v{m7#<@BIc`GL+C!Sp6q5ltB1e#FMqNCVXriwpLno;%?u
ze#z*piI1vs-uk%Hd*eCh;7r}`vIZM(W2nR88wk%%TRo~%d!nN6OuGf*E$Z&vo=#sx
zlIE+*R$6vg;$%kHl^&ATa6ston)8<JE!|H%5-_KyS~D+L+KTU>Ebf<<Fp;|J9zN7W
zPP;c&wVI_-9yLzlgMsQe=IK;#gc+{OCNuSPNW6o!S4T$TL4#}zg*;GV-@!u!>WK%)
z<QsK;M=_%C_#q?ISo~0`VSfZ-6NM42oH<MEEZOo*v|C}PDSxa=#WW`?*oqqQ*mCE7
zt;qpwaNBxHzma?Cd0$Fe#Ck8ezP|0$*#!IUiSpp*%rl-78Zpn@e@}Dwy6jdyYj?b&
znwxeLYZCfk=TW(a>G=%#26A_SM1w@>P6NftTALZAX`uSKC^d$GNcCS~S3aL8Hg(^d
z>2}ud$VIg)*Sogbx=-5fNy$eAmypxhhIpW6lOfqAUCGd-$f3qHtN2odKhmc(4TUqY
z)4!@nL3(Wt|EESfc2-X1<Vf0QN9A=CDz=)-2?Fz8*hyaAKWeVjrZnUzZjFmyKg2v@
zH*xuMr+R&ip0J@y{K%H&d3TD(V*gTBau#J|%x&=(4a?`>G4{t_@_Nn}SqZY;2ibq0
zMj;c8DRu7enzmH48Y1f4zHrnuHKFYzs1qt1T7tIdF?CY4iLc^LN>U!G6e=hM=%kF<
zXcy%E=HcdXPoYxX=`eY-IBrVl@Fz!a?poC7we9?U@79ZkSZn`cu3aj0K_@Q`wo8Y+
z1<@i%JOZ8@;J^M)a!XEoL<d1Rbw;45+hQFpj8Ql~muuZpAE`Pz$cqSPm4e4w-y&&z
z!XZTOyX48{SC?D(Na7RM?53ZUx92bm#zmTp>p;#A5#094tGMPuClB9IK_SUS_m7Hi
z*SyM4i<#V@GJ{u{|NgVTCD&Q2kVTTmLp$u==Jv{A%9JtvO`|<sF+Rem(64*sx5v4>
z${8pPO3t34==YV^o;w-I2Qkc(h^&flg<eWo`4Jm(UtgH%8=W5JM*EEYSe1x<p41lo
z025d|B{EI8Y2=V&_BHyT8BW~jKJn%Y>ltp|U!)OE%>o_}q2YT$t^-)D*Rgzvl_e9J
zYh`%g>+9p=YmKuBiC}ogdWP`Wp&wqbW$b51CZMshHMDfmSho5RYM@c+XfqO-iUz1O
zIvM~7^iN}hOi5!A8%PxCPb^{syp7cKX>Nd&IAr1?i}EGFQl*cNp@X%(4V)x+Sq1`&
z1o4slKpvh6F%Y}hM}>0W)4>8P5LK}uj)qcZ1xdGrNOJ|Zt}BlPB^F`@tUU{2BQ`(b
z!;)t#mm3+fA0{R#4=o6nEH4>aT2$FJ@Im_#Ag-W7Q#L&P2Lvn;HF@6=@L4_<Dr}dL
z*C%6qPc&%r_s2H-Mp_72u~-ZU)}XAxm?O=yhfL>X$<!|snORbdg_d1BAgm=5Bi+=I
z&4PD=+t=fRVoWTQ7hp+*!26FLQoS2UuMFZ{G?vL=^EhY#50LNx8XL~%b1lpv9URCQ
z#$#~9*~}O`hZkln>s>b<hZW0&^w0i(x5WU6FvlzhQELW-#GEi`4AN*8^l%oP;&~yV
zim(!nBH_b>rcx**3OW=mvjKD}jB@aUj#{u005TnF-sLs|0f3Ua+y;<oD{TOkKwo8}
zlc~^WpcTGUBGff2Y#<p5bD0gGkjPMFR@!JQ_X8*(arJv348hBNsT4>XXoU@=5}@0=
z{5_DmdOwg(TJ1}xtn#H3iK}x!1y`Mq4ld0ppT~f{K=7o}d7kXV5V&seCHCNO_-JYU
hWfu{?_-JWAO1U<M&*1T;l?Lb_6;w1d-05no_&?42*m?i}

literal 0
HcmV?d00001

diff --git a/press/report/report.tex b/press/report/report.tex
index 94635af..c6fcbbd 100644
--- a/press/report/report.tex
+++ b/press/report/report.tex
@@ -5,7 +5,13 @@
 \usepackage{graphicx}
 \usepackage[margin=1.in,centering]{geometry}
 \usepackage{hyperref}
+\usepackage{caption}
+\usepackage[export]{adjustbox}
+\usepackage{float}
 
+\newcommand{\mnras}{MNRAS}
+\newcommand{\apj}{ApJ}
+\newcommand{\aapr}{A\&ARv}
 
 \begin{document}
 
@@ -33,7 +39,7 @@ Power spectral densities and time delays of 19 wavelength bands are recovered as
 The local Type-I Seyfert galaxy NGC 5548, while perhaps the best-studied active galaxy, remains an object of intense interest and study to modern astronomy. Direct observation of active galactic nuclei (AGN) such as that thought to exist at the center of NGC 5548 is rarely possible. The astronomer may infer the properties of AGN from the dynamics of their variable spectra. \cite{2016ApJ...821...56F} published the most complete set of time-dependent light curves yet collected from an active galactic nucleus as part III of STORM, an extensive optical/UV observational campaign carried out on NGC 5548. We now attempt to use frequency-domain analyses to map the reverberation in the observed light curves.
 
 \section{Reverberation Mapping}
-	One model for AGN suggests that a hot accretion disk is incident upon a central super-massive black hole (SMBH). Electromagnetic emission emergent from the gas surrounding the SMBH is reprocessed by the disk, resulting in observed response delays between emission peaks. If the temperature of the disk decreases radially from the SMBH, the observed time delays can be expected to increase with decreasing wavelength. A transfer function encodes the geometry of the system by describing the time-dependent response of each light curve against the others. Recovering the function from the observed light curves is a primary goal of reverberation mapping.
+	One model for AGN suggests that a hot accretion disk is incident upon a central super-massive black hole (SMBH). Electromagnetic emission emergent from the gas surrounding the SMBH is reprocessed by the disk, resulting in observed response delays between emission peaks. If the temperature of the disk decreases radially from the SMBH, the observed time delays can be expected to increase with decreasing wavelength. Furthermore, as the emissions are reprocessed in the disk, they become blurred, so the variability can be expected to decrease with wavelength. A transfer function encodes the geometry of the system by describing the time-dependent response of each light curve against the others. Recovering the function from the observed light curves is a primary goal of reverberation mapping.
 
 	This technique has become a standard for calculating the black hole mass of AGN. It is well-described by \cite{2007MNRAS.380..669C} and \cite{2014A&ARv..22...72U} and many others. It continues to be refined, and may also become a tool to measure the black hole spin of these systems (\cite{2016arXiv160606736K}).
 
@@ -45,7 +51,7 @@ The local Type-I Seyfert galaxy NGC 5548, while perhaps the best-studied active
 
 \subsection{Frequency-domain Analysis}
 
-	The transfer function is related to two light curves as being convoluted against the reference light curve (\ref{time_transfunc}). The convolution theorem 
+	The transfer function $g(t)$ is related to two light curves as being convoluted against the reference light curve (\ref{time_transfunc}). The convolution theorem provides that convolution in the time domain is equivalent to pointwise multiplication in the frequency domain (\ref{freq_transfunc}). The transfer function therefore lends itself well to frequency-domain analyses.
 
 	\begin{equation}
 		\label{time_transfunc}
@@ -57,60 +63,144 @@ The local Type-I Seyfert galaxy NGC 5548, while perhaps the best-studied active
 		Y(\nu) = G(\nu) X(\nu)
 	\end{equation}
 
+    The power spectral density (PSD) provides a measure of the variability of the flux in a band. It is defined as $|X(\nu)|^2 = X^*(\nu)X(\nu)$, where $^*$ denotes the complex conjugate, and can be produced using Fourier transforms. The transfer function is therefore related to the PSD of two light curves:
 
+    \begin{equation}
+        |Y(\nu)|^2 = |G(\nu)|^2 |X(\nu)|^2
+    \end{equation}
 
-	The power spectral density (PSD) for a given emission band can be produced using Fourier transforms, providing a measure of the time-scale of variability in that band. Given two bands, typically a reference or "driving" band and a delayed or "response" band, a cross spectrum can also be constructed. From the complex argument of the cross-correlation function, one can derive the frequency-dependent time delay between those bands; an important step toward constituting the transfer function of a system. Very good explanations of these techniques and the associated mathematics are available from \cite{2014A&ARv..22...72U}.
+    Given two bands, a cross spectrum can also be computed. The cross spectrum is defined as $C(\nu) = X^*(\nu) Y(\nu)$. The argument $\phi$ of the cross spectrum is the phase lag between the two signals. The time delay $\tau$ can therefore be computed from the cross spectrum using
+    
+    \begin{equation}
+        \tau(\nu) = \frac{\phi(\nu)}{2\pi\nu}
+    \end{equation}
 
-	A top-hat function provides a simple model of the impulse response of a delayed light curve. A fast Fourier transform method of this impulse response provides the time delay spectrum as a function of temporal frequency. This simple model provides a guideline for how the computed time delays are expected to be distributed as a function of temporal frequency.
+    Given equation \ref{freq_transfunc}, the cross spectrum can be written as
+        
+    \begin{equation}
+        C(\nu) = X^*(\nu) G(\nu) X(\nu) =  G(\nu) |X(\nu)|^2 
+    \end{equation}
+
+    The time delays are therefore trivially predicted as a function of frequency from the cross spectrum. The frequency dependence of these lags in turn relates directly to the transfer function. Very good explanations of these techniques and the associated mathematics are available from \cite{2014A&ARv..22...72U}.
+
+    \subsection{Tophat Transfer Function}
+
+	A tophat function provides a simple model of the impulse response of a delayed light curve. A fast Fourier transform method of this impulse response provides the time delay spectrum as a function of temporal frequency. This simple model provides a guideline for how the computed time delays are expected to be distributed as a function of frequency. Once the time delays are extracted from the observational datasets, they might be fitted with a tophat function, however, a more complicated function such as a log-Gaussian function is probably more appropriate.
+
+    \begin{figure}
+        \centering
+        \begin{minipage}{.475\textwidth}
+            \centering
+            \includegraphics[width=1\linewidth]{../img/tophat_timedomain.pdf}
+            \captionof{figure}{Tophat functions in the time domain show an average time delay of the reverberating curve and a constant distribution in time over an interval. An area of unity indicates no loss of signal in the response.}
+            \label{fig:th_time}
+        \end{minipage}
+        \hfill
+        \begin{minipage}{.475\textwidth}
+            \centering
+            \includegraphics[width=1\linewidth]{../img/tophat_freqdomain.pdf}
+            \captionof{figure}{The time delays associated with each tophat function. Distinct features related to the average time lag are present (maximum, value of $\nu$ at steepest change), and complicated relationships with higher frequency waves can be noted.}
+            \label{fig:th_freq}
+        \end{minipage}
+    \end{figure}
 
-	(Side-by-side graphic of top-hat impulse response function and FFT of top-hat giving time delays.)
 
 	\subsection{Unevenly-Sampled Data}
 
+    Traditional frequency-domain analyses require data that is evenly sampled. Due primarily to weather, optical reverberation mapping datasets generally contain unevenly sampling. Because of this, until now, optical reverberation mapping has been limited mainly to time-domain analyses, e.g., cross-correlation. While it can handle datasets with significant sampling-variability, cross-correlation is only able to determine the average time lag for a given light curve; however, more information is contained within the light curves than just their average time lag.
 
-	Optical reverberation mapping techniques involve time-domain analyses, such as cross-correlation. Time-domain techniques have limitations: cross-correlation, for instance, provides only the average time delay between two light curves; they also require data that is evenly-sampled across the time domain. X-ray reverberation mapping in particular has developed a body of techniques based on frequency-domain techniques, primarily Fourier analysis; these techniques still require evenly-sampled data and have been enabled by the relatively good data coverage in X-ray bands. They provide the astronomer with more detailed information about the variability and response delay within the system compared to time-domain techniques.
-
-
- 	Some X-ray datasets contain gaps due to orbital mechanics, which motivated the work in \cite{2013ApJ...777...24Z}, where a maximum likelihood method is used to perform Fourier analysis on light curves with gaps. Since its development, this technique has found success among studies of observations captured by low-orbit X-ray telescopes that exceed the telescopes' orbital periods, such as the analysis performed by \cite{2016arXiv160606736K}. Until now, reverberation mapping in the optical bands has been limited to time-domain techniques. Many datasets available for these bands have uneven sampling across the time domain, however, and so do not lend themselves well to time-domain or traditional frequency-domain analyses. The maximum likelihood method is well-suited to extracting useful information from the data available in those datasets.
-
+    Some X-ray datasets contain gaps due to orbital mechanics, which motivated the work by \cite{2013ApJ...777...24Z}, where a maximum likelihood method is used to perform Fourier analysis on light curves with gaps. Since its development, this technique has found success among studies of observations captured by low-orbit X-ray telescopes that exceed the telescopes' orbital periods, such as the analysis performed by \cite{2016arXiv160606736K}. This technique is now being applied to the optical datasets published in STORM III. If successful, it may provide new insight into the reverberations present in the accretion disk and other structures of the nucleus in NGC 5548.
+ 	
 \section{Analysis}
-The 1367\AA$ $ light curve, obtained from observations made with the Hubble Space Telescope, is chosen as the reference curve. The power spectral densities and time delays as a function of temporal frequency are computed for each band in the dataset -- 18 bands not including the reference band.
 
-The light curves analysed here are unevenly distributed along the time axis, which suggests that the maximum likelihood method developed by \cite{2013ApJ...777...24Z} is a reasonable candidate for producing the PSD and time delays in the frequency domain. The latest version (CHECK THIS) of the C++ program psdlag associated with that work is used to directly produce the PSD and cross spectra. The time delay spectrum is produced from the cross spectrum by dividing it by $2 \pi f$, with $f$ the mean frequency for a given bin.
+\cite{2016ApJ...821...56F} published the best dynamic data yet collected from NGC 5548 over a 260-day period, for 19 bands throughout the optical and into the UV domains. These data were collected from a variety of observatories, including both space and ground-based telescopes, and thus have significantly uneven and variable sampling rates. The 1367\AA$ $ light curve, obtained from observations made with the Hubble Space Telescope, is chosen as the reference curve. The power spectral densities and time delays as a function of temporal frequency are computed for each band in the dataset -- 18 bands not including the reference band.
 
-	\subsection{Dataset}
-	\cite{2016ApJ...821...56F} published the best dynamic data yet collected from NGC 5548 over a 200-day (CHECK THIS) period, for 19 bands throughout the optical and into the UV spectra. These data were collected from a variety of observatories, including both space and ground-based telescopes, and thus have significantly variable sampling rates.
+In STORM III, a reverberation mapping analysis is performed using cross-correlation to find the average time delay for each wavelength. These results are compared to a number of possible models, however, the average time lag leaves a lot of uncertainty when trying to constrain an appropriate model. More information is contained in the light curves, and a frequency-domain analysis should provide better constraints.
 
-	(Include picture of Fausnaugh data here)
+\begin{figure}
+    \centering
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\linewidth]{../img/lightcurves.pdf}
+        \captionof{figure}{Data published by \cite{2016ApJ...821...56F} in STORM III shows highly-variable light curves with significantly uneven sampling.}
+        \label{fig:lightcurves}
+    \end{minipage}
+    \hfill
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\textwidth]{../img/TCCF_fausnaugh.pdf}
+        \captionof{figure}{Cross-correlation analysis of lightcurves presented by \cite{2016ApJ...821...56F} compared average time lags to AGN structure models.}
+        \label{fig:cc_analysis}
+    \end{minipage}
+\end{figure}
 
-	\subsection{Error Analysis}
-	For the presented set of resultant data, the error estimates are extracted from the covariance matrix. This method assumes that the errors between frequency bins are not correlated, so these values only represent a lower limit of the true variability. Scanning the likelihood function can provide better error estimates at the cost of computation time, as can running Monte Carlo simulations. All of these methods are built into the psdlag program provided by \cite{2013ApJ...777...24Z}, however, some issues have prevented proper error analysis using the latter two methods. This is discussed in more detail in section \ref{results}.
 
-\section{Results}
-\label{results}
+The light curves analysed here are unevenly distributed along the time axis, which suggests that the maximum likelihood method developed by \cite{2013ApJ...777...24Z} is a reasonable candidate for analysing these data by producing the PSD and time delays. The latest version as of July, 2016, of the C++ program "psdlag" associated with that work is used to directly produce the PSD and cross spectra. The time delay spectrum is produced from the cross spectrum by dividing it by $2 \pi f$, with $f$ the mean frequency for a given bin.
 
-An atlas of the power spectral densities as functions of temporal frequency for all 18 delayed bands is provided in this section. One is also provided of the time delay spectra for each band. The reference band PSD is also provided separately. Errors presented in these atlases are obtained from the covariance matrix.
+\begin{figure}
+    \centering
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\linewidth]{../img/PSD_1367Å_{σ∊CM}.pdf}
+        \captionof{figure}{Power spectral density for 1367\AA, the chosen reference band.}
+        \label{fig:psd_1367}
+    \end{minipage}
+    \hfill
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\linewidth]{../img/timelag_1367Å_≺_7647Å_{σ∊CM}.pdf}
+        \captionof{figure}{The time delay computed from the cross spectrum of 7647\AA and 1367\AA.}
+        \label{fig:timelag_7647}
+    \end{minipage}
+\end{figure}
 
-(Atlas of PSD)
 
-(Atlas of Time Delays)
+	\subsection{Errors}
+	For the presented set of resultant data, the error estimates are extracted from the covariance matrix. This method assumes that the errors between frequency bins are not correlated, so these values only represent a lower limit of the true variability. Scanning the likelihood function can provide better error estimates at the cost of computation time, as can running Monte Carlo simulations. All of these methods are built into the "psdlag" program provided by \cite{2013ApJ...777...24Z}.
 
-	\subsection{Dubious Error Computations}
-	The errors obtained from the covariance matrix are only a lower estimate of the true error. An error analysis by scanning the likelihood function was attempted, but dubious values led to their exclusion from these results. In the case of
+    An error analysis by scanning the likelihood function was attempted, but this analysis was plagued by computational issues that prevented some models from properly optimizing. Monte Carlo simulations were also attempted as a way of estimating the variability of the resultant values. Many errors obtained from this method were larger than the expected accurate values. Therefore, this analysis was also excluded. The ultimate goal is to use one of these more accurate methods, and the errors currently reported should be considered temporary and unreliable.
 
-	(Example of bad LF error)
+    \subsection{Results}
+    \label{results}
 
-	Monte Carlo simulations were also attempted as a way of estimating the variability of the resultant values. Many errors obtained from this method were much larger than the expected accurate values. Therefore, this analysis was also excluded.
+    In figure \ref{psd_atlas}, an atlas of the power spectral densities for all 18 reverberated bands is provided. Figure \ref{timelag_atlas} provides the time delay spectra for each band. The reference band PSD is provided separately. Producing the time delay map is a significant step toward recovering the transfer function.
 
-	(Example of bad MC error)
+    \begin{figure}
+        \centering
+        \includegraphics[width=.9\textwidth]{../img/psd_atlas.pdf}
+        \caption{Power spectral densities for all observed light curves, excluding the reference curve. The decrease in variability with increasing wavelength is in agreement with that prediction.}
+        \label{psd_atlas}
+    \end{figure}
 
+    \begin{figure}
+        \centering
+        \includegraphics[width=.9\textwidth]{../img/timelag_atlas.pdf}
+        \caption{Time delays for all observed light curves relative to 1367\AA. The delay increases with wavelength, as predicted by the accretion disk model. These functions can ultimately be used to reconstruct the transfer function.}
+        \label{timelag_atlas}
+    \end{figure}
 
 \section{Discussion}
-Frequency-dependent power spectral densities confirm time-dependent variability in the emission strengths for each band. This behaviour is expected for any active galactic nucleus and has been long-confirmed in NGC 5548, so it comes as no surprise to find those results here.
+
+Frequency-dependent power spectral densities confirm time-dependent variability in the emission strengths for each band. This  This behaviour is expected for any active galactic nucleus and has been long-confirmed in NGC 5548, so it comes as no surprise to find those results here.
 
 Analysis of the top-hat impulse response model predicted frequency-dependent time delays, which have been recovered from the light curves in this analysis. Furthermore, the distribution of time delays indicates a wavelength-dependent nature. This warrants further study and analysis.
 
-(Maybe a graph comparing the top-hat time delays to one band's time delays.)
+\begin{figure}
+    \centering
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\linewidth]{../img/tophat_freqdomain.pdf}
+        \captionof{figure}{Time delays modeled from tophat impulse responses.}
+        \label{fig:top_freq_comp}
+    \end{minipage}
+    \hfill
+    \begin{minipage}{.475\textwidth}
+        \centering
+        \includegraphics[width=1\linewidth]{../img/timelag_1367Å_≺_7647Å_{σ∊CM}.pdf}
+        \captionof{figure}{The time delay computed from the cross spectrum of 7647\AA and 1367\AA.}
+        \label{fig:timelag_comp}
+    \end{minipage}
+\end{figure}
 
 The analyses performed on these data have elucidated clear trends in the PSD and time delays. With reverberation mapping, the goal is to recover the transfer function, which encodes the geometry of the system. Recovering the time delays is a significant step toward that goal.
 
@@ -120,48 +210,4 @@ The analyses performed on these data have elucidated clear trends in the PSD and
 
 
 
-
-
-Consider two lightcurves $x(t)$ and $y(t)$, where $x(t)$ is the driving lightcurve and $y(t)$ is the reprocessed lightcurve.  If they are related by a linear impulse response, $g(\tau)$, then: 
-
-
-
-
-
-
-
-
-
-So, $y(t)$ is a delayed and blurred version of $x(t)$, with the amount of delay and blurring encoded in $g(\tau)$.
-
-The power spectral density (PSD) of $x(t)$ is calculated from the Fourier transform of $x(t)$, which we denote $X(\nu)$.  The PSD is $|X(\nu)|^2 = X^*(\nu)X(\nu)$, where the $^*$ denotes the complex conjugate.  From the convolution theorem of Fourier transforms we can write:
-
-
-
-This means it is easy to relate the PSD of the reprocessed lightcurve to the PSD of the driving lightcurve and the impulse response function:
-
-\begin{equation}
-|Y(\nu)|^2 = |G(\nu)|^2 |X(\nu)|^2
-\end{equation}
-
-The cross spectrum is defined as
-\begin{equation}
-C(\nu) = X^*(\nu) Y(\nu)
-\end{equation}
-the phase, $\phi$, of which gives the phase lag between X and Y at each Fourier frequency, $\nu$.  This can be converted to a time lag through: 
-\begin{equation}
-\tau(\nu) = \frac{\phi(\nu)}{2\pi\nu}
-\end{equation}
-Since $Y(\nu) = G(\nu) X(\nu)$,  the cross spectrum can be written as:
-\begin{equation}
-C(\nu) = X^*(\nu) G(\nu) X(\nu) =  G(\nu) |X(\nu)|^2 
-\end{equation}
-thus, for a given impulse response function, one can trivially predict the time lags as a function of frequency, $\tau(\nu)$, by calculating the phase of $G(\nu)$, and the frequency dependence of the lags directly relates to the shape of the response function.
-
-
-
-
-
-
-
 \end{document}
\ No newline at end of file