From 647d192435d36d1155adc8f1a3e56197323bd036 Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Wed, 5 Aug 2020 02:54:30 -0400 Subject: [PATCH] Initial commit --- Gallavotti_Jauslin_2020.tex | 666 +++++++++++++++++++++ Makefile | 46 ++ README | 31 + fig1.pdf | Bin 0 -> 5399 bytes fig2.pdf | Bin 0 -> 98309 bytes figs/action-angle.fig/Makefile | 33 + figs/action-angle.fig/action-angle.gnuplot | 25 + figs/action-angle.fig/etas.dat | 49 ++ figs/ellipse.fig/Makefile | 24 + figs/ellipse.fig/ellipse.tikz.tex | 86 +++ figs/trajectory.fig/Makefile | 30 + figs/trajectory.fig/trajectory.tikz.tex | 26 + 12 files changed, 1016 insertions(+) create mode 100644 Gallavotti_Jauslin_2020.tex create mode 100644 Makefile create mode 100644 README create mode 100644 fig1.pdf create mode 100644 fig2.pdf create mode 100644 figs/action-angle.fig/Makefile create mode 100644 figs/action-angle.fig/action-angle.gnuplot create mode 100644 figs/action-angle.fig/etas.dat create mode 100644 figs/ellipse.fig/Makefile create mode 100644 figs/ellipse.fig/ellipse.tikz.tex create mode 100644 figs/trajectory.fig/Makefile create mode 100644 figs/trajectory.fig/trajectory.tikz.tex diff --git a/Gallavotti_Jauslin_2020.tex b/Gallavotti_Jauslin_2020.tex new file mode 100644 index 0000000..3d00ce2 --- /dev/null +++ b/Gallavotti_Jauslin_2020.tex @@ -0,0 +1,666 @@ +\documentclass[10pt]{article} + +\usepackage{color} +\usepackage{graphicx} +\usepackage{amsfonts} +\usepackage[hidelinks]{hyperref} +\usepackage{natbib} + +\def\Eq#1{\label{#1}} + +% colors +\definecolor{iblue}{RGB}{65,105,225} +\definecolor{ired}{RGB}{220,20,60} +\definecolor{igreen}{RGB}{50,205,50} +\definecolor{ipurple}{RGB}{75,0,130} +\definecolor{iochre}{RGB}{218,165,32} +\definecolor{iteal}{RGB}{51,204,204} +\definecolor{imauve}{RGB}{204,51,153} +\definecolor{RED}{RGB}{255,0,0} + +\def\alertv#1{{\color{green}#1}} +\def\alertm#1{{\color{magenta}#1}} +\def\alertb#1{{\color{blue}#1}} +\def\alertr#1{{\color{red}#1}} +\def\alertr#1{{\color{RED}#1}} +\def\alertn#1{{\color{black}#1}} + + +%% symbols +\let\a=\alpha \let\b=\beta \let\g=\gamma \let\d=\delta \let\e=\varepsilon +\let\z=\zeta \let\h=\eta \let\th=\vartheta \let\k=\kappa \let\l=\lambda +\let\m=\mu \let\n=\nu \let\x=\xi \let\p=\pi \let\r=\rho +\let\s=\sigma \let\t=\tau \let\f=\varphi \let\ph=\varphi\let\ch=\chi +\let\ch=\chi \let\ps=\psi \let\y=\upsilon \let\o=\omega \let\si=\varsigma +\let\G=\Gamma \let\D=\Delta \let\Th=\Theta \let\L=\Lambda \let\X=\Xi +\let\P=\Pi \let\Si=\Sigma \let\F=\Phi \let\Ps=\Psi +\let\O=\Omega \let\Y=\Upsilon +\def\V#1{{\bf#1}}\def\lhs{{\it l.h.s.}\ }\def\rhs{{\it r.h.s.}\ } +\def\*{\vskip 3mm} +\def\0{\noindent} +\def\be{\begin{equation}} +\def\ee{\end{equation}} +\def\bea{\begin{eqnarray}} +\def\eea{\end{eqnarray}} +%\renewcommand{\theequation}{\arabic{section}.\arabic{equation}} +\def\nn{\nonumber} +\def\xx{{\V x}}\def\pp{{\V p}}\def\kk{{\V k}} +\def\AA{{\mathcal A}} +\def\BB{{\cal B}}\def\CC{{\cal C}}\def\DD{{\mathcal D}} +\def\EE{{\cal E}}\def\HH{{\cal H}}\def\KK{{\cal K}}\def\LL{{\mathcal L}} +\def\NN{{\cal N}}\def\FF{{\cal F}}\def\PP{{\mathcal P}} +\def\QQ{{\mathcal Q}}\def\RR{{\cal R}}\def\TT{{\mathcal T}} +\let\dpr=\partial\let\fra=\frac +\def\ie{{\it i.e.}\ } +\def\eg{{\it e.g.}\ } +\def\equ#1{(\ref{#1})} +\def\lis#1{\overline{#1}} +\def\defi{{\buildrel def\over=}} +\def\Media#1{{\Big\langle\,#1\,\Big\rangle}} +\def\media#1{{\Blangle\,#1\,\Brangle}} +\def\bra#1{{\Blangle#1\Bvert}}\def\ket#1{{\Bvert#1\Brangle}} +\def\braket#1#2{\Blangle#1\Bvert#2\Brangle} +\def\otto{\,{\kern-1.truept\leftarrow\kern-5.truept\to\kern-1.truept}\,} +\def\wt#1{\widetilde{#1}} +\def\wh#1{\widehat{#1}} +\def\tende#1{\,\vtop{\ialign{##\crcr\rightarrowfill\crcr + \noalign{\kern-1pt\nointerlineskip} \hskip3.pt${\scriptstyle + #1}$\hskip3.pt\crcr}}\,} +\def\ie{{\it i.e.\ }}\def\etc{{\it etc.\ }} +\def\FF{{\mathcal F}} +\def\tto{\Rightarrow} +\def\ap{{\it a priori}} +\def\Ba {{\mbox{\boldmath$ \alpha$}}} +\def\Bb {{\mbox{\boldmath$ \beta$}}} +\def\Bg {{\mbox{\boldmath$ \gamma$}}} +\def\Bd {{\mbox{\boldmath$ \delta$}}} +\def\Be {{\mbox{\boldmath$ \varepsilon$}}} +\def\Bee {{\mbox{\boldmath$ \epsilon$}}} +\def\Bz {{\mbox{\boldmath$\zeta$}}} +\def\Bh {{\mbox{\boldmath$ \eta$}}} +\def\Bthh {{\mbox{\boldmath$ \theta$}}} +\def\Bth {{\mbox{\boldmath$ \vartheta$}}} +\def\Bi {{\mbox{\boldmath$ \iota$}}} +\def\Bk {{\mbox{\boldmath$ \kappa$}}} +\def\Bl {{\mbox{\boldmath$ \lambda$}}} +\def\Bm {{\mbox{\boldmath$ \mu$}}} +\def\Bn {{\mbox{\boldmath$ \nu$}}} +\def\Bx {{\mbox{\boldmath$ \xi$}}} +\def\Bom {{\mbox{\boldmath$ \omega$}}} +\def\Bp {{\mbox{\boldmath$ \pi$}}} +\def\Br {{\mbox{\boldmath$ \rho$}}} +\def\Bro {{\mbox{\boldmath$ \varrho$}}} +\def\Bs {{\mbox{\boldmath$ \sigma$}}} +\def\Bsi {{\mbox{\boldmath$ \varsigma$}}} +\def\Bt {{\mbox{\boldmath$ \tau$}}} +\def\Bu {{\mbox{\boldmath$ \upsilon$}}} +\def\Bf {{\mbox{\boldmath$ \phi$}}} +\def\Bff {{\mbox{\boldmath$ \varphi$}}} +\def\Bch {{\mbox{\boldmath$ \chi$}}} +\def\Bps {{\mbox{\boldmath$ \psi$}}} +\def\Bo {{\mbox{\boldmath$ \omega$}}} +\def\Bome {{\mbox{\boldmath$ \varomega$}}} +\def\BG {{\mbox{\boldmath$ \Gamma$}}} +\def\BD {{\mbox{\boldmath$ \Delta$}}} +\def\BTh {{\mbox{\boldmath$ \Theta$}}} +\def\BL {{\mbox{\boldmath$ \Lambda$}}} +\def\BX {{\mbox{\boldmath$ \Xi$}}} +\def\BP {{\mbox{\boldmath$ \Pi$}}} +\def\BS {{\mbox{\boldmath$ \Sigma$}}} +\def\BU {{\mbox{\boldmath$ \Upsilon$}}} +\def\BF {{\mbox{\boldmath$ \Phi$}}} +\def\BPs {{\mbox{\boldmath$ \Psi$}}} +\def\BO {{\mbox{\boldmath$ \Omega$}}} +\def\BDpr {{\mbox{\boldmath$ \partial$}}} +\def\Bstl {{\mbox{\boldmath$ *$}}} +\def\Brangle {{\mbox{\boldmath$ \rangle$}}} +\def\Blangle {{\mbox{\boldmath$ \langle$}}} +\def\Bvert{{\mbox{\boldmath$|$}}} +\def\Bell {{\mbox{\boldmath$\ell$}}} +\let\up\uparrow +\let\down\downarrow +\def\({\left(} +\def\){\right)} +\let\mc\mathcal +\let\mrm\mathrm + + +\def\iniz{\setcounter{equation}{0}} +\renewcommand{\theequation}{\arabic{section}.\arabic{equation}} +\def\inizA{\setcounter{equation}{0} +\renewcommand{\theequation}{\Alph{section}.\arabic{equation}} +} +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +\newdimen\xshift \newdimen\xwidth \newdimen\yshift \newdimen\ywidth% +\def\ins#1#2#3{\vbox to0pt{\kern-#2pt\hbox{\kern#1pt #3}\vss}\nointerlineskip} + +\def\eqfig#1#2#3#4#5{ +\par\xwidth=#1pt \xshift=\hsize \advance\xshift +by-\xwidth \divide\xshift by 2 +\yshift=#2pt \divide\yshift by 2% +{\hglue\xshift \vbox to #2pt{\vfil +#3 \special{psfile=#4.eps} +}\hfill\raise\yshift\hbox{#5}}} +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% +\def\Eqfig#1#2#3#4#5#6{ +\par\xwidth=#1pt \xshift=\hsize \advance\xshift +by-\xwidth \divide\xshift by 2 +\yshift=#2pt \divide\yshift by 2% +{\hglue\xshift \vbox to #2pt{\vfil +#3 \special{psfile=#4.eps}\kern200pt\special{psfile=#5.eps} +}\hfill\raise\yshift\hbox{#6}}} +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% + +\def\eqalign#1{\null\,\vcenter{\openup\jot + \ialign{\strut\hfil$\displaystyle{##}$&$\displaystyle{{}##}$\hfil + \crcr#1\crcr}}\,} + +\def\qedsymbol{$\square$} +\def\qed{\penalty10000\hfill\penalty10000\qedsymbol} + +\date{} +\author{\alertb{Giovanni Gallavotti${}^1$ and +Ian Jauslin${}^2$}} + +\title{\alertr{\bf A Theorem on Ellipses, an Integrable System and a Theorem + of Boltzmann} +} + +\begin{document} +\maketitle +\kern-8mm +\centerline{${}^1$ INFN-Roma1 \& Universit\`a ``La Sapienza'', email: giovanni.gallavotti@roma1.infn.it} +\centerline{${}^2$ Department of Physics, Princeton University, email: ijauslin@princeton.edu} + + +%\kern-1cm +\begin{abstract} +% +We study a mechanical system that was considered by Boltzmann in 1868 in the +context of the derivation of the canonical and microcanonical ensembles. This +system was introduced as an example of ergodic dynamics, which was central to +Boltzmann's derivation. It consists of a single particle in two dimensions, +which is subjected to a gravitational attraction to a fixed center. In +addition, an infinite plane is fixed at some finite distance from the center, +which acts as a hard wall on which the particle collides elastically. Finally, +an extra centrifugal force is added. We will show that, in the absence of this +extra centrifugal force, there are two independent integrals of motion. +Therefore the extra centrifugal force is necessary for Boltzmann's claim of +ergodicity to hold. +% +\end{abstract} +\* +% +\0Keywords: {\small Ergodicity, Chaotic hypothesis, Gibbs +distributions, Boltzmann, Integrable systems} +\* +%\end{document} + +In 1868, \cite{Bo868a} laid the foundations for our modern understanding of the +behavior of many-particle systems by introducing the ``microcanonical +ensemble'' (for more details on this history, see \cite{Ga016}). The principal +idea behind this ensemble is that one can achieve a good understanding of +many-particle systems by focusing not on the dynamics of each individual +particle, but on the statistical properties of the whole. More precisely, the +state of the system becomes a random variable, chosen according to a +probability distribution on phase space, which came to be called the +``microcanonical ensemble''. An important assumption that was made implicitly +by Boltzmann is that the dynamics of the system be ergodic. In this case, +time-averages of the dynamics can be rewritten as averages over phase space, +and the qualitative properties of the dynamics can be formulated as statistical +properties of the microcanonical ensemble. + +To support this assumption, Boltzmann presented a mechanical system that +very same year (\cite{Bo868b}) as an example of an ergodic system. This +system consists of a particle in two dimensions that is attracted to a +fixed center via a gravitational potential $-\frac{\alpha}{2r}$. In +addition, he added an extra centrifugal potential $\frac g{2r^2}$. As was +known since at least the times of Kepler, this system is subjected to a +central force, and is therefore integrable. In order to break the +integrability, Boltzmann added an extra ingredient: a rigid infinite planar +wall, located a finite distance away from the center (see figure +\ref{trajectory}). Whenever the particle hits the wall, it undergoes an elastic +collision and is reflected back. Boltzmann's argument was, roughly, that in +the absence of the wall, the dynamics is quasi-periodic, so the particle +should intersect the plane of the wall at points which should fill up a +segment of the wall densely as the dynamics evolves, and concluded that the +region of phase space in which the energy is constant must also be filled +densely. As we will show, this is not the whole story; following a +conjectured integrability for $g=0$, \cite[p.150]{Ga013b}, +and first tests +%(by (GG) proposing a relation with KAM theory for $g>0$ and by (IJ) proposing chaotic motions at large $g$) +in +\cite[p.225--228]{Ga016}, we have found that, in the absence of the +centrifugal term $g=0$, the dynamics (which has two degrees of freedom) still admits +two constants of motion even in presence of the hard wall. This +suggests that, if a suitable KAM analysis could be carried out, the system +would not be ergodic for small values of $g$. + +\begin{figure}[ht] +\hfil\includegraphics[width=6cm]{trajectory.pdf} + +\caption{A trajectory. The large dot is the attraction center $O$, and the line + is the hard wall $\LL$. In between collisions, the trajectories are ellipses. + The ellipses are drawn in full, but the part that is not covered by the + particle is dashed. + } +\label{trajectory} +\end{figure} + + +\setcounter{section}{0} +\section{Definition of the model and main result} +\label{sec1} +\iniz + + +Let us now specify the model formally, and state our main result more +precisely. We fix the gravitational center to the origin of the $x,y$-plane and +let $\LL$ be the line $y=h$. The Hamiltonian for the system in between +collisions is +% +\be H=\frac{p_x^2+p_y^2}2 -\frac{\a}{2r}+\frac{g}{2r^2} +\Eq{e1.1}\ee +% +where $\a>0,g\ge0,r=\sqrt{x^2+y^2}$ and the particle moves following +Hamilton's equations as long as it stays away from the obstacle $\LL$. When +an encounter with $\LL$ occurs the particle is reflected elastically and +continues on. + +\cite{Bo868b}, considered this system on the hyper-surface $A={\V +p}^2-\frac\a r+\frac{g}{r^2}$. The intersection of this hyper-surface with +$y=h$ is the region $\FF_A$ enclosed within the curves +% +\be \pm\sqrt{(A-\frac{g}{x^2+h^2} +\frac{\a}{\sqrt{x^2+h^2}})},\qquad +x_{min}0$. +\bigskip + +From now on, unless it is explicitly stated otherwise, we will assume that +$g=0$. +\bigskip + +In this case, the motion between collisions takes place at constant +energy $\frac12A$ and constant angular momentum $a$, and traces out an +ellipse. One of the foci of the ellipse is located at the origin, and we +will denote the angle that the aphelion of the ellipse makes with the +$x$-axis by $\theta_0$. Thus, the ellipse is entirely determined by the triplet +$(A,a,\theta_0)$. When a collision occurs, $A$ remains unchanged, but $a$ and +$\theta_0$ change discontinuously to values $(a',\theta_0')=\FF(a,\theta_0)$, and thus +the Kepler ellipse of the trajectory changes. In addition, the semi-major +axis $a_M$ of the ellipse is also fixed to $a_M=-\frac\a{2A}$ (Kepler's +law): so the successive ellipses have the same semi-major axis, while the +eccentricity varies because at each collision the angular momentum changes: +$e^2=1+ \frac{4 A a^2}{\a^2}$. Thus, the motion will take place on arcs of +various ellipses $\EE$, which all share the same focus and the same semi-major +axis, but whose angle and eccentricity changes at each collision. + +Our main result is that the (canonical) map $(a',\f'_0)=\FF(a,\theta_0)$, which +maps the angular momentum and angle of the aphelion before a collision to their +values after the collision, admits a constant of motion. This follows from the +following geometric lemma about ellipses. +\bigskip + +\0{\bf Lemma 1:} {\it Given an ellipse $\mathcal E$ with a focus at $O$ that +intersects $\LL$ at a point $P$. Let $Q$ denote the orthogonal projection of +$O$ onto $\LL$ (see figure \ref{fig1}). The distance $R_0$ between $Q$ and +the center of $\mathcal E$ depends solely on the semi-major axis $a_M$, the +distance $r$ from $O$ to $P$, and $\cos(2\lambda)$ where $\lambda$ is the angle +between the tangent of the ellipse at $P$ and $\LL$ (to define the direction of +the tangent, we parametrize the ellipse in the counter-clockwise direction): +\be +R_0=\sqrt{\frac14 r^2+\frac14(2a_M-r)^2 +\frac12 r (2a_M-r)\cos(2\l)} +. +\Eq{e1.3}\ee +} +\* + +\begin{figure}[ht] +\hfil\includegraphics[width=100pt]{fig1.pdf} + +\caption{The attractive center is $O$, hence it is the focus of the + ellipse in absence of centrifugal force $g=0$. $Q$ is the projection of + $O$ on the line $\LL$ and $P$ is a collision point. The arrow + represents the velocity of the particle after the collision.} +\label{fig1} +\end{figure} + +\underline{Proof}: We switch to polar coordinates +$p=(r\cos\f,r\sin\f)$. + +Let $O'$ denote the other focus of the ellipse, and $C$ denote its center. The +first step is to compute the vector $\protect\overrightarrow{O'P}$, which in polar +coordinates is +\be +\protect\overrightarrow{O'P}=((2a_M-r)\cos\f',(2a_M-r)\sin\f')\Eq{e1.4}\ee +% +Let $\psi:=\pi+\f-\lambda$ denote the angle between the tangent of the +ellipse at $P$ and the vector $\protect\overrightarrow{PO}$ (see figure \ref{ellipse}), +and $\psi':=\pi+\f'-\lambda$ denote the angle between the tangent of the +ellipse at $P$ and the vector $\protect\overrightarrow{PO'}$. + +\begin{figure}[ht] + \hfil\includegraphics[width=8cm]{ellipse.pdf} + +\caption{An ellipse with foci $O$ and $O'$ and center $C$. The thick line is +$\LL$, which intersects the ellipse at $P$, and $Q$ is the projection of $O$ +onto $\LL$. The dashed line is the tangent at $P$. $\lambda$ is the angle +between $\LL$ and the tangent, $\f$ is the polar coordinate, $\f'$ is the angle +between $\LL$ and $\protect\overrightarrow{O'P}$. $\psi$ is the angle between the +tangent and $\protect\overrightarrow{PO}$, which is equal to the angle between the +tangent and $\protect\overrightarrow{PO'}$. $R_0$ is the distance between $Q$ and $C$.} +\label{ellipse} +\end{figure} + +By the focus-to-focus +reflection property of ellipses, we have $\psi'=\pi-\psi$. Thus +$\f'=2\lambda-\pi-\f$ and we find; + + +\begin{figure}[ht] + \hfil\includegraphics[width=8cm]{fig2.pdf} + +\caption{Two ellipses, before and after a collision. The collision line $\LL$ +is the line at $y=1$, $P$ is the collision point; $Q$ is the projection of $O$ +onto $\LL$; the two ellipses $\EE$ and $\EE'$ have a common focus $O$, and +$O,O'$ are the foci of $\EE$, whereas $O,O''$ are the foci of $\EE'$; $C$ and +$C''$ are the centers of $\EE$ and $\EE'$ respectively; the ellipses are drawn +completely although the trajectory is restricted to the parts above $y=h=1$. +The distance from $C''$ to $Q$ is the same as that from $C$ to $Q$. The upper +ellipse $\EE$ contains the trajectory that starts at the collision point $P$ +following the other ellipse $\EE'$ which has undergone reflection.} +\label{fig2} +\end{figure} + +\be + R_0^2=|Q-C|^2 + = + \frac14\left(r^2+(2a_M-r)^2+2r(2a_M-r)\cos(2\lambda)\right) + . +\Eq{e1.5}\ee +See figures \ref{ellipse} and \ref{fig2}.\qed + +\* +\0{\bf Theorem 1}: {\it The quantity +% +\be R= a^2+h\a e \sin\theta_0\Eq{e1.6} +\equiv\frac\alpha{2a_M}(h^2+a_M^2-R_0^2) +\Eq{e1.7}\ee +% +where $e$ is the eccentricity $e=\sqrt{1+\frac{4 A a^2}{\a^2}}$, is a constant +of motion.} +\* + +\underline{Proof}: +During a collision, the value of $\l$ changes from $\l$ to $\p-\l$, while +$r$ and $a_M$ stay the same. By lemma 1, this implies that the distance $R_0$ +between $Q$ and the center of the ellipse is preserved during a collision. +Furthermore, the position of the center $C$ of the ellipse is given by +$C=a_Me(\cos\theta_0,\sin\theta_0)$ +so +\be + R_0^2=|Q-C|^2=a_M^2e^2-2a_Meh\sin\theta_0+h^2.\Eq{e1.8} +\ee +Furthermore, the angular momentum is equal to +$a^2=\frac12a_M\alpha(1-e^2)$ +so +\be + -R_0^2+h^2+a_M^2 + = + \frac{2a_M}\alpha(a^2+e\alpha h\sin\theta_0)\Eq{e1.9} +\ee +is a conserved quantity. \qed +\bigskip + +\0{\bf Remark:} Some useful inequalities are +% +\be +\eqalign{ + &r_{max}<{2}{a_M}; \ x_{max}=\sqrt{r_{max}^2-h^2};\ + R_0^2\in ((a_M-r)^2,a_M^2);\cr + &\frac{\a h^2}{2a_M}\,<\,R\,<\, + (1+\frac{a_M^2}{h^2}-\Big(\frac{a_M}{h}- + \frac{r}h\Big)^2)\frac{\a h^2}{2a_M}\cr} +\Eq{e1.10}\ee +% +hence in the plane $(x,\l)$ the rectangle $(-x_{max},x_{max})\times(0,\p)$ +(recall that $x_{max}$ is the largest $x$ accessible at energy $\frac12A$) +is the surface of energy $\frac12A$ and the trajectories are the curves of +constant $R$ inside this rectangle. + +\def\SEC{Conjectures on action angle variables} +\section{\SEC} +\label{sec2} +\iniz + +In the previous section, we exhibited a constant of motion, which, along +with the conservation of energy, brings the number of independent conserved +quantities to two. In a continuous Hamiltonian system, this would imply the +existence of action-angle variables, which are canonically conjugate to the +position and momentum of the particle, in terms of which the dynamics +reduces to a linear evolution on a torus. In this case, the collision +with the wall introduces some discreteness into the problem, and the +existence of the action angle variables is not guaranteed by standard +theorems. Indeed, in the presence of the collisions, we no longer have a +Hamiltonian system, but rather a discrete symplectic map (or a +non-differentiable Hamiltonian), which describes the change in the state of +the particle during a collision. In this section, we present some +conjectures pertaining to the existence of action angle variables for this +problem. +\bigskip + +The first step is to change to variables which are action-angle variables for +the motion in between collision. We choose the {\it Delaunay} variables, whose +angles are the argument of the aphelion $\theta_0$ defined above, the {\it mean +anomaly} $M$, and whose actions are the angular momentum $a$, and another +momentum usually denoted by $L$ and related to the semi-major axis $a_M$ and +to the energy $E=\frac12 A$: +\be L:=-\sqrt{\frac{\alpha}2a_M},\quad a_M:=-\frac\alpha{2A} +,\quad +A:=p^2+\frac{a^2}{r^2}-\frac\alpha{r}\equiv-\frac{\alpha^2}{4L^2} +\Eq{e2.1}\ee +It is well known that this change of variables is canonical. In between +collisions, the dynamics of the particle in the variables +$(M,\theta_0;L,a)$ is, simply, +\be + \dot M=\frac{\alpha^2}{4L^3} + ,\quad + \dot\theta_0=0 + ,\quad + \dot L=0 + ,\quad + \dot a=0 + .\Eq{e2.2} +\ee +These variables are thus action-angle variables in between collisions, but when +a collision occurs, $\theta_0$ and $a$ will change. +\bigskip + +The following conjecture states that there exists an action-angle variable +during the collisions. +\bigskip + +\0{\bf Conjecture 1:} {\it There exists a variable $\gamma$ and an integer + $k$ such that, every $k$ collisions, the change in $\gamma$ is + \begin{equation} + \gamma'=\gamma+\o(L,R)\Eq{e2.3} + \end{equation} + in which case $\gamma$ is an angle that rotates on a circle of radius depending + on $L,R$. The function $\o(L,R)$ has a non zero derivative with respect to + $R$ at constant $L$, \ie the motion on the energy surface is quasi periodic + and anisochronous.} +\* + +We will now sketch a construction of this variable $\gamma$, which we obtain +using a generating function $F(L,R,M,\theta_0)$. +\bigskip + +First of all, by theorem 1, the angular momentum $a(\theta_0)$ is a solution of +\begin{equation} + a^2=R-h\a\sin\theta_0\sqrt{1-\fra{a^2}{L^2}}\Eq{e2.4} +\end{equation} +that is, if $\e=\pm$, +% +\be +a^2=R-\frac{h^2\alpha^2}{2L^2}\sin^2\theta_0+ +\e\sqrt{\frac{h^4\alpha^2}{4L^4}\sin^4\theta_0+h^2\a^2\sin^2\theta_0- +\frac{R\alpha^2h^2}{L^2}\sin^2\theta_0}\Eq{e2.5}\ee +% +and $a=\h\sqrt{a^2}$, so that there may be four possibilities for the value of +$a$ denoted $a=a_{\e,\h}(\theta_0,R,L)$ with $\e=\pm,\h=\pm$. The choice of the +signs $\e=\pm1$, and $\h$ must be examined carefully. +\bigskip + +We then define the generating function +\be F(L,R,M,\theta_0)=LM+\int_0^{\theta_0} a(L,R,\ps)d\ps +%-\int_0^{2\pi} a(L,R,\ps)d\ps +\Eq{e2.6}\ee +% +which yields the following canonical transformation: + +\be\eqalign{ + \g=&\dpr_R\int^{\theta_0}_0 a_{\e,\h}(L,R,\ps)\,d\ps + %-\dpr_R\int^{2\pi}_0 a_{\e,\h}(L,R,\ps)\,d\ps + \cr + M'=&M+\dpr_L\int^{\theta_0}_0 a_{\e,\h}(L,R,\ps)\,d\ps + %-\dpr_L\int^{2\pi}_0 a_{\e,\h}(L,R,\ps)\,d\ps + \cr +}\Eq{e2.7}\ee +% +It is natural, if Boltzmann's system is integrable (at $g=0$), that the new +variables are its action angle variables and $M',\gamma$ rotate uniformly in spite of +the collisions. +\bigskip + +However, in this case, the signs $\e$ and $\h$ may change from one collision to +the next, complicating the situation. A careful numerical study of the system +has led us to the following conjecture (see figure \ref{action_angle}). +\bigskip + +% +\0{\bf Conjecture 2:} {\it If $R>h\a$ (which is the case in which the + circle, of radius $R_0$, of the centers encloses the focus $O$), when the + motion collides for the $n$-th time, the angular momentum is proportional + to $(-1)^n$, and, thus, $\epsilon=(-1)^n$. The sign $\eta$ is fixed to + $+$. The increment $\Delta_2\gamma$ in $\gamma$ between the $n$-th and + the $n+2$-th collision is independent of $n$. } \* +\bigskip + +\begin{figure} + \hfil\includegraphics[width=8cm]{action-angle.pdf} + \caption{A plot of the increment in $\gamma$ between the $n$-th and the $n+2$-nd collision as a function of $n$. The blue `+' signs correspond to even $n$, and the red `$\times$' to odd $n$. The variation of $\Delta_2\gamma$ is as small as 1 part per million, thus supporting conjecture 2.} + \label{action_angle} +\end{figure} + +\0{\it Remark:} The change of variables over the variables $a,\theta_0$ to +$R,\g$ at fixed $L$ is {\it remarkably} essentially the same as the one +(\ap\ unrelated) to find action-angle variable for the auxiliary +Hamiltonian $R=R(a,\theta_0)$. This might remain true even when $R0$, provided $g$ is sufficiently small. However it +may still have invariant regions of positive volume where the motion is ergodic. + +\* + + +\0{\bf Acknowledgements}: The authors thank G. Felder for giving us the impetus +to write this note up in its current form, and to publish it. I.J. gratefully +acknowledges support from NSF grants 31128155 and 1802170. + +\bibliographystyle{plainnat} +%\bibliographystyle{alpha} +%\bibliographystyle{amsref} +%\bibliographystyle{apsrmp} +%\bibliographystyle{spmpsci} +%\bibliographystyle{annotate} +%\bibliography{0Bib} +\begin{thebibliography}{3} +\providecommand{\natexlab}[1]{#1} +\providecommand{\url}[1]{\texttt{#1}} +\expandafter\ifx\csname urlstyle\endcsname\relax + \providecommand{\doi}[1]{doi: #1}\else + \providecommand{\doi}{doi: \begingroup \urlstyle{rm}\Url}\fi + +\bibitem[Boltzmann(1868a)]{Bo868a} +L.~Boltzmann. +\newblock Studien \"uber das Gleichgewicht der lebendingen Kraft zwischen bewegten materiellen Punkten. +\newblock \emph{Wiener Berichte}, {\bf 58}, 517--560, (49--96), 1868. + +\bibitem[Boltzmann(1868b)]{Bo868b} +L.~Boltzmann. +\newblock L{\"o}sung eines mechanischen problems. +\newblock \emph{Wiener Berichte}, {\bf 58}, (W.A.,\#6):\penalty0 1035--1044, + (97--105), 1868. + +\bibitem[Gallavotti(2014)]{Ga013b} +G.~Gallavotti. +\newblock \emph{Nonequilibrium and irreversibility}. +\newblock Theoretical and Mathematical Physics. Springer-Verlag, 2014. + +\bibitem[Gallavotti(2016)]{Ga016} +G.~Gallavotti. +\newblock Ergodicity: a historical perspective. equilibrium and nonequilibrium. +\newblock \emph{European Physics Journal H}, {\bf 41}, 181--259, 2016. +\newblock \doi{DOI: 10.1140/epjh/e2016-70030-8}. + +\end{thebibliography} + +\*\* +\end{document} + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%################### + +%plot "lambdak-0.25-3.4--0.3333" u 1:2 every 2::1:::2000 +% plot 'file' every {} +% {:{} +% {:{} +% {:{} +% {:{} +% {:}}}}} +%plot "grafk-0.326753-2--0.3333" u 3:4 every 2:1:0:0:5:0 w l +%plot "gammak-0.3-1--0.2" u 1:3 every 2:1:0:0:2:0 w l + +% Syntax: +% plot 'file' every {} +% {:{} +% {:{} +% {:{} +% {:{} +% {:}}}}} +% +% every 3:1:2::1024:1 +% 2 significa che inizia dal #3 diff --git a/Makefile b/Makefile new file mode 100644 index 0000000..dc76955 --- /dev/null +++ b/Makefile @@ -0,0 +1,46 @@ +PROJECTNAME=$(basename $(wildcard *.tex)) +FIGS=$(notdir $(wildcard figs/*.fig)) + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +SYNCTEXS=$(addsuffix .synctex.gz, $(PROJECTNAME)) + +all: $(PROJECTNAME) + +$(PROJECTNAME): $(FIGS) + pdflatex -file-line-error $@.tex + pdflatex -file-line-error $@.tex + pdflatex -synctex=1 $@.tex + +$(PROJECTNAME).aux: $(FIGS) + pdflatex -file-line-error -draftmode $(PROJECTNAME).tex + + +$(SYNCTEXS): $(FIGS) + pdflatex -synctex=1 $(patsubst %.synctex.gz, %.tex, $@) + + +figs: $(FIGS) + +$(FIGS): + make -C figs/$@ + for pdf in $$(find figs/$@/ -name '*.pdf'); do ln -fs "$$pdf" ./ ; done + + +clean-aux: clean-figs-aux + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f $(addsuffix .out, $(PROJECTNAME)) + rm -f $(addsuffix .toc, $(PROJECTNAME)) + +clean-figs: + $(foreach fig,$(addprefix figs/, $(FIGS)), make -C $(fig) clean; ) + rm -f $(notdir $(wildcard figs/*.fig/*.pdf)) + +clean-figs-aux: + $(foreach fig,$(addprefix figs/, $(FIGS)), make -C $(fig) clean-aux; ) + + +clean-tex: + rm -f $(PDFS) $(SYNCTEXS) + +clean: clean-aux clean-tex clean-figs diff --git a/README b/README new file mode 100644 index 0000000..6ad6333 --- /dev/null +++ b/README @@ -0,0 +1,31 @@ +This directory contains the source files to typeset the article, and generate +the figures. This can be accomplished by running + make + +The figures trajectory.pdf and action-angle.pdf use the data from a computation +which has not been released yet. Instead, the relevant data files are included. + + +* Dependencies: + + pdflatex + TeXlive packages: + amsfonts + color + graphics + hyperref + latex + natbib + pgf + standalone + GNU make + gnuplot + +* Files: + + Gallavotti_Jauslin_2020.tex: + main LaTeX file + + figs: + source code for the figures + diff --git a/fig1.pdf b/fig1.pdf new file mode 100644 index 0000000000000000000000000000000000000000..928aad8cf791d8e0370aa52cb6ded6aa57a495d6 GIT binary patch literal 5399 zcmb_g3s@6Z7Vg@97N=OGVq3MX1Bwx8n3=p2NH7EfQ35eR#Zt7BWCDR?CT1o)R7%xq zYei9@kM&X4)wSCyv{Gv;sGuxLEvQI89~DHar9PoxD;34=oe&;^cDvta;QMCox#ygF z&OPTp_nryVMXI@2h!6@qawGqrP+*{j%%Ds-W2HzLoI=Cg1u)laP^zQ!;29UMQ_eC( zfwu@{f+_|)Rg~QUn2lDOmC24JG61o)@Ha%L@_WpEj3g3*~6IfkNu)H7C7Bty#?rwb@EN7;>p zS!uI9KvcRIN>4F3L+J1nf}V)PQ6yl7gfI$I#xy7_42dPk6eblG0g4_0wlIVfM$|SS z5=odSGoesGG~*<2JCwCC*|3Wok9jx!Hs+;geUO|4sY zzceFc&Gh}2^5Y-l^In-799lQ-Ti$G9RZIHqgC?H}RrK%$ajkVB#f3@rt;MD1=qF3%gc9zTDhw6tO8#f`IxsPhXf6TF=p5rqnfBD}(eAiT3!Y%NNnPpTK zM+VL|))bk((*%cIy@ZYSf3|by!4;=Y?g(nv2k$;*LbiSV)$KE{Ee{bIKOJYCo@ZDQ z_u{(psFu55>VI9buUgpU_oh;LSH3K+^SNg`jxC-1Owxf@eQP+rHGWGLk6IGzTeI5B zymFdi<+SZymzQ~&-)UtT^_>#{@el%xhh4y4O!6aob{#9n) zvafKP?@#0AcFvu#YIWgzpF{`b%q&Vl?%z4D9hR~EEb@!1Ay@xp<@RA?Km13UVeyE} z@B=nVx$ac+yhBa#%Q)3jU#Rj~Gd3>!Limcd2_sVEBR8kMzW?m|U11l_Uig|DcIA#w zMSA5Ff;4-=2j19_vw>!x0=^jXkOckNkdW;l5H?p3En&9e5mXkOgo0@{?WE#l2+TPN zl7abbxX-1WbJ7KXDiN!td7%D*KSj|i-odk=;5=ew$>=Ot-@(&t3BNU3* zkdmM*SjiNij#QGgwQmfH62qnsYYW(+-Z;8_9qe38Ld5hLhG0Xs-ALK!hdKrfWVh>~ zo?)nHplgyo*F-2ATcJQYz)%?LDQ=IluNM&@huw?yImCvdZ!p-~NUOy{04#vyPXh52 zA&yil23{nH4YwQ3cyj*^up#nnK%Y1L5p)MM4!47i`;&-A04M?^axt+877O^WNQx!1 zWf~kv7#{=Cg%BleVxb9U0IG0D3}H=4^$gk7*f~}lLNovjR+ExUu@Nu|AqIx9&xSD! zh-s}f4MOWq6qqLjlWZaZR8K+T5gjrUNeyzF%br<$7;{6g-%t10<7m%GGVmY)3k^&V zRo}@Y5GoXXN6eRbnKUZFAR&_KxCoiey!c-8QV5+*bMfFII!(dPMg7iqOzcPBbcwNzWc9>W&6Kvy0N=_pZcAKg>_rI8;YB*-9LS1O;?5E?X^EP-YFgac85+eby|w9 z;Hf;nD*>+^Qt$K1G*mCheZ6t~kr7jt&eu(Pfmdx5?l8@bJ#u5`^Gxj1UllL~HeFRI zykS=Q($#5ei@)%*m}XC`%0o;1(2urU6^@EIwGEoKH6ZT7W{20BgPJ=rvc);u!hBmt zJUx5Fi_#jeI`7*2lKA6Dy>Qs1pShvOzRc$y?ac8Z8-pfoP=_S`@MdvWM)*yr(R-8^ zlz#2Z72Q3$vG&B-xVy8Ewszrpo-V*2t)ko8Ma!-jzfFst6I~gN)L+AnC!}RXw?I?aG-Edz4Mkq-FYudD~|CcfVW|+5Ys2Cw#{3MdhE8 zUSH+cONaTBTRW=jkILrH9pl5Vy?pWcvbN3@X*DhI3p=k~6?KmC`8Ce@+mkO(WImWy zwDrkl>tf7CqjQQ{W8A$rqN?kgf9&5{ylcMf=4E~NrKN;dZPn4=pS_TIxb4$~%IJh3 zukPVfVRbHBxE?$+-KFTSAox*btrxvK9uF2u!~nr)FM9id3y1@O>mK$27f3z=xDv^q z;bO8D2YD-iMqVAlk2WPHPHzC5y|R(te#&?(o^pDz4q>j&4nLrC@M%?yUtouSy@?BSY3Ih`b=GAZB5O|d$$%%XwKi{%73}#55K+K zNgEwEk4Wn^(X)4+t>3?Wk+pv3HyXlFx_#UIlKTyl+Y4@%H+9b6k^9cEjP)A=g7(4q z#m@>>pG+w`J)y1pSk27Cr(T*@@ZtGQ=aMC#FHCMa`9a`Ur=bI@rZb%xzh!92PzJ6>?>hE=R2Rkoy`{a+$yI#T1`B2)Kd>y;hd49s^NK>eD zT*G_&e&jWkYs+$XfpS)}NrOfudZPhoLtGf>H7yo?qe#xQs+Is1^*AAA++y)^Zg^F(cIZ!O380zqFw zgoS50a8o+Lz(yj)O3FEHx16YPES(e>l2X{9(%WflSS(cpR7LkE>mW*&^Fcu=>0F_R&MWzSW)*gjmBg zBE1Xs9{r%9gxQ60BjsY`9GA;#mWjz_mm7ru}S772xXA#8zrX;KmB zzStj3_R<96p)^#Ag3jSVK2$6mO5^cBgF2WGtT>2d_7wc>ws YB(q;pXa;vO?tlYrBJc$?Dqao!7daMkNB{r; literal 0 HcmV?d00001 diff --git a/fig2.pdf b/fig2.pdf new file mode 100644 index 0000000000000000000000000000000000000000..03d11b6d5ee0253cac1ff74c48f10a23cf46498d GIT binary patch literal 98309 zcmZ^~2|SeD`#(O6r3guOp;V(}C;On3Qno0Swd`akV;xJLDCmXWe06k+T`wrpdk z64Q)r?915382>Y#&-3|yzpvlx@8xyRx$m=G*ZX>(>shKC`1t= z1(AiwfvfVrg`GXzbsqVG+B&LoGIBC9N^&Z)mlYIb3JU*Eg=+|3T@zo>A&az$wy3YSA1$FT zDkCK$$0BX}NK{rz256!E&#(4P6H$jpPkau`Jo0fBl{WBjxBnk`6#vt_Iq2Tk*WTM* zR9eUWsk5#99o;LI{|bHi|B-c%# zZxhc)w)Ubl^BG+=@f4-$0-B{I|3`_#d&hqt=z2d2pfxG|Ul0$g>zX_gr2%~i{C|-E z4bzHV2YMZrx~?Ru^xs1F4IjJM+a6|sLjMEkp>i}S%|-u<9}SRy1;1_Y?&#|zs`8Jl z*PWjL&_$)MKLPgCvA6ZGv!}J?;|oaW#)8Nguh%bK>VBFzibS@$Gk0*l+-y2?oFST9 zO-q1|{^$*{8)Bc682++r@f1U926-hZW*JXe%qnkNOe*~?`r8jBzlF<4u6R)HBK8>vgXDC7YL@{(p{P2Iul!e-+^MdUzgDRK{sq_meF?6eMQJ(#I&CN}Pq+mQ#lyIwkX zcqv_hl&(f<`{*heD*`hqXjF~%o)l*0+);>Lgvmtbp>{N2`6au$5BD0rPWm?=v}GZG zYo{EXk8~==M&1Y()}b6_KXVW>>Eo^UAT*^u_#pMarA(H@S^MjUS?#lMjx$sS=A?s`6&00*Ez81? z-{uD)Fp9_G?&v8UWU_>DCpJ#;f6E73BDOya`n|ED_6M}txGYLeC+Ym!bQaLT4T_eb zFrsfm@k*H)xd%?J4+cHzTqtoG2Q`s%9%ZVxHUo}_W6g5N>NHLD2X*$pY}cpO?Y9-q zdMwbo`!)Ia_%58Yy5jqv&sY6lh5iyH3#Un*N|##NTYOUSy0H$qJBXnAnK|w^j*GhN zZHS)P*8YH8(>_bxW9ldG3pTti4w3c?{dM$U0~BNKM^;PpzeR1M;(m!XH;#$YAB04T zBTF{|ty3v9M1--(3u+?^R$0uP)f{K`UvivTn^cY_S2?@b@5Y(MX(r4zrkq-^KvD6V z*j~5ef2^Mz4?k$oHvPFJHf+C_x!;9q->~kC_q)08b)3IS`x)|!HUp5vct75GL2NA(UhG_Sog~`(+WZRxy=>v&AAgzDOdgI>kosxy^GMk`%Rpio= zu37`fS88s$E>ZwHY21nQOOXq)me{I$z$#eVS!xHKsEikAK&FSkpYssCPpWI-qB zMT)*v?8+WvpoHOee_&XRN}oNIh9Ho*$q6fnr`qxsRazQUqn{87sKNMQrWheux+>`M zw=s?z2B;=8X^1k55xbFN*quf?cHmIdqrON12p2x5=jF5jG5kU&{M2PQ42hd^m z#8#LKN)5XyZ0pBxmWuPLK8h0}yS)}gOeo9dSTn-<7ed4~ETjgbSTpeWLh@{wIo(C& zP29aOSDE@6)(uY361_zvG5wcUQ!!*Td)I0Qi(!ShGor6|DOJ&j5yF<(>{jX1E>y|# zGz73C*p%2RrF8RzaW@A zc+tn+@W4TLgiR*bjuA>TsxZn))+Pv^>(AraH`K<5lDvQEL z5bqmaZ0ox7*67cAZqp|YIwmJ@L&jM=Gs`yF>mZ}m4|zElgc#?TwPrlRQudV1iY$L- zn4nsGFoGfL+-a6SF{4~5bn3@fZ(FsxV$yb%%@*bv`P5O_Y`C;Y^O=2Bo2_nt&-o_K zj;^slw31e5ZEuXCN1I^4iIs%9ESrqj&J_RQiJHp4FJkJq?QYd4YdFRE= z<+U(#_4h*{55S$<>yc=M;%gGjo^u{j-Z(QWbRl8kUa!Twnf{N{8SY`#dB%ipS4&|H z%iix$x+<*_0ml-fbDaSd-Ek(!VYyOatklFDG@kyGhiD|PP3rc3>H80h>nD3(j)p0% z8~=F1Sx^S3Rp!2qZH4lmlN7$e-zJ%Ch4mY)fhMgTF9(bQK2LK)JkIRqUQ&sfj_bxUwGl2! z_dYV_G@B+kp3M6l#yrp>nP|-EPUj$*^TJYlkM4N*$aVW#y6f53jK8?%@o2Wz+$i>Y zVJTmccOjiwPHsgB7p~yNYL}7|BwDLXA0Q}7FR9TP3XRorjeneKAgS@fQnjK*e8kY@ zg7H(151k_KO>g~GgHXruVDBcm@ZOhXyI6IiOVBsWIIL5msARI*`<>+w&TRR#K?eL< zMtE%V(A1jgp(^(^vW&xq2|t>qli9Ue>I4`XP+XtjZa4=G2r{`(nf9T3L z(YynfsyNpOm!u>c$rd~{K@TEjc zq3{%?#PGKRowwYKO&(S!j((uv)9Rv?$fln2MHI=bq&ddxf7LRVy%b9fgk(r;#=rVq zv|+ISDxeng%G-AP;gH&8tDu*G?Ue{oKfZFE`BT5};l8`R!rPs7KhB{~LBb1opFCZ< zD?PUOv&ic6muks*#&pY|;SI&wt&g*=PuyPyHdF?xTf$$#YnAao*u(S+t6X;(#IDxz z6$Q91K7ILC2@aDiEDl(6u%PnXfa7&nV6!Y<6R-3Ee-L66K=Jcnc=HsaI;@hIY%t~{H&S#cBcS{Em-|yukJOLPn@7CG%RB zt$q!EYd|PU9bBJi;mvXE!kDx!S?j~1LBdve0)>gW2VTaV0;Jv-@^`=H<=uE>GXQ$hu$ zFR$N~-*W(2i~DluD@!k#(Jv)1Pa{s7?j7LSnmN5@eftz5{hzKroaBLINg?U&;9{Ct+8HVoK!qSapP$(hkg)a=0+rsD06;8wBvkH z@WlWwk10zy&n$Zh!KGnolR?Z=itX)rTrcsxGyQS~{O$Rly~L}muOm$+f#KpY3E(;7 zq)Ld>=?#2ZFaxbv&eHg_^IIb|$E3V~3ZGQng^WlOP5M~}!@;_rD&G})iZYX9M>Myr z!g;0fadn;f{f#mH>zoI&qruX_4Vto#dr+VNj}CL{`@H_%bAOFolb$Fj1?E%ooY%Fl z)}JH`-Y9rb?!!KCh1e<)kmMy+iX(yiDYKi~T{a?NoGZ*)c5fB^(<2qsmM#8P=!viG z(SFFW@}H#s+q3@ZFBEPrQ&{KohL?LU=q(gQnxLi!3#4&FX1l)gx)d<$;4%sucrBO{ z-_NvUnvCYDc+H-ya5fmP{`-f5;&MjQW&JocUWN2Lj0VakB}IO0`^qAOe6uG+EniBq z#Y%|NHK4}Fx!jp$T*c4rExo^i<-_dh{hiRR*VZV$kLT+%lV9C)0Q#M;%1nN9&p~=3 zwktnbLhkyf?{1pszZBQOzr21s{Iu69-Jv40d^GQY+w5P5*Tvl2hn40(zM2BpG1=9! zUGJ8jCtq9THt1XFqxwt1?T3a(`p9#&{Br%lAHBreu3w*QH@m&9T(m;3RG3dxZNsj$ z#pu7ReB$5g_Hbp%;+)&tf<v{X7_l&Gz17ki&;03s%5;Fa!hruMWUXisU%?M;mGpxc{7+6AdR=<$8SDpqeyp_WW;^`UbH0|8Mm`lhw$kLxR4 z;!W$5{EnEFT+1)!^g(y*D9yRApSHSlFWTZ<WJP%im=2)M$#1HNkHWEJM77UY0yo z;&O2Yq0vi3ZbN6jacKwN+Q;Oc!u6mivx%xRTc0|Vb!nHo|8sG1Hd`FICb`5eh4O}c zFya%vBhK8^@#jXhWj*VsPFE% zGUu-WDHgU}ipQ4zuD>AHTllb;yEd*&NvADZ5V}Ybp{wMI?LgmI{3SRCQ@YkxeQ)9O zoTo*@X0 zPhR=(d11ay<^DqPdc$zJ?RA{cxU85b`^c4rR{71^;FRxFCcQ`Pv6!0y;8Go6&1aJ@ zO!suYeGY%J-u5x>WIZfX&!AszhS@!O`3DUe(%<}-n&6$7%lv2{>WdS zIi_el2Uu`rF9kq5r(nCl>z^G-LUrY#)Wy5^0uq2>EgsM zOS6L(9<7&IkGH2>`Zw^e$Wj_7+(2a521$B%p04(}|LEDo^`qBg#sr+ad^`nWPpok~ zn{-Ws-ar@XEOLK+9k6!Ys)-m0QQNd68|)uhb>e?&xxr=j`(e(>+BC$dqUM)ZeJ66q z^xGaL)d5-n@h7j?q2;yg3S1k0$JMQ{e9p2-iI)89wGeZDlQ`$Y*iHIW&*8T@tMHk- zt}l2!rnocUM?RZO=7Jdo@{7**i=mUoJq?`?tD$qeBIjJXjd+75VD&EctN(~@IjO}Y z0oD0yem~&qW(_WDkiSFxOERt4WLki&Xk0LJDz*)gs3NtdCU2(MXw|F=;kqk^HUDK+(ZU71 z;!I6tDJ#gq2I0*7z8-p=e&%rBd4Er$-v#TzqXWmTyPdf8E;MU`Eee*~l*H8r z8zpT*#G_Qb;pcI!oHz0tBESUj9Yoy$krZ-^z?^pv^RZ9!#a8Mj)=i+f!XXLzyifPW z10Jb^O6!W{Cu?qO&@oBF@fSdy3zx!5dV3cHJB(B9;YS|bW%dCv9U@maTwz^G0}%w} zUzqQPUHIgpz;3WFZmiK%pl{z1amXP~AmW{)cKKxRqTu`+ zj<*jshh7Qliqhg^m)HCWYYb8c-3Zmq=yn7&9~(V8*RyHI8O});^i<=_2fMAsz&jmr zO1pLdn0ucmAW-d_5>V$UMck1O{f$?pC)mRtUS^7TIG>}o!Fi(^V~=syB)GT+2SHWv zKh@SQL5C1m87>jqv|_~)^-z{V=aypImGRm)y_PbS@kElgdM0^4iSX0hFwsHg4DRAp z4t??9vHrJS{h2#bB6lsY>gY8FcrfW$WV0izm~nI>Ml9=kJOVQ)I2RTZ?t6qk+KXWT z7&tb;^f*B*af!k9jc5G_!5<##SEX$&f?$U>(#jAM&WLG4#KOEX-nAhFRPeGB z7-9kfqtq73RVFJ8>^7Ls6x=dsC_^6*NhP`Es-$AO)k)8qSob&6Xl?#a^X} zy944Uj<`NN(;F?ncr8!%;U(2c z_ORqC+2e#GEz0ywH-Cd%&WXfq-~)uy5HW>8%HZE?K1HJd*+;emp?kLOkYSM4n(F3+K9%fcVU1;_r;I3 zAdo|&P#aIvEIDAr$ZWae1WBA&qOEkUV#vTdX_@%%Zaf-^Gi~PCE4UZA&Tsf8j~~%( zZ9`;q0n&XZdGsZ3>rU`_Q}Lor?TUZ@7X-)V?oyR6wJV42F3%OkuYQo_l?H_q{8J0e?Go$N@* zLi@hdYc))U8??0ewCKj}+@0iQeO9MS>4!Z~Z==C^4zfMkQ^!cAc8mgYv>D4-N_dL= z?&d!dfE=0?w8@PXJ8pS(t6PsBr~|6y1IP!I)`dJ<07TomRvw~9K9t$iT>rCNohE;L zom$|bgwf+ivVZt;r{>L4+$*_>f8?GGr>bO0MmyRz`MQlOz0|Fn3}d_!AHiQr2g;f^ zwOBV(sUJ>QicVstL&7gh{QNY>5b%H^&hOvMmi>ONEW?X+yc_T^f$R&7#sq|3K0>ek zOKnZ`_|WyypgP^XHlG$Kif7I|8yOFq4C9Was9>Hjis_0zi28ns4^a(}M8xk9=vXDD z52cPoTx4~isAP3U@H?tsw5y1SsS^AQ`R;>p#h7CvwLaP4&5vUp7{#!3$ao)CmP$si zEdb74$DDE-Flfy&K`kN)T`x&3&=B?=pNzUNHDs*(?~@LtLl?Il@!PKi*ap&Tc~D1W$IEk;F*cSopjkW_G3ub!&sbjNT%wu==^bS2ukBcd3^ zf*}EhH<1z3qP;>VyB+8@@4AM%caKoU-xN1}qDW($s+5I&M=BWENn! z9|@7b7Hnu8jJcG~KGe$3(~#C1*{VZ443X5Lq7e^&ZCrGB!@kV-3$ju#Hb0(U%GiM- zjA!xsEbz4a@W2IF#@Et$STHM4(NflX8#Tcd&UBRX(G?!!=t!y+SMKQA) z{&izi;3SbhU$*Zp!6w`>RBV^ z#M0sFF%kQ?ZQX{6u!CuJUaI02gCElplS9HXPEiw=uva~{Z52x}=kA4uu#zPEHq4J* zyjmUJUTR`DK8qlPQ;DuT)nmMNzOOEAG4ihlLbvKquII1hPTk?pE(y$)(*mQL_)@<5 zt001czb&=Zn&v{Y5~(LER`87DgvjQv)DYeasd03wsbALeSG1ErzhBG$=eIrefL(@thu*Ir z*5V!UMAKqFC`qjL#f8T+w3%Q7+)!Q&Hv?z)==KEt{-DwgdG|hFpm#{?uyU*)<`tbl zyV5N%!%jjqE*4p@zU9x~&_V%6IJDx^bLsP1eyNOzMLQ!1Xds&}Em)zx!ZfEsK7lBS z3XcS1*FH`&L~tLcEE>QDfcmvP#Ub;I8FN-1^9kd|wxiMVRc?)ERd{V7EUNM8FSao^ zhN)wHXwpu_FTU&ztx1v1{Vw;uj=D(S8JY9Q8ricOy#vjYMZK7~?dUEM3ft;lxx5yk z?%Q8=OhR<{ZAV97G}wd#bRyzKP_&$fc4l>K+^ALpgezgsx7>85xUfx@`xWc4KS$R; z70B~lYpJ(NtI9Q?5YLy5j|k$Ho{}}3?|<$VYMZ)CzK)!5*PVcdAcq{=3awcuhE&&t z<~Wh(bi*^5$!+Hw?or`tsb9PpqWi7x+kjX=1HazxR0^Podtb&k^D0~r%mIvn zwa%9l(j7K}2cNjVNxc{Ea@&xrl+FY9g0bjb44~L=$cjvx=LBE255J%;E{^(7GoP1~ zFr5OrxZd1$&*fu%0ol5^G^(M0^=xH120WYl&EL2t;MszjOV|M_H}o(GI+SkG6JUmk zjxvv1%8y#rF@jl3uC4sUx>r>xv{iZSQi*ZRkE`6Vqgs>M`I_|p)B@|qMdHbTZ)aTa zN5-p30&%}OF|<$$Bfrn7jBp6gyj9xYQ5)kZJSQu6VxO~djaq>0aUt#><+QfBpJ~&A z_8fXdpN@FBV_)aF`slH|m50_lDIveDuu}w;KC#v72RvRo$zq66-;oUYYDkFf& z>=GjOXN!5u^x)j2$nGf_&tpZI;eyx!vvq-09QNe%PVR*5_AZ8l+gWf+93pQWp zIr6(Lxu)X2g{j|4(2iu51n*!j!WJ`Xu4?m2IyUnM2B?f+@1H4t}%|YHgXGWi-*g}QHJqdjaHxORU|26TJ zed%8vh;>X{SN7jxr)Pn_e}+7cjU%3C!EYP^xZF}&pGZse;m3_XJ#SaXWZ{UX2lfzS zbgy1{gNj3GJ;yxL0J$F?Kh8qDN?2nuIA9X-K5H6xOB10IMfkv@e)rJ~zYx>y4&@Vi zBxlal1h2Cxa|zEd)(cac`$y98`qzZ*^J zWD}2yf^l$BTI7{OE(Ui1KAU_`O7tOmxi7HEoX!~Pd=~3SVG?ia;SMM_;3htX@)RaSNC{0Sp_f}0;CEYEUl%ET7d9!PTBqq|d z8Jw~eiaAdfHC~Oc1MR^0(tNY@+z{BhD_F6L=&t0TyTNgpS&iMNh8|TML;FbUw@Ig3NYW{hAM+Gt52)_I+H{ z%@FllEovn7JI}4u_sj!Wh0l8vz$mx_0kaML2g1T4I+OJ*d1?)qe+BtLwo~$$Tqb(R^iq| z*{(3l`OFdD>lhwdw%aEn%XS7(*FSqm4tq2(KE`p5nEVb_3LFz~OLoG%F{;l)cRqE2 z?lZa6PdUuoQ0Qll5irYvIVpMWpS}DSOO}Wajr&B1DhyH54d~vaRoG)e*0ai0nrM+v zcu5oaw#ZqQ24TpooqbNPeV9*~Tz{Cg4IkiWcgQLA9}d&ho);7LcuVqi6!2WsD-Srb zQyirt^|EsPr8IhPI^?t^IpuMjP>k2%3P}`4>^7q{AY6X%V+JyEQ(LTA#qt|08gYld zbuiZ-fHlu1e$~DFlat%|cpWd7i+~O~@*#;%Il?glx;Zg$JrYD8b=@VV7KW(7cFU7g zpS?{>1B48DNORAFxutEn%##f=<)^}8R8#3(8fA@V2Vj7PgnxDCML!#-AT+GSz5JnA z-Wo6d?4fx-ccDPC>D+f3@Tl7dQuW|{-G}-1$=r#w0K5xk5%i2GYgzs zzUFkW+M8?Ey5SA&Obp^gKnA+YNi?TZL;{S-Neevda>)Ju#4gxx2#aUyaQ$(YE=&x6}N|slc^OR*4#@c*Pq+M!_|`alJl~iDe72j>t)jGE3rIBt4%Q z9c>`)xJRAXKE9e*&>{&{KYP*3RST)r0N7yi4(5R$MhDC7*OLCj)fWQ1=ie6S=JYkT#WmB9f{n+}n0AvwcZ)#NIwmd6!U=INVu|z~SQQgiyyv0dW0yr7JR^nD5wXiE z(4C>v)Eb24cIT=Irzmb5t7i|-ijmli&2q-GWx(0sbfsWd4#WT##u`4C2SRu-)yItX z;ut-_;5P~jL)n-k< z6cyfH#Q7i0S#y!i$^gKN-A8}2S!jNb^3gWm3g5V43BS0xaKVww=mJ`+Z4m2_MWaLl zV#e!mJs7`!6LgnxRT^cs@ibH*P1=Hk@Wn+Q+xN}OPI=>E$a_O*e#7_tzx|0|y6l{h zrYMpjw-xZ>mJW2W?;$bpi~mKCGi(0>{3`(JF~+Psq)N?e!lw~O^SnqU^IlQ zv~Ak9?iPwTnQN9iZ(eX>I4b7+ku)caVRpB)=LPgck;^zw$@Z#U690vlhe3`xp_7Z)W5b zMj-%}%yoa8?eclx)2RCJG%|Y#hAX5d7}Le;&Z`Vz@d;mmRmV4P2zUqu3PdV_J?2;? zPqM-eL*3v0V6oJn*tEdwzI$A&r%$BUB3|6tK9@c;yozr5VZeqwO3(wR;J)rWAHP_!| zHky$wZL@y@!ZQ~4k4I8<{Y+*ogoy!Z==Q~_^`kpSO_NwHSW5mrFV!9jy~voV8#v%F zQGM*AWiC$ch9RzQ?P5rtrpo4FU`ZT+&EE}!?!ZcX=1T-TGv7ZCm_WR4=zl&=T9>8G zK1;Le8)$v%e*tLubwP+JFnUg@5llQ5zvl{b(6U=o0{{V19C^#fDh|B9rouPSQ+cTz z2vvw)iYaHwW75#TXmOWE_KWVYqV(6MXTXz`K|v3E92^vDT$2g`v8rBG=pU&yKrD;j zU+?DPE?(!*-^WPH0QaufG~Cld96J-r@53$cHD4Q?aVue9+A-{m!*Fa7pGLtZmm%r< z`pqH7ITRRCM08DmvegPj4mW?@FtF8@*8sF_*08Gn99nS3qIK@*XXbqPRgqO4)FqUT z2FvkPuGklzO9L-1u7?kzyYh`EbD2xQzTHPXul^j}etf@(M^-L&2*fj+cgZ^Szz27Q zxyJw1wgt`iO68FApzAZt8JxcO{@q}lirzE!H!tJz=zVp zJ|z4;e>T&6lv*ll;HMRs)%SAXo6AY3^9^4a<}~cn=VQe!HN*9kpi*Ka$e#oY-9TuH zY4IkE9r`Mz>xf494i}b$UH?90HV4g`13vjcM`Vqwpp51Mp4|>scJGomlH4@1Qx9I6 z@{UEE14T*Lf#(^tC3JoDLk4gF0>)L7`^3^Uz4`N{)g4wK9~YU~p?6<*v1^cb58FSg z67n(l#O4yOa>#Z2g`=iF_&qCfp5GR+h+p6pu0}v!#e=fJ0$A_hjY}ad^7}vu8j>N< z{Lqw~3>&-N*1xju#I@T@fJw?ji~EmDbAX~0A=;7bxJ&cOxES(95dxSk@VtC&^0=HU z97Anl#)e=x8vj_2a{Mf;b#D6c)dWIYxy7XVEmR>E5W0vILpFZl$7FGKZtb_y#JhZ1 zgnWJQ6)998#8l_y3l5=I1CQVP7{6yp(9M!^Au0j*se9m(GyWKX5E&2Jn;6d@hN!Rg5646SA` z(G@yCErb`>VC)jZax-y2lS6Kt^U$0HP791#DGgT#Fp!V#u#!d%%lt;USHwi(YrOGn zT!7TzOeHv&8iU9~z5@H<`x!8@!Y&LPMp-!ag4ZF11WolKH{a4>&+ZJx$r!cT{!70f zBs?$)3LFPDt+TM8zo83;zn28?qC{pav@HO6;G8)!*5xXAMA|tB5Rwg%XZYH%gJ?3H zTmZ)3OAwG^v&}DD=QQ7b;Y%z&Gc{1<81~aidvq+t)rU)(TDu1H0p;ag#tmg&blnTa zD>!M+#=+kYrQg-bov4$+-FqX2O6+n}gIvtVd!*%+uw>*5Wgkz;ame&Dw8NE56GM9}Kfa`d_Rd)=O~p2r)YGY0!R&~K$jh$pGSjWvIREj93G&7c zHmAx*P@q^kIMDmOG0ke6WpJ9yOX7ZN0yh)bmYrn!7)^b9(Bp`qKCouq>2S^f$A+*2 zN`o)W)-R&KKmJYsbVo$onA?bnGDS#B2!UZpKpChKofT=Rck%5;oWc zGuFJA<$#g8Z@snK;jC`ZgpVQ)m_Nc4LoDIJ^)iHvV4PeD*j4{3RJs_GX-1Iioq-wv z`?qj~lm-7yWCjO~6&JlWB`Fp!?G4-+`SSA+bcc#^RyF(X-lqfSthbm>xX4;}9UCj3 z?~S(jQLgGSVHfB*-(Hxe{e|alE!!b&Afqaj7T`}zpzcw`b?rwJ5N;K^lZ2s^ z5T=n7m-Kh{JzVOHQ8e0x0auZLxQwKoJQeHJBbq*yezU_^0(H0g>{W7Enh(vL2PA;P zS&R%LlXf9SbEa+TO&BzhTN~pT##1pru>vgo0Nz8Q9Y-4qE5D=qyoVL^>a1DgjTbXmdj7RT})D$&!9v^G;D(|AKs;jlNY~~^wa|GC1AGu*|)9G z1Fwj?_c@R75O52a`5t5qoL*0hDe=fg1w%jNWxn}4?d_V2g>L6RHp|HxCw+5dD$33C zDSPz+s0TjyJIT+itO6iX7% zWc7$-kS`rX!^ny0R#d>X$GL;JJ&f})_oyehp-A(lxA-p+Nwo8WQKyymSJSj6(wqXI z`%{9Vic9*_nT7noe}8xa#{uh~o&T<{8SKy}S>aXUU3KicWIvOE2KwAjUyEBQ8DuGu zTPcT<3UU?yN0pwBb+U?tVRU>-4>J$~ge|+}RQhiSrV&Zie1dYJY}E^(KxnTkrudG& z*I(noEyWSQA>=231;e#tAI)A{#4ka_e;>{66EC|wr0kNtU-@Jy@6eOwSOLPoU7%9t zdrlerE}kb1PmZ5(!@m})?J zdYaX9gFF$?2Eyg=0c||>knr6`?;mgKe23Pb`xzQOs7s+2!|0QK_S4$Y|<;H?VQtLUS#*$6rT zWj2g?b`^bWKlg?JAX#c7H$l-nHC zPoR50B}}CUc+>jkR*@3N9N1OYIVlI@AiWNt@fA>>o5 z+t5>yPj}6bPhu|*snbWmldgbreL#7g#ieNP$Eam{z;I@g~ z8)k)M(bt32iXYQ=RxPTNjiq@A^yp-EuP9bC!j9ya%^`GC6*+>kp}_B>W-^#?l6^{p z*4MsNb|Wr}c+Quj&*@P=p)h4ZIioAzz!l^J*kdYDJ`xZBD)=<$^L`3|>=A(MqTP0L zA~Tt0>;qt^)su6BamMTzCW1I7>1Ttc+~C&M*e`Ecx<>KJL$C(Q$BD9-?`2q6>enB?OdZ`h~Sf~EU1Wt7FQL0qUM6*X9Qh` z^XL5xU{cU0f<5>Xgushjqvf6Yw6_ebuYCZ%t)T(2DzuXg*@tvc7?K|cm$9!f^jT@{ z#E7YZV;fBB%<%5zn$N-b+Y>wGp-dwyVC)xB(s{u^91WYrw4dk!9%vWj^w4B+%5th_ z*#NlI(JHF+-*Sn;jvM$UsyLqT{9cnrEdgq?mgtsxxGeOgI(jdm~qWOeUJ zdy>Ndm$EkHJFFny{0xv>spX~OI5AXsg;%c#+ow1mgZ+zE(B7fg ze;N!0n4B^j<`|D{HO~ecVKg!Bg|dn<^hg4>0eYt32nf&rtuxHk zLwidBTYiUjoT3pYa0DyL7?UU3KNS;yR5>6A_pd^LtEeV?%-@;I!}EgZmmY3FqkbgjrMjQj=ZWvAW%Ychqd^h zXy=$G&hebBVT^It&iR!x7vP4wWqtn8U0}2IV|SYCEoP3Vs#`j8r0lWX@NfvnEy8;;xHpU|NkfBpBipi|EtEx?DCgl>~C9hRBd}Dyz0AK zhMv#@+YYf$quXEGP7}b1vLaIONdDvu_95g*J@jaH9?S_J zgiyTv^R{Mxr~O3vJGvkzG7I+u{?`Lk)4x3Whf*3I%At4?9dbPS)31%|cE25D+#C`C zpJwAaRNWlHG~IPZc+UQe(D}77vJL~2F8`3o-uusR50z1JT@^L34hc>iomg$J-vl&U zNQeqP70lGi#fYT1-po9>bk<5tcqb)HlUnhi7s+!1Y2ZsOH|jCT!6zUOK5n3*T3r9B z-H06C#a_X8$!lt}G|i_)YEom$4J3IHaR{6VJ#`1xt6E`!>!vWa<=YjW%En?2$|!+$ zRJLx-I}VACt2^6k9KiuI^7j)HHrgB0H%{+zm2z?fv+SOy?&wFWhUhYZTrTli7z^!Q zE`*n9{6zi3SF}VOgF2XDs9o|XZPW>FBxL~0Ox>~19u-^qNAOVog4>e;V1^w+t@A&z z;#(D7-Eyi603KgN0;6dAnpi7?OPCDWo1wL^q0=PNxb494shxKNzAE5Fj!X@k4a>(M z1He-A#Zc(51pNebIql_WG&B5lzzJ}wisJ2Bm{q-V>;+A<>3c!!`wI&dhS?L4PwYlO zvqT4usrM@jB)|^eKvP+o)o;O$zBeHVy3k;)Xowa+6M;fw2=uz6z@t;-1 zhZ`P6uVNGhiKTU11pC<% zgpo9+S;chi#Hx_zt&>1n7+>ROr9^cfVAlNSZ|V?nz)|1ANN%j9r-FymX|xHgNwS{= zW-Ul5=X2~~mv0v6677ape#wOo3rG+&aE1w&$biU14l8>1giAG{9Gad%DV#g&-G-j~ z^d-Mpp%YM0y^7ux%(?&$jUHFae8{82RmVC&iqFrF=F%K{; zp509Wj-edFeE<9^`9r&ZY0U484PnpRy!MboCEtQ!>Vx3|3u0o2#BcfvmQsQM%C?Iw zKndytDvc+!#-gRM6#D$NvYLLc2$4~UAy3>4wdyK;_GvFaFPJnxB?&<73@Fls{-^h+X$S}D`_5Y8l z_YS1`i~q-MAr#pwAxTCg`y$Cmt`f?YnZ4&VZ<4Zwloc`(S7fhi?|X%;?2BvfbB+Mq+M}=i*&* zu+7fxf0`ZH^!pzhul&ol2Wj!K^Ue{WhFR(mkj!sGNK>LYbClC4CcHMD>WvxU&egYq z>+jsm(BJ|*bTB%<3iUDq{XGAOzV2_{&YUv!iP|ZwYp3999yA_GQN?h*tK**skSu{Z zVD_d${&2Q!uvfDp+7TAxY-LeQpsIJ+V1CZZE`|SOLDGd^PGNyxU(^$I5qy!<$Z#3F zP0S|%J5<|wnT}(F1(j}WQpbKi?|nsU@JLJ1U)3tP<8A?jT&b1$n(7Y5%{Zl=D zZI(CV-Tzc`rYA1||!?ogt5{xG3v-#cQ2-(&$X=9_k1Xj(X z?ZdhTfN7Y=!+zkMl<&ah3%-GwJuUSMib?iO;YN25m)nql7d=`O4T92{+a*-%SuH=Aq+{yD>?0TY>hNe zU)%r|)|EGqr*2A=-c{(66_j)VJfRDVCDda`3cNNa*~Zh-2M{x}6BIBc-V>L=s1@bs zq~OodW>$78l|33R*XmQ^iCroe8Y+-%CO>#)~Cel zv;i~!TEmMO`hQ|h_nZg@%fv70%rKUBOn?Pw#?w;ml8H(YJR9E7tW46a6`eHD!%+Dq zK#?LYtPaD~tef-u*t+dt{V2&*{qo19FB|Rf24fMXj$f^Q8O$UFdA(~<_V)P~_Ae4S z?Qvi4i{H0bD}Vj-rSkW_h(1;*CPl7C;kUx$M*KKE;F!L2>-S6GMb%f5;T-0^`yY4Z$ zN!J(H0~F>yaJ4o-e(63FKQNDAcOft>QLRS@4?+#OJGwgcN{BOW7{GTYFR4(iF@ZVG z{1eDksGL&*U^P@*8Hv($PC99=b;UibyIabe@y~TXOih$0(8wznpKjLvWf`;`2K(rj zMtEXa|66$b$z>J6<5ngy&ANCSTzKj^3xrT)9(cjZP(ECV8F2S5uz`3_3)RAJL33ZO zzfU@f5P5UU4?Be)8}rUnnRY-Q<=B#;ujFQ^eB*(1L@(d$_bQwr)Rqx?iOYJQRTRnu zbP$>VZT$c`C6_RgvD+$D7G2BCb<9%ob$F_k-ziLU$DErsn9{DkyV6b*D26Lg4DbLw z=Jiel@Bz0$`maWxma!pE1Y-BrIIB9Wu%ore4DdeuLR@kA7g!!CbLF7m6KjZ-#F3m5 zB@?oFIih$QtaFz7-bXmI35>t--Tsdgcpag4Wn9BY>pJe+_m-PDu9IY24}Y@%I{xFX z>Sbm%(W}xhU_&{_(^P_2`d!Kj{&vaFxayQ4}8TtUG``S*iU?n`|{Ho;yTE3KzzYZxHr&)=)3wVzd&FUtiw}C#~8Lz`6dP37Y+br z1Hq8bOeX%KC=^W~F<81otlIbMhl9g@IU@1MEjdfyw84iPF$a{DAFr&jG2=e~a%pWX z?LSSWIYkz_atO@JWDbXL`P>DlC6hM#m@mKyw^hN20z;qr7hl=Ck7^z1f$b85_)@!s zF!fDgVxdjQ)`O{ejL|W^JciJ;v+R^OFvMImr2AXQK_6t>$p~`Qrj0iEnSJ76M9+)D zM%p^S+lT*2S{KHRZov{I$37-=9O7%B$Xim$j&N?qE#gBSY8XM6_g>2;p!{+ffT5ff z0tt$(iJKK7)+?ZO3H&&zWKISFe9;$kwQe_oQAmM(F}eU8C}R8hc9;?A7ed)UM|8kk znjHRSz8OTZp71P(Q323b%){b5V1pUR2C>%X#Uzm#ir>BvCXsAeBZ3xc3`HFC;{rXkPw zn*T!ovj(5i=vAZRe}s=P`W0;OM|U5v*a|f1b2B)O@ikm6u zG-n_GehT&FNG~TVsuBw0Ry7Cc6!R^z>5M8LDSQAINgnP38pC-yf{dnQlMdvz}|!mt$D3=WCcNSn6xd zz^l|dn;l+2Tz1a>zcHJXuLil$yQ`Ks+2n?J{$;^4dJqpT*cDfF{c}ED z#s2gpG0`BB_Kf0@K-Jvu51;6#gLQxDA4W)WisJR**77l#n4bGZ2os3s%|G#buJPo8 zydicZxX=t3hS7yB)doEf5orA5b9YHl8pPt@R;WjQTM^ z>J6336F+Bog6my)H)z|h0T3E2Qs0GBvL}x@83J6$KXdZ&CChiU{+3GHL2vYJMsR^9kWyhw3#=0_Of#z` z{S!z4-9Ld`9=f}0?5QGa6tVeTaJ(Ok5}2L4#Q?j&kV425Ki*vPpUr91Lx`^#*Y8I4 zTvgHEJOy$A(4#hS>D%$lw6_60?1^rWJ!i9wD(fFBNw{YXSw(qJb`^c501UvV9#f@u z^Z|JUa6|ESoxr82A|gvw#+SO)ITrtbOrU0ry(|*01tKShLIUXGI}=BD-hWIus$SKUH`L&s#6ERLDOZhTDzbV z`8cENDFJo>a}5M(quBlzg?ghr-|7r0>r3sE6l8JzjfCl z9yWE-2D>M|ECJ6-;Zy!UEQe}IE~bIC@t(%Jr2t|h=nICACD#d?;Dk3xK8yd;XuKXs zDZE`CoBa$=VfHm6IFknw1*!^tn{2-o`qq2y$*~>C;$#5q+`0?cshTwg7)7pFz;v^Q zIXT3Y&oAV0=QkfEHks`>?&Kv>N84UDV)Y$YVfovfA+9-*fu4(QcYrwoC+7D2Vczx4 zAMf0kncS=bL%3n^ZS^SrX+wiw!kRGW*T2Q^px5&9niJ61q(=88-wBd>*X5@tz}y2Y z2m!!aW%@;M8-T$kE>V1(ylU<1o8K#`Pvtz z>?yLzkOZFHcl_<}_pw|A`~{^zl1%_&Nh(arbY}2i;z-T6)?~oTZCP@Key<%sPV>dJ zxjIK$7(I_85aB{AOPfZ=iRbvJj5Yzl>65-#iXr;v<2>D9N{&ybVVEa=V2>r|45Kb(sI>tcoPHU98G>mS^P?h#?T*!aYB>Tp~8<{HIDEzp>o86YS31>0?gFxbl<$^A8~kx?d?Ugn`Ym z#eAGshJ@&z{WSnN#+?&B2@UuOcL9jq>0o`<72?VkE^6uIU=I6FjKQL;8mnf-rVc7vZyROcnv^YdHEc zD}9@Qa!(S#w~L4hKGK@VeepXDfVrc>wFaDr_?Mb!gA+o89Tvqzr!tf8`jviS+F^CM zX13r8I(-YN!F4V8lx4<&QLh!}A{V0J8uGX)DiRDtt*XW1$$=ImpdZI1@<)U&JN=P; z0PXK<+iuOd3??+xQj{fPoOrZy$pAU&^^B%&qM3=}JBt5n$f{NDh6S&`koz5Uk$9wknl69unN}E-tmmL@)G%{}1Zv&Tfek<60nM#^+ z1_aOvE{C=gNhK4ae1UG0`|;NFb48E=RL~3wkCl@gpQkx}Mo3LbE6446pVRjLB=Vnb z*enCkyapGqMj%1zTJys+pkZ~QPl0kq6z>5{_BT@iCr_h2V~d2c6!Zlf0^>ab&djy+ zk8U6E1})*|EPgw@0pFnU=`b=u`nEAsJO#>yhVD=L0DE)3ET9fOClk~D>wy8g+zwb=$@l*(O;@z11ezE5muSIE zBL2eb<_bGmi_cg}KCt-XpVS;Sv=s`2a;Q-!X{D$@1ey~F{ec|mVu zKptd3{2q!&Uh34Q#{>5z zX4E{sD$q%QpfFtke4FA4iXb^j_LGm*TGZp1-=Q1Ub1xyNuX#bCASiiGz^t<5WGzgL zQM;CHYlY7}03H!!q0@!w{A)-sP>Ikrlzo{VDV^jTDI0?cdJdN zlhv5?%1}`+2nUkMpMl^$EE0@&F%Q!eR;Lx)Q$Pv`p){Qka5d?G3(fgGnwmH}_mJD&8j1r2i% zw78gq-k7mdnBNoD>;yYn{#P}TJYAb@D&Z89ggbA}?{_1LcvU8csM*fC#>E1GQ?X~r zh*cwx;^-q&OMYHl3>+NP46k-q7zd`mf;5Gxy;pdjCS4Wms@xLC!BQ2%)7UV=uKdb5 z?NwhnJ;sdKyxWS9r9s{wECg$UNeCu46;wiJAxhGRN*BISH*+c&_g-JfD)r(jhGjS@O48 zM+Gy67MmjS zCwX3AN^Fv;oPI-?L++0iX$qlqG65&HH#n920PGHOYqzG#y?rn}2;6H8PJN8@p;}|; zr(U7F)=bR#p*xRNnHVdG>^EBgeTe?r)J=TCx<$lc!bHdd^*<-{>wS8RY2!1iBR#mJ z79*|U^@AP!a5Wi4(Gc1coa#(UMm&It@naPPlR?FlV}ZGTLQ^FhMZeRh3$GT0N}NOL z*@KYx>XZzG23~EWO^@nywO`D&Wdv1kF^=kwtsVw9(~5A8<{=?XoFm3cR*LP0t$dNw z>fq$+KNjhNj>I8m=F>!@gS9ZFk|x9^L>weK9*cu#} zSARy?4`({hq#?5egbJ0YDkhEy5r$(>hX+*1Y8*M_HDhW_;|?PLolH=iO}VU{b?caA zR2;z1Wy9-FS?ah^fNL&CL0dIj)A;t&2&FTKBy9!r1+_Q0N(?lew|`g6oO)xIKIwCg zN&>PYS|;@Iq=0br(ex~u@;0$s14I6zlr;OuCW~rrqmNm}!GWOEvQ!|NpF?tr=!9`J z8Oa%)G4C>guRkoBEgW+c!X)DaHmD$x0bCDJ5bUg@M^zGyvGaqUW^H6PbUJ7;F-L2Pf;1>1vY9QJxp>j4f-z$lDLhLsh#jj;* zi_CtvA=uO4ai%Gezm&-=mXeIayMckZj=tUAjJD%X>&)e)L_FVWq@7@b6U#@hjSku$ zo0$@wTpP_R3dp{&i+4HQ={7W*CIb3I=W>*>%zw&nswC0T^;fUzGiKh!NQe<7V+fBGMG3@90X)^q8{7c-oNN3)59XuX&6xYkeU61?7R z*y`OZ_10^FuV4K?JyrsZzjF>F=aBp-5UJwf(?gdRa2w6$I%x zwmF$uocV)wnY7_z-j5MI2_r3Y_TOgQF4|iAYXa7cNss+Ib}1L-Hrp1CDxK>%A?s<$ z$Oe*w*UpWh#3{x{dPkSo6=Gjv{X1fi+r%jn=Y-vvH03aN3T9`3YZ3_we1BwOL(1dg zB3cHtvK|gkk!bLSA-BIrt_F8W(ykc-Qvtn9MnbmSNrhbMNNMi#j^d`a8KtITXTml; z^|ax=dNR?n`pZ?@|B9`z>$`}Rp(>kgW@fs2S*K+K$!OAh@8BpkX3QwW42e1(|Ie|9 z(k*t3T`HZJ+#P>JMQ`3qKY)=^M9Vr~ou=fM+JJ3dm#qqOUXMXS58P3Q=UgEfYuMDy zGr1cIOtFZ1xXk=!AMr^MT&8!ktuh+R#C+1bcvT8QMBV`@^CgLnK{owFJ8%jvg}~)I z9cbA%|Mw_pt7L6B<3q*nfIpi3JKE3KP@n!#3dJAkjEFDn=dqM1%n2Zp( zP3(+2DW6LTNs3Si{Qr)f4lk8x*n5?%-Dk}DofqmhUEHT)F1eW~jeA!~ud9+suggxP z5|Qi7n2m$m22CJ8d`*=^#&^A~k|z{K63rC?2fbY$KS?RB&sZMkf6?)+{%9EOu}+^_ z+=;e&p4J}t3w%Du{`ZjqN6v&$oF1#-rJH5~pKQNzN-=99pNam~2vCSoUwtS@N?$Y- z+Dk@X)Ot=KJ~(tBKK9>7XY2VPg7VoUx>X1N0mPo9_I}Rj0v)Ih zqDVlUJ^tjeVcMMIPHJ_LbSv@9FzfjNqVgF9`<;K3|G4hs`YDttTmT$KX@%rUrW>R* zSE-sm{-cwUeZB+1_L1S}@5skR(~g}ACUX5J{d(H^+Iq2iPwQ$%zN})#8woqUJ`;2% zoNlO5C_XDZ0e7+#ijnB6N74rwYSBM0uDOkN@MHc$c{@gfD5cxabuOXY{kD;TCj!_@ z4HB~xxrxq=5*$sHx!DM0qePsuLj3HNn-@VddQX|B0|z@O z9tlGI8$vb7qnw#|*qiJkcwPvFB!uhXm9&eUcUpK?%0c(da149J?Pg!#fMKL-vZ?m+ zXlmFbs%atQKd*Fx2r>v;7t1njFnsh;WRHp+0ho zuxxNpNq1kmMVlBX18f!P{aAX}ifXIo^e3qgs`7yU%!6}ulM!dlsZ<66fs;3;PoFQe8#UieU~dH=mGMyuO(+Ei?2p2C~WE-ygxQN zsc$el4@1oou9I>;UEA-LgIg|y1dQ);{Vh-O&FyUVuAPIhAJ7$h^y9gH+EGDe&BYlNj zxjm^h-s>($nX!4sIAhxImnc=$gs|q#F4;P1`eQkQFOAz$+54tG__p7q`J&pJy}#6n zqX(sMqnsZ7KJSmT!aHotlKmIAS*WJhv(R!kiEFomd?}8cPJ%WhW8rWbnQ&2Q_92~a z3lT;Wq6S$P?x5&O(SCfx5m!<=Ufii8Iq>?dVnU1P%xDp82r}=F;Sm#q=DqNtUh}pI zu}pvP(``A}qTduZqJ(Q;T$$iXw zB}X1i_Hws87-Veg$PQe#vI5^%c_4j14H9#BmzmDWQ^}WuPV4-@Y}!IGohZEo{~ntizdE!{P0PYx^$u zWPkJ@ZD==79E8(ZbHYp9Mtitq;r&vlYoIEuK-`Ri7ptX{S1#bs-K3Nb);mV3E|8hK zFhv7o3`bdt?~SHr?~XU;JGcEYax@?D6Pc?Wv4LCMcrkCV5BI{!7^5-jf(R;wuAsGt zW4*+rj8iv?sphP0ZuVBMi2&wY;5b*l_=+Kky~rgPfVMN54*V2S37h9-|U zn=^%pAMQdpW(D^?>(4Bh@`6ui;Tm)$l~8r^vqqo;?yKhZ6A$^NJXfpt-bV@A?lnyw5t|n1#{vdYD-u>F%ds@L7>1AC{O_~3_HP}eBk+coL8v;L z-@VFdI2Cteg5E*EX1T=FgqS8sc;#a-GzcZ~Dd=Y`t2^4yuJr#uL&hzQ`O^v&$vumjeB@NEObbCT!r|z$_!qaW`7ict|O?GpNB+ zb+_x;x8A#gNiv91)$hMQi^!?->tL_^J7o{l+qunQd61dW{=p}Y=9Ks1r(I38k%zW` z5EH#MOlH(sRJ^T{lN7XPlY2CrJ)!v4c8F<;H}fOoy%m_?Q!POLm~geeStH_zD7pdb}mt!R%sullfG&%CGmzWBCRxK8MVo3&`_#o@&ocXad5%| zPNMjiJywuP^Or#vxfFJZ(>1M79_3BAwhunJF7Y&CpA-mEG`5hw35tB-Je`@W`-v@a zT4}^a2X;gHGzxJ_5-Za3Gk53Z%JL;laUQXQeLE7*S4s_b6D28b!EvEBoo0I3>C#Vz zqUV_uat6Qdl8QYC*V8?48rF=IZplySt_^uY9U_761U+9y`7Q2g$;JfX)K6}6U`rA6 zq^ZYS^eHQXwjNhlp3`SY3-vsqNn-2OlsBtsUDn4}+}FKMWtl$=TX;g$8D+5ByZpR( zE3*2jeS@3SyK4oG=+yfQ34tcRG6ZC}wwNJix$nw&*D22z)PCbg#qPMho_#{3_B72#TchM8?m)xYalH_7UekPHnKfhe<73vN@V9g=$>2;bSQ} zZ7%fq?mDKw!nCdBe-sg!!@%xwMPPt`aKcGc9+UGRZAb)t?G~Xz4qP*UI+4%` z+c%48WyDc}ChRm`E4<4d<1I>74=sv6o>1q)+4t z1^+~3qh<d~1E-s_KT z_{gS-mO5|E{hP?N{V+WT-q}HZ-O`gN3 zcxCGslubO;0Q2(_y9dg)s1}|1ZQn5IT43RE<_eA#P+y1J zQD%=Dp;ffCI|K8TQB^}Jc5t6#_z+JpIVe5YgkFO7)alGF(Vp1uv#$5#!K1Eb(a)Gc zsm^p}%{CSJOSXuOM*8;3;##1*1ObHP{-auk$9PE(aeTh=5k->ezN8s^lOzl-c{|ZyS?-;HHzVVcna(rgGMl)8@ zQQ`8&E3*)JNywP#EH&ijTB8v2*KbtDZ9JkJt(Vr50hb^gt!j5BXRax7g>b)UZvPv~ zJt~rH6Y0|5q{JcR10Y7O3_n9;nfX55y-919^#jUX(blBOsA?h;F=dmIsb`>ZZf z;Q&X5CAs`uzAJ?#-!EQVT>A+>s#o#P*VL2{I~iANGQVaeGQR<|633Bm98{r)JYQLw ziuhT$8zWf`&BC3_Y8>4!nPX)>mWf{Xyk%pqNJ=A?xkLEwp?sI`l$oE`P7-Jl`Y1nr zUXJB!E(^Bn(w~G>ex6>_>$4-WRebrfEmDc>9)ut+6Urj_ zNzz}G&f_LL&h~~I&7)8)>;JCyL}Eq#n4)i$SRKxHnr9~XUe2>&mp1F!=!$Dm+-#0( z*-V;D;0X_L3gaPnSdf=~z`Ud<^*>)L#w)~}=-cF!%99!kinHFd&I^|a{e1*giAGxt zrRk8IQ{8b$saY1D-}Z%mqICQG5&=z{rXIQ+nhRIV@_7o_i7wY8zOxgVNjfiFA*l5g z3?~@vGaUOv>nXZF_=on1?(JH?wZ)>fmGI+~?qW_*p!WD{U~cK-4μ;Q&jiPQzFQ zg1G#Qf(5cc@q-0r4oxdl$PoJZ!NTaS;aD(%FJm))_FZr9kf7b)*H!kbr$U0+GQ?LRf$*P3W7d~i+Qt9n$& zXspF+|8e=Cg(uW=!pyHXo#{)Yp#SJ}{6V6zy?*2KvGP~mu2!@ko*gE-yRQisubImp ze^UhoIV-YFUqO#+C#y>N*FRr8Kgh|qGLPp_*qxN>%-{b#3T|>WZAnpou9#Cj(cVD5X z#&!ALO6nZ^KwX9j=OOncLt&{QFyv-PS^zy|w|7IHu4>rE-eM`&FNHRd9bf)#YHQv$ z4!Ks>)jtJZ2Tw-tPmM{MANtr#4-Ngh)}2DL^{Popw0P|V*1*H1#J&WIj-QaOSk-!J z;J+X%R{r%(9)D!GJQ0&JT>GW}zUuKsIZm58Pq{L$u54a8r^+o{^Srm$HZilPLy-R3 zHcGt%IQvKGInFdFHAwvjYr;=A18`Tu;VaIv{!D!k(U{t54_2H&;8Pp2dtioQZL?GR zyWSm@!@`{V%M{W&n<6uVno<#0=&g1>Uy^$(PI7SLmg2OC_tq7}jDk<$7S*2Lud)-Q z77k@XBD>H*yh5_4$+2rDlI1Hm)d=&AkCyg%M>~2kL)|Z;-61z>XeM%I<{Uz|53eUZ z7^OX69JNz8veR~czE!-|^PvhSP29bRB}XcGj)DRjBlrSe6{HfHXwI6Zh#1yXtD7R< zh_g71ajh|E&jelr-B8U%an;41k*;)-*AhD{aAB#fFRjs(R0fsEBB#!q1R_i{1EoCU z_sG|2U1g5om5Ysv@+W49U}XX_!=U;PZL23>t`sGsD+50e@ArbUn7fnP@D5YVbHQeh zNLN(cI-cn#c_r&N;QG#{@pP+{l^dL>*Ty+n-60 zFb|K~XUR(?BDYCeSpB1}kl9?5sxSjAZ0_NZZjCEIH@+d4>M$$$Iby>D+d192Is4%> zVr%ZQQbys{`hE{|J3&Te;#3pFtfgC*vFJ_g;&kJ^k5!}jn5s7)5M?e#&25j2Gbp!1 zP8SDQ{%IKtH=8Fr!o})v7jUJWHi_<(DC(Z@{gvb6;~tT5v|I>w z!YUrpR2lVSKzK8l5S56|GubY={@iZ&rD5vWkwVYD!W2Qb8tYh+Q7UOH=AJ?CeQso2 z)jI?)Il1HH$g5GB!uO}7uFTmQV%7ku7WC$B#jYxtpm$Gpfzud|Co8W<#tuUC-u{4k}B zO8O+1%keH@mA*PfiJp`IA-lRg?)tDykS~KmfzCIb*Fop?eDT~@4xV5IGJ6y5tg60j;7H%IdHdhuh(X z*ed9_))UAxYml<(3*AP`Rb2hhA5<;)m3j`WVJy_rNvCPzz)Nly#f66qI;ja2O;}1lvB?^ zXB*zg8;&-e&hg`}Acj&X%Gw8xR3P-y6B#Kbv*)QG8E6Rqi00@%(MBgX;D=l3{jys+ z!`7@^zkgLodX@%l-`$`@Ea1BYtJD{7QO8Y?J=!Rb7>WIDFWQ6*12*<@vq}Cyw)lr;H8KyS)3t zGvhH&>_T;K2Q_RKI>Ofb(EwiX)u)1N(C1qKgQV*By}r1vqNbrM^6oE|@j46W`w=DF zbV6;dQTwk|w+i!bZ;!{`CbpBW+q8Fp`;ISOpXMksAP6N@t-njZ9H}OhHvC}w?_pY> zZ)ckO3lu&+zR0^F(mM$s%HNr0#(H)>SQH?o0yMi_NPi^gkhA0zE7?i_f+3q=FW3$3 zl$SPq$amk}`J%Zmw_HpYkxFA4&-h=SgTJ`kwY~=>sN$Qsf*{1E9T7Y_PpCF$;~%wj z>NPo1pUMz)eZB(3mD8Qpkqxm+H>!_` zfnd8(eHqxPz{)Y!EfP7&KNW*Ah>&_?(Z?6#JA?G1o9oLI1U5v0E|Tpl00e`7uzID}g z(6?CD{t%z00d|ZqgPIwfi%b$4l6i-T7p108&~BSBi`Cc-dj6%t`&P~C-Q zdO>&HglOuZxrua@>P{R;SsV%8klQ_W4a_Z1&=(O36-xg ziw(ZYYequBtVHxE_UNk7A!V^Y=}3L3oRJb|aR6sK+Zw$>^ciW#&&C<40}AF22qDip zz32dhu;NKO?z~2z%8919w!(MrL~i~%t+X1wEK0oE(P&D({a|kuVvt$!-VCJ5GO%so3mSou=h|^#48%?E{NX zASOt9f(8!tD6^&epCYD`;BMd-?L3GCQf#s7kejfaU*Ux^u=z_`9;wd`eL=Q+@m}@rtEb0@MVlW z`$G3FS5DBm(trXobm-vub=kmktc{YLzeY^JK=x+RZWzjS_daROwhT@*?WDMubt5Jq zMUJg;^B9KMU1Cp=9KQ&qX_3}W6Y=0YmrC!8!}%6jJw~pHFH)s-oTt{YL3KLLeAyfI z^Uim}Rpc6(`n=^w`a0(O*0{ny7~7(F4SpiAOkaF0+gp0($ByVbUx$pCBHv*@Ecw3Z zY~kuyJ)4Fu7Lr05Z}jSD-V8cFb99|N1=R+zI~@iousAbr;i#CGVHxu1Fli^q&CT|| z*BT-fLOg0`*q6#1j4CxZMxN;8eLlM@O)R1?jdA+rQbOIc#omn4apXKF5nFpu)@a`% zd`~v<;PlX)zj(n$E1EHkU>sk_z2n>L&aJL`JbXJt-B2}z$PH*FUcOET(W^`iR{UQSIQ z4oP6bBmP9ifZUMb1wSa#^S{TKwjS&jx(Dgql@GaJ+YdvC#z1|2(&|1trO5>JoPhc> z1FB6@@Xt>quW}^(-48r2(_5w(Vr3EXS{d$S#)d6^&U1>29KOGgo35NJH`+{{KqO}b zH5x%D5Z>rd9lM*GB#|wML-81&v(Rc6TTXem$vS9A!veZFTb8ipEyQ3+9Lk&m8$a1^ zJ@fMpU5wzSZKS{JkdwVH)$&Njz+2oQ$Kv$HMGAQcEDYOeeXVxJ=#IE)7`w<+8Z@ov zD51o?dL!YJdv$S@>0xT!m$QVDn)AxZs%kPTA2u~0oeth046x0|T@kMzEXKj97Nh&&h$7t$)tT$E%JMQMpws_@GCVBO} zK2|5s8MiC5SIrH6%xKSGdzQQ!;V*}ddi_{#juY)}F&?k#+Jiy=!E!M+u_u^9Vz1kf z%ugP^2(e4`#(>K|=f4m?+ zlf13kp?AijKzyQcf*pEK{H-er%Ty8!366s)ni9fBUTctRIEb8HmV#w0EC%^b@9dH; zo%DK?CQzn=1Y<=sj!ktx7>Nd_GmP2+#r9P1qQL(~8|C>-@0=5HUa2CzmoczG7y^dL zS;7$Ti{fqpXL?E3G?JL8GI((?&^H%};DCmclQjH+hV!iv*T8cc?+9=q9YFIT6Phv$ zf>t@e3|y+0VF>iMDwI1?&@-LOVNT@BDU6Wz5=^Amv|`u3r;wZ^knA9-ByT~@lVWMs z@i!C5imz%tNje1H@Ma*H+!YhLgRQG1_lVA}@;Z=E_`|R3=*Je=b7=>{19EqM#_)2* zZ9#OT*ToYDA=m%9Zos}2q?%s8f;!P$U?+KQa=r68g_!P*+3Pwsu?6D&3UJLdrQBhz z>LZA`=X4c#cWM4hwPq;vw|`yk%2cI6MbKZIN#o`h-}+9&qfyD-6v1NfUraGVK@@+F zf43BlDr6QhCcpe>@%z@4r2nM1*YLuq_c7nP?Afc;jL6M>$ESH@{CthURB-2T)i%f( zJlP7SyGGx*z3ziPlxC~cl*hEY`Q;)3K#17CXGnZTXrI)G8B7F-{2t98o?8;3uT@#C zycUFUm2TVHPs{7f6kixsZekxboG`;!d=(l+ryU#szV>)oAmjAxgMh~4dJ9c~qT-*s z_BPIcBJNg}-)=?On{B$Us0O9!d`L!4?PcT*NfN`cOfIuv&lx8Qf6fMfeUT6w$_cT3 zIb*%O;cBZD-}kV?;%)H~&^>6}<``*Zl%sC}uRpQj`n`-eE56zJbr2IG9Jr(T-=xW4OP#zbIeyZ0)lZnR)`vN)wSzaV7&J_ka5WrFTt;;&y_@FuL8 zaFnRReZ4xofsYZCcK2A}Q%KV0AVJDjXSC&I#SS38EL9nRy7z!Y1-64QWwIq{|m z{ixf;x7;%p&=5g9FDV=-nuNla9oePS8G@cSf@<-@h)26>K za*my#t?}`2@Rs==HTL!v#tZp~p_2H%H44SnYbnHU3pTXVK@!(kHp(iuJzidD80zlH zU$Sqg?)#M34_Y0V$+Pq?jhKE-8b4ZC?U z17jC|rQrQ>r4&OKlQ;EJ>5qsk|5e@QHMxfvuWM7i?>Ii#{^{}|N)vmPcioQ6;ZGN3 z&sRtEY3+wNmaSl6SCWrO-(Qh3Os>u7KMyvH7Eg`)D~-tKwZYb9&8McByrc8Pke(vj zZ=%$$;r=Thr2X!M_*ylG^1G}6ugl!uPsZ-P$*0~6mT7W3#tJ#`Ri{dtlthTrB_^wL z_VYytcR$l8%$$e^)2(d=w>oGMoJ=8C;(B1RmXDY`+Rqa=jZplV+$WzOzZ& z7)Yr!rrs|;f2Q7VV*5SQ>2t=>kIwb*m>LJKXn~cgG1uJu6YGwYEG?tMRn(>Ry=T{Q z%LX4m?L^VI>G`t!&C%DsvUxX^&^WM(B{uDEWIW8vscdEGg>GJGjjR{z>xYGS@o368{b38|E zZR*H)i&7lilUlcjhGC z-ecTmRc8HQq_-u6D1W0((P2mNH7oiL-2+n2!lkzwkZ@y{jSgeMOA@ZN8+Uk1pHtSf z#x+f|KEjKuu03wz(49V6sLHgfW>0m`ugMr0d77w*Kx89WO-xl6$jgKZWQ$g) zi)E8!Wb*!u6wIr$GYzE-cCj+1d8Yj^|D7`S+VuKBA$Cdpy*^sH@h#Xbo*t^_uSp(n zn8+UBwS8Oqcx|rD%_I1Ag|6Q*dxB$PI!bbXTx$Qt`?)2Uj>!9Lw4}*3Rc!is(syZF z(>7lv>A+*pFIC2e@ifI>AbhV~lo&+S90ENQSMuhmal4pb@shVSSNWXx_B$U><9AQ z2#CIw;Z|Pui+bBAh=uQOR*p75rlQ~8MbqPAUQNrlr$dHnZKEz1)le7bXo)GzxOTlC zTq+rM5G2?)zk_WN$^R~?SuE4CX4R#j^2$Z>4~un$;&E(}xZ(9bOoun_9U~W2Z1+O? z?0+xuj+#1+)VI1C?+2t1Uk^y$HMn;yoZKi#p-XsPsiv8Q?ER5Qtv@Crct4xg$HXD` zAg^0f7h$by?eX@9e#&ot2tq6y3NvwdWYgpKlO^D|r((Dynifg=r*Q0sa+>t2iq-eq zB9?-OMJDUDMd6uUqCN*Gxsy)pN#i)sZ~rzpD%y0s<8S%-S(QZCFR+DT_K@` zSFDt)F1OLPa>UEnLckLL63-9r!GDBrPQ@_ZIewxwE0U zCd!YnbAHq*x$0?uVqdQW{dgM*`RPZ^67XuDRZQmDh<)eMVBtvY>(8Yx!OeJZ^O)EF z-}Tu0ZZF|B3l9>+NbK)QLlM4g=SmZe+=5+PQvB-`&mSKnA&Y(jDSFA^k(yNrOAj@C zLGF6cB3DrTTTLB=h5gD!2w^0#hU?8VFQ1LH@z1=A%9UyqyD ztt|J#{`Pe4SqSiCL#J?-g?ZWVO*=6e@lhZgHlGU89OJW=lwM-rc0EMix5!8o9cr54 z+wqHZXhZ2h#}2GCZS`}>sD@6!!^0XeV~&u}kGMwaiz8_`Q2bw- zildwKa=kn+SBU9yVRmt^$@Uo_84hH%ek6w9$0cZgJ`k5l8(r;N@|c5_N(||r)E;Xp zB$tXC?_Lp<9sINNcVM`9xw9+IRKLgB%4RQ7At&pY7{J#IqFvgb2PSz5V zEhI~pktCEQ*<~mpvM-TshD7!fO4g7)`x;|kB1_g#cG-6(+wi-i_viche80cXAM?ze z=ef(d=e+K@ckXk~IXR0q`}t;fhmRVThFsl<_<|xHY{LamX@L|Ydh<4^H_DZ(AM^)z zq&}*PeSXcTescGztEKa2@MconR(EiG%x3&E4r`Q~gTof{n^q?`p4h_oC~H zMTb1URe|d$l}Jh&&pA^SD-t(XQky@e;7rK~GvdDZw1F;Ma=zo5p}I^Bwp;6w!>-${ z2NcWJ;XzD8-_J#wjOZaGa4~IF;Hu|nmDN?@e_hEAn;*onnq8VEaXwk?dSR5+(_FX4 z68*{VMl=QW;sd|O4xZ$n-@pHxM*jnf^9uXyk#AN}MTe}>{plFb(yUpX_d(w69MlmR z_j|M+O+>KD^4mNuD_d!{{bN1yYIu#sW-i~or#EX{&nB|qylj6Kb~DT)k!4lFr}R#X z?E#epuDsZ5;X~ly!S^kOg59(^0n_FAdmkLie1&U?i|3|{er;t#ZjPtxWQ$uA7<*&Src_Qg{VZ~6}8mjLZ5Nyo(>VdwoYMeT;w}7T=c1Fe4 z&Ecv6p1y}{1?uVIc?u>>_tj72I`u+bYfN7=E?O8f`z@Y{6bh>w=rd;L&Y$`k-b85W zQDkbF+%q3;p{31(TX9%U1D~{99}NBORLZ7-w-9dAU&Ux;o8RN+;rnc-fJMlcv-!tg z!%_Ld5}SGu)m(q)nX5=}eZ=GpD2O?XGdgJftehJt4|`*w5F4V(Wc1u!i0Pd7|GD6< z#GN}jW~hvfSBo=xo|~PntHLR)T(+ZUyLeCCp)6LoY$x$SgdVry8D2rPgSRHnSNuGo znfPdad+9wvy*EBQ4%6XBZ#1HPx8m>C1ZYW+R6rBwU3kMDG=8e9u|oyt~!?x1zGIEAOCbyIwqMA>PP!h_{5&q(xCe56*b zK7kkPJsO|xhd)76j-@EjlIlZNM^1{O1aZV@gT zCxgu;c?^0^J{?Biv8tK9zU>^xks?v4@gDZ#4aLH_e!Jkr(hTTZ07tIfd^x?~p99B{ z{G?#^iLBO4BPgz~(f^7-MxdK{T}|IHy~C`oOt{Pz z6^Cu@kBS{VjYVS@^$;r|?IbA%qApf`O`LhF3$gkL`&M=bWQnFARP&3rhU{JCgMq{Tb31tx4>j-4mA` z#BK`oaG2!cw#*P~gVX0cf#|8(*k;9bP+B#2*kU&yuiZJE1BY@}F1gG9C`qf9c~)%h zIgy#i1}x0kYlAqM)(2G;;H$}Z9z>caqA@mHenm4qHY-&+!t+akONDG#T$b7NU`CKH zc}kL7Ew9nk-wPYv4TSclz{x}3ThJ;WfnSTLoS6Bf8kD@;OfRw%o-f7`i})R?6p+`N z#IRkg5nmYMS=UJ+F1d>rK2)a#C+JoTZy-7;wp(cYsb3WKJ+Jd_$Mzo6svvDvM5wu3 z?E_+_n{t}GorqFrnDcI;hKeXfTdv)0eyA+48`Ki=-ixuw z;JX@;WukkO2(?tllmw*c`DT-O*7fJ$P}s}iE`xGw36c^G?!t5ahCx)`l`GLq49(81 zns$oLo~^yK>};h$0Rg!FhUY*CD!e>?nB$C>JA&rHut-jg!^@^mP#OL0K`X3olUx@D zh)$-XZZ4lkc^{0}17*(&vuikUq@WsM8bwVO>!a#_r5OnCce3vDC*ES;leV5n7q?cx zuMkK_>h$tBBGY`pYxUmaw-K8Z8i~A#YF^UqB&8xz$h3f9HaTJyd6-;NB{-DCo$N1+ z;pVi{>S2N-GyE`AyRlYt%22YZCh?bH;K_yjj&DExKB#8~1bZJ5uhq@6sHhozdM^_3 zmO>^dP=Aqr>8j%>3LJYzDZTLit_VI2ie6Gu)_&IwKA9UgtlpEBij{G=l2O|8#Zgmq z$i-J`RdZ9#tLjBm+0PP3g!2$~H82yQMP2H1*kpbQ-hbwk72BglqVLX*VHvy0fnVZ_ zs9dM^DKAFNvA`dUy5l=prP-c>v+MylN5@_cYX3MtJX@XO-Bs%L$Tbo# z115R4hTWG54l6|D6cR3oS*PA_Gfdo6P;I&a-dq)GDUC@*R) zq-%AGjzpCN;lS+!1E0vhjI3rATeSc{qKyhH!2)&zDF{ zcQM6DD)>%nv}gU!3Z$Qh+DD+8f!fLTkZ9KO_&KTutu(nG*;sO^&>z<_WyT#85PqP9 zyRDs>+-)wBxKDR^Q8Z7$n(m6*k{+$fH#hkvN={DnK(?_e+YQ^tc2n8s3=9va9Sd3F zeA49E;3wxUSZ%BBJ2Lk<&Ga|82E13Nt|}YAwZNqUt<{1>psCADc51zpG%}T!el1 zmMoOi))`gvN^%Aqli@Tr1*g*UU4Kfk-L%XV^8(}k2)#&~$!wHvOq)q-l>TJ{2Wfw- zBg|xaF^o;Jvcu%sobz0l74Fy=@xaTfmgt5bv;4tKLuTRt;N0nT@)U3JxLZDU%X zRr>7Q*VikTZ_*ej!+@K;i;IzQ1P5c}_LLP~+w9z<3_1ZPm;HLN4dZ){CV);jK_Bl^ zO0O*xNG&fe6fl}LQ@Z82f*4y%?YimJX?0>r@+lBtAPO^Kx!0oi>yZYI_*t_bz-&a8pbPE;BJ!& zd;2h2sm01|6%R1y9iiXf>zPdMMMT8oQJ|TLfU4 zsC>QKFC66d^R$_69`sL18GCP;V>AxhQ#aeJMd``}OitsFk>$?hlSwHm#+4>>a*A!` z8|O{4iQ8TuZbDtzB#H539XzsG#BJeCn~I4vo5l)$u#I6?2^Wf-m?{VP3=hq#oW`DX zA_Gqdjb(OAnW7_>d7og`LYp>kQ4>nFE;l;OIGX&iNe4iwRy9gFobKNQr-?U-KjkKD z+rWU#9Ut)5p-tvy4FnufHT6j{5|fTK_uvt%oW*Ummqn0mPSF(XlGJT*x3^yNj^0{? zCm4JUunXAJLR`8cz^G~M@PmaKh$=s5kvLrUYoflbmv3+}{HfRUR%H=ttDnt77Ugw!8IQJDKb z=?o&4TkkM0u>BbHz|%~@0+l5)@pX)ni&ma39Kr&s1Bc~+p*)xmX`)AF~pyBR5c6`>0RNLs$vSP5u_{O4y)U6dh#?ls){6xI|XVkcZU*$&h^cZhdN`+L^#nw zc~`jRX@MpGMij&N`5A(peZavfQ;wkI$yDvqB<)s6I$$c5m8)zl1dcRuD!ER#E>pHI zhW2ni#{ttWQ|z{AbHei}q3Pb#;~!f%C>2=`>JbuHo=2k3NkhfPFK{eOc1s*YA{(yH zI178pCM7yaBUWgJQ+z|ET-#O5PMUSGOl$|XFL|%n?HAfyECenbNQF^5-`(f93ptGs zMe6&#&c6v%ua|;0D&@Vd2S(i!7lac#mM>_{U-uGXIh^aNTnRrl$PQLVROtEFO1rj$ z>gA`Qnj-bcMBt9#4%Gk<8u8sz+;+mLD%$J(rbLHnb0m}Rb%XapW+Iv53GM_+FVP-R zlzx0NH3wtP;&90~d46he&4JM~txz_RX;m?cL~ibZlOV{%br|zmUj!fU_W+7 zz;GZ5Ngigs0(@CONSD}f&a)3w00xN@S*$(lAMHW*dr=`z{9E+2C({_QTXn$R(YQHP zJguEM0$7A5n-Zl8nLrQBe>Wou7WIV=G(psO#S~!Cx=)~ld;$JJoY;-y2S6|88bL2Y z4~xn~K+kt6y8|)jL@++&r)0uGA;1@fa@;)S6)k-6_E+RzVMr7vW0DXaBpg@tLZ=vY z{rE%3KM8#xqB77UREot22!C8J=b#Dg+zfg}RqE!#x4yj=~UNR0^a@KnZao4WPyi zUN`SA^eP1y8U5_S;7n8BtYa9jLEq7XIj%BG;p=3Qk&U=CYjA zsoy7X5L-*+Z@5Q5D=>nGj|yu2`eKQBbK0rHe|Hy4iaSJloy~IE9gvy>hX1%3*gGOH zE_mw+L>w}1BCR+V3Xw5|lBawX)SCHsPaTN!bIs*&m}#x*7;dd_>7@o_H~^@!Y%dEz z6q$d!qnN@^b`^&C#MRr3WR-@xVas~{Q3F}q^`LnI#^ih1NQ_Lm$Vy2tv0S|mF-zJa zjF&wBC2<7;n6U)G!U3^8T)pCCJ5*4)0OT?y31*0^_Z&$_D($W8!CQAuM={Nh$b28d z$R@xC4ADcR@mXt1sGVYp5Lsy*@~kxix~iBWM^*~Q;JA7@NJawy0~S0Wl=mY)hKGfn zD#ZUlsl!>4=P-<^Vv7FR1%`ho%MhjcFIRuSM;enPb5cyw|4Qz{{E@b2t1)yqgS>i?fCfMIYfHf-`4=ckle7N($LE<}n*~s^C^H4^b+USigB%J0|B0EN zOIpsgzN@=1+7VtWFX4p=aG5ctsajG@Ibf2VVv60#le96i z&OKcJtqPjhZKTy9;WuSXUHx(Yt@T4q*4lqyWWz1!=+{3QS=kFSpX?~tR>*VnbTYzGYP=LW zFNS^KDE5^6HBj^p`ai?Wo_r^ZWeor4%k4jZoiR2_Ki^lFp?o@uRv#9btB2P2F#Ql= zWGyefG)plM%LdjT6}fY*Psp!uEQAqe*yYFqPPUj+B{;SlDMkKl?)2+VOGa7*Bi}4` z52zt@zAG5LKMJf{zE%lL2rsrkszs>Tw^rYbe;3r+Jcb_f)Jma-hc3)`|6(G{;Dou# z@%1(Q&6LR*7%umFkVqfz#OGU_*xOz`4wId$@HkrqGd$m;GG=>B9qy3Wk!`)>AC>MP zd1JEIbW94PiYY>qhx7+6tm^3(CR_m{G{;#DiaV45G%{lH*Swzl5aq37{f zA7Rkn#7)O|f2mq5UbQ*R5t)y#eQ=j$XVs>a$ISc7gH?+1+y@QrWr0nKW385K2eYS( ztS#ta_^UX(N8;FD?>t7XU7Hy@{i(hK|7Su<25Nh4pZL@CnGPFmz0hi-Zwrc z^QYwxpojPdsVd;+`TU-#k2@Z zi3q4q{naD%$GWrq6pjf$UWsR%QVyC=GHEcIEO%arj;!aX*{oKQnFMPhQ7{*_&b5Di z0Z<3lbQRr7rSSB;*hF9)Cix_pj&Fb6e?R>nm_1r6$A`}~RmZy-( zub#eA@qv&L8Li(fW}MLJbUNX;^6oobkX4-&NqgwZW=Z{q{S;96Cy;r!*F*%yzpjOi3c%(z4evC9-^WImsvkTCsJzb@7FXK z2w84Oc^}+qqr9}Iu;4n+0V#DaQ!g-25GKfnIZgx~%T@X^`aZ2B zEXGdI+ts(jfP}e44;Z{7xSrKFOHEKR55vUiW@HQXAKLetnjLqt$`%zP8(!(OC&m%A zliwFJl&1fxZia|xV`0!>G$lSH)xeXoZ-Gj1nde}JddFl(4U2u~h|PSTF=E;7tc7A{{+Dcx z9bE*^YCk{_2e31vQq${a(Yrr{^;*oNYbyslg+{cMi<@^7dtKz!u2W`{@$w~pJPWbk z2VdDCZ+=z&(VCR4jRHZn!Vgf$qEq%czhBj;n~QXmNZLUbn_L}_LRi*mTgPA1qc3zWGHJhdz*#Szf|>*;-l;Qcdi zU&CtzIqEf~vThV^k#?&~E|An{TitH+v_J^scc2}H-*TPEfnPGYN)9OfbpvTy!>s<} zM7|}i*FYU;YWQ5FJ$iOA5$kteM_|&5;;pKIC;7@(rzaWqQ@Xft8NzPD081g_RIe1n zD@Ov~I}0oSTxkruDXR^C58NnvIWJWNYvNfiDO29|=8vQ})7`^KAAN{^|4~88KIi36 z0noJcE~P8*8zSI&_{`Mf7q1T(#eA9av@A|Prw^caF0R=NKY4Y`$9 zgW2voBLi_5VCi!LTR5e1^}hz*;af78oqpa8A32#GP$UiK4tOQAsxV1NM%J4iL3GSu z-cE6*T8bB!uypq~oqbk2)!=R&HFTHG?WVQ+1)HZIH#T%o4n>U`6I z!})#W(r!l0sUYyY&qev!0mW(JW60p~bgtO_LeEW=8N_GXhw_Zuw-PUa@tS#q4+Ocu zUl*yd?{^Cad$v~f@{^T8LPQ5~;mE{$Xlmb1T30|!>3+zW} zG*cE_Cy>3vq@}jcO94aEPDrH^X?^3o0%LUC2&Ypc!MAv3;MQxjGZ(8*h2buNN=6yTiJF#egY{C0z@!%!J^w-Iwwic`C;>a5C#0^;bwqz|mNW(w zJkwfPn6zc1zkU$Fw?q%rn#^s-dWYT{^%dfml6RQqV@O@0<5}YAbX(avk#bo?)qJAE zwKr-%#W4NQ$!wSZj_r{DZu;H_n1#miTpDoQ$RMMTVE5kNG|pO_?4+W%k=a{=2F6ok z(~zm3{<1VFODBWQd8yQI4)Ko31TvwWa&e$RpY{^@E?gzYg&=GC5~ zP=Fae&?Q}&SGZLaSe9&1?Ab#Ga|wcuqPADuhqBlhd2aAPa*Wu$`OrCfa1=yJ6WCM{ z1QF!@Y0U&;g}W5G%apl+)~J84_!x%}PIN&;+1!;pnXH7x+aI@N>!i?a-26v&Rz2W& z9}C~_OB-&E91&0ye142l(rfW18~JB1^jn*ueWoC^z2Y}Iw-DDG8i~+Exw?jLcuL2l zFe3~2d{dUPyW{v5$5F)MP2lWbqxtnFBN1M`OW|a)y~r}n&uNK%VP&6~1EOK2&9ge` zI^>z=bBK^~i~r2JsV1LJ=VWP8@f*;03~66EW0q=6w*~*EF4fL7pK*Bs$2T7X5DA7a zfqhpu&u!dD{IvWsL55fcGB4!uW{WSbdk?lP(k-s*zotpe|Eu_&oP}YDv9f3qcL})$vY$>q6oUALkN4JssG!2c7inJRs&U z|KKK*Szp5^T5*$l@rCzNX#w&9=3b_lzn?SN!ceoUcIIuq?|(H;f_=`FspJHGo9S@e ztYyLStpxtT+K{Tp?zh{-$Remd5iOv2f8|CgOjeOt1yX`aRRJ*h`^c*nfzS2J_+X}` z(=oqT{dJ|`2^hp1@H-d=`W~q&`fU!bI77)@KO8W4#YVt-z0hl=0~(bJ>Ht}x_B>X8K$;KuCnewgiHfcb)G-@x0&Kbu zUdV&R0d76vy_h|2i*Gu%{X77W)p$oKRAX+Gg-A8AIw+w==;w-C2B7hO7!{*^x+qd7f+G70{Gy|))av#U9v5=^IG6K}s5ve|@ zUKLd-j<>(Yr_fEe8VZENz6rlG#o%}?{X&JO3b_U|ub;-h?< zg>;l;V)-`_3nxd9UE)|5l3l^X!L8igbzJ58&9AE z!q15!!1%&^$p&1K-AR3xP|xdGOA&l9cf$;&#U5m$ZYQh;eODtvNe!X$?E`@mw-rEC z@4y4)gY%%s?>jvMoCz}1lGf+jlMJ0NcV}Mu9c$_Fe!6c#9n4P_T71WuyJF1m+1jb#mzUJ@_7G!U<8gp=W`Ul<#M@%LrOWXqY>y8xTfuy^#4`U!#ZZ#9r#(qI}!4xsTEOHU3NQ zFT_1OCQ6?s&4u|7hmI>qy_x>GGIvkRYoiu_=JyN^%JxAdA!OSr)Ki?@aa4bC7;Sv3 z{mH~}qVcEpKRCp@h=LT4Pwm-KoHf%d=MgJH+yKuNcV2g5y@~EQ1++{b#wFD;O{REG zMO|>PJdyeXl-3C5p3KXY?pZ(jr;?oXF4x^{wd#xJ@QW^G(Jd;vRCxq6HAEkgyN6 zpKV{ApBO7hc@F(Oa{|x20(>Y7K~^U`fko&`*+S5v7FOr}_ujBfKFi}U4Na=k=Sl|> zld(iwV0W|}=o{bqqt06ilIVP2eNmMG!oQ!5hccS3evbV#nix5klnG-)%o(Lb2Uszm z+41=CBp8$gU&rdm7zLV&=BV{m{U}*cu&|;hb-?(5=0&&vu(+r*1DA#( zBn4a2yM5x5fDaL;E_0*aVd1yGpJ5E0FR+=>^!b=~?fB|A!JV-F%rpAI*7tx8;IJ2< zpu^Py*2K-1c>$1aZ07JGV%S+M-5+1=auCfWd)N(GbiJ=l;YN_U@k_b;$p6^WQm^gZT@%xi?qM6#n4<^OMK>2mc=* zhEnjW;=gS>%X6n$-T1KVdm}0CHO-#;jU+-pF@uIPfL7N>}oROvXLcbzk32A(a=kiLb}I+iCLFHR*;^A42_n zo>qTekG%Uo>1gDzYRdBA_LMZI;&VUy8tNGOy=CGpiCP6NRFE3Q(Jo|P=Kj6O0z~`p z&arxZLZq4JNuU|~Wx);H<(irk`eJGZ#=T|D^UGT1n7U*7V&U~VSbir_-|?}r);W-C zc*4I=%lx&btEu0Pno=NaO<^cwv#mP4u9b7SZ|~8AbAEa1Gm&mzDL>KG9CurfMacX? z8MMSblZWQM4XZm&vH9R3Ci;l}TfXtd?P#k)EFdAq_w|KN$7J5IZy`}R-9TxTG zC{bW*E-})~W0#^vb8k8H>JyLMB@V6=pnXhvk9NL7_LaOqZF^nuU+BD z+9RumQyd>WuUx(K`zN&4A<%Dn04_pzeBfJSNB{dLZB$ba%~LSX%1;$glEiiq*E*PY zZg9vFYRRA>r1rMD9nMlmqn?HH*%>gb5&CLn> zL+gV!8-*Z-$6)OSr0zdRceVm;?sgIAJ5wv~kB>0C`{Fa6jXs`P zfHuF$!#G`Igo1^z&@$wG=3=aPyAT1GDjJHdOBk!S0x^v<_x^#7m#=y|BCc5Gf ztryEuATIyC2>Bu61E>yGft;}&xTPhh1(?rh?d_KXSJPkz&Q#PDG(lIfHQSWUe^NbC z0N;dqE7wBdYGD2JxRGUnoOy~o7SfE;vwQsDkIOh8n<9u8>yVj0DrU+n8!frWw-g|< z?oE;lKR1l;r9A^$DiGI5JgD$|zacvni_WWz1x%eK4_@}Qb^`xBhzL*9XS#vN36Ilf zgh64wS98w*PuUWP8Vah9{#MNz7x{IEFb>i7y1FJ!UvYbmWRZGMR$%55t`DewG_}D~+sxQrZ zd(JZj=XQr&o$+a84V=j-xzGicx5m&M&^{&oPiT0kNTk^=vH`GtNe@8}WDeY7S=HPD z09Bo!B?6M+EUe4ZIK;?9mGe13yi2!Mr5~Q6zFyboJsasF$I|O+wrL8HI&^~yn}EW1 zmw_6R0V^Uwg>3@L0JIP27g;}oI+YcE1`(7epiOI* zU^B{J2M#(4i=XQW1Ej`1^9_)iBnz=3`6)&lH&_T_tvtS{+9n3=4ol#QL>DCyDPDUZ zq$3YyL=)(2JEc=?JROw{VwkG)X?lqHIH@D?9{?Kpc!iS(>4@yf+xxNs;4UT-^n(yY z^+24!&%NUq!7Vzj*&N&EoEtMh-t+Q1gi$nm_uy3<3>6mCO_p~cZcm=&_W>+1g+!fJ z1qnV-kg*y`r=DcERQ(#^z3!Ude7^s1hjqV418rJgCkuHv*kk24IZ{a-&c=*tHwG{E zauh%qVWFqGN#LJP10a+4C$wwJsuc)Drt^9oGGV2QjtGo)+Otdmp4*L*6z?KJZ!E`LMqxu~a%Eu%lbZ3We z;7HIw&A2Xw>F}v$-jO1fgK*12ZuBdiOyGz}^lBE*W9<-=HQJS`;IyDyJeU{=K_#;g z_Lx?FW0fe`UMc;r#6ZVH6Ka41?bVm;#5+k6c0@vE7>V^eb4Z?dyT8QF_6OJ?{4aH6 zn~-WfL+r>04W&&~!Q#_(0EAU#{-u8C)N})ktsMadCfv-t|nTa;7>@6W8UIa1qddxRFGUG>>~(4oNR_@E~?j) zACR3b=J3E{^k-yf@0k+-q%<0XIqMGgRfM+piWEgo!7Nc40Ui^AL*3!V!;O58})-vUSE1@U52BbBa# z1n6lrT=V~pxwZCzF#HMl>!(}xS+8Lpel6x1?ekH7@1ZO>s-*mukLlc$w4a%Eq`IUW2Fk3$Y;6jE0T!x$oGikA_ zASl9Uo&KC1!=SB@LKez##$dv0ez5XI?qV>AmzpR$_3GfR45g>VwFg}H+Ty=f!OrnN zknO#A2sENl=27Yj9ymS8{%TZO$8u#5M4OM{-dRD96|L^wgB;ES9a; zftKhiLPoC$osQ}tOcD`|O}uUXO8CU2p2e$lXFfv#ulexCy99yj5Tf?_g1O`>(8l3Z zOS zo89tf1h^Xzbod);F!nDDnDobw=GLg!eZrCC1>axiv9&>6eP8+s2f`4;GBsj--)}rM zkpa2!V8P}{na8+5fN4B$&`_&s%){bqxNZ;>TjQGWAqxls&46m>c3UPhpvm62Gc;}* zXjP2EbpvHA185}ZyWj3Zy{0kN#V>KKB4?{g@7;Ehvb8SgyK1Q9H6J6ZUU!eL6HAyS zPLT~}5xD$UnGvS0!Mx~NwVxJ=E4tIX%ND4se~L8-%`y;1^Q(i&(9Pbe`os&%7J!?z zIv`t}gmWEPMe?B8w$zVv4eMw>O7KZn*Yo$eCRX}RBWT*2*!NAu((Ihj?#EV|)=skq@$-^^m zCCc`lZ+>R!Zs;rB`eb}$qDE-YCnc32wK9n`*`y8fVdIrLv3#VtG+#w0p5@%f2 zBapSBmdJY;gS(3h9f8*&(}Od#hQ7^qN>W&iI9Q7Om69rsYKJ6j18$MZ|7tgLxRe`z z&~J-kc} zw82dW)`So`fqIhsxUN@T?9Hu0kCfLRTcCIdi;ayL`)`TYszc`Eo*;6uVRpYoT^}|? zq%<;2khqX~JARYRI}Hi=6>rgl^dox7hqTVwt7A*eE>@6Q{jkC`T2aUev1Pv-PL0*( zy7!$W{<6oid_lH$G4D<{SVOUs9uSM&Z?M0+AmO4`#|yj<@Yvp12Q^x}D2uSX|24kX zisARrQ&aIwi>}`iE;4&+?BERmfN1@G9EIbuSV8tc15GN$hDBM#rVXwyyR=|BP*&Y4 zLn9Wic6jkH5H}6x$h;S-3S7zz14Hbl;)NE!rX+~fhuGzF9X0g6NM-sk0CBKnOIDuA zch|-w8=R0lZF2Zt59e@_rLB!1Sl_(sAJrCV=A z{AWF3jVM?v(B|RyX>ZZPvsh%hzeGn4g+$SdOS}kRq>2kJWS>0l zJSZ;65oC(x9{sX1Tk23;aQgXAUeV(K2fn~up*KIOIhpUK-tBxNVZaGQ=o;cQQZSG` z32!7PT%hs`09*;%e{HV)z^5z z8^$CH6n}z64_X6#`+lU5-qS1&I`hySZ%ID*mcx}C~Ig>`d}{+ZJ};%C-%Pul-iljVTk* zR8S&G&#^JLOMUS5%`X4&Z;Qp@Oxf(`#qM+Q92@|>H{0WtjAM9~%|GXNyxF~?nUjv$u zd+o+LY-+(zRX~=O?XfWtgcrVQkGX7*lXlJ7E);3i*@!vdLz};q{Vw5Mw=d-AgC$;M zy!}?js{%QFT=R_bqFSOgweQ$`VS#0GBUTFGMQWf))C5k>KERe)8SETjSJ#LnDKB1* z_bEY6`?8;@cQ_QZ$MNrkZWmf{&5>LQgYwGwtxMLyKKDH1>Vk$fK6YR<`x*Pk5^`&b z_6#;4$VW7((N_Xbj|=1Lcw@a$J8@KACqJ=yKa{WIb*+)Clm2*a3!Sd&&L8Zm%Mbk$ zO>-YSo4`r0#Rl|T6(ny~Rc6o=EYD=h7RV;}64UDoi^ z&$#xFy+`t?Zbn!6%)K}I`0|H7Y1ZBHkCK1v03&*jiGs~u8JVL%B@7?G;_Uq3#pB5$ zUo0>;vkUnE*8=4XxwL6-My`&@2Wn*C`QoilkZ(!6iiy%L23w9u))A}zneVhOLOz_G zV4}Qu>)3~ojPxC;@S*S`E%y1SjYBlI*NHF&lOD8m-w&7x9YI;87@HS1BPa5&^pfiT zq99&%VZaxN9BNjmnZLqdjWbg0@~d?zMu&2Y8w#7z6ZvKaIc#+bD0jUp)?+#LuqImK z=M{~U*Qb{BXaq%5p5>C2p3$ucF?$86F*Z*j*Wby-uC5~w@cYFj{-DAy#tmSR=9Wtz z%Ut;dRI8C=tej$(R;^1Cuiz#?70aCUg@JUGV`(_uXza%6jTCT0ovN2W8I++m%i)j_ zH@IT_gK0*3MStGB?w#+RpC5K5@5)~H8&Jz=Q#PnDKLI4wsoH;ChIYLq+^r(qP5kh@ zj3}Dz#OK%*P)V#)5rVufgu46k#$5w2YSbLB%E#>fD`Q>Dupb3KTdNdB9RWp@zzs>%(BllgX9(KZlq z5~jIwM4$3vYZul?R@5FhI6X3Zvb%r0#uLid+Lg5OqNQ>qx~0nJc`@>#_sUo3tM^Wt zjJLx+0J)#zjD~&O6#WJ80MJP8h!Tu)hcerkn5^7+VqzkAHry{1MoTx2HQ+kO#6)>> zX43ANbj>Q`Z8;+DVDgo6V=mESyDu=>nThhT_aRR9V?y=J>*c>dO@wE2CHUiDC;Mea z#V{A+3gp8Z)oy^us(etY+YV%MvRejZ*cw-e$0!(Akm*T#_+Z7W47o(3>^{I~iAT!+ zS8~~EOwixsbFCetB9--+W#*-ZhL)A}uxpf;EZ-QHpdU_-#>41|hld1Hc1GfrX^DS; z6mxyKB}B;4QNl~Q@*O^L#N7wWp3B&dxLw$$^E+J0cgRwlHRvch$Wb(@$yXwg4M!M< zDC-mHG59=(cx`wFa>(=wt2!p2Vsg-fK{Ao8hFR0;=b-DoAL3DGe8r;f`kE#?oD{dv zMLJ66l0B!o{Q=-2@^Q)Qy0#Z~gzH|(NPB^u-ZcTD(yk6m7`ptf^z}myP$rM4xBSwW z&1UmniJ^J#qso!IE4{*(J^-A;9+yCoUlO0{5$#YYFqI$_)_aEQDxfKCD>mT4hdA)S z8WAmIC@cD9G*5ovoLioJug1j>gYepV&DOLAULecsxx}SY?UrunYm6g#YzyZs^4NMw zFN{6hw2mqgOKqEdQX-~!82D6=*EL*7kJtJn((*9|w{NY-Zk-%yIj9)<^l4y8L|LZ3+#fIqY3(2=0zx}Ic9U8TsNh^`#s)py_XWEMo>HhSzk%s1k7 z;iD{RhHkx&E=EL~U5OQGLayFy+F z-HNbO_IHo5*1Pyr)cL9aZK-ebk?Ivib}*USZAgLUP=>{@&Ctb6q+K?b{)UmzJvLC7mb zNPg-jB;OJ;Ft%EijOOM0eHbJg-ArK!qa(#Tv()<8xSeEh&8m~(^x?2pAUA$>koB7cIOG(Yy+pm5%DV(j?E(v-%B=7%S z(F2$SVOe^nKY3aG)~ngeH1tUsKluLsS)B(ZX0VdY@a}wf?{fgq+l;sGK;j4Cqfhfd z7Cqee?rR_Q zn`T^UfhS)me;SAxt(O#_8d#Gv${v?xtCC0Y2_B@IR>RFmw6pbd&knqlJNIF{=B1BB zQj`1#ESoucjGXIi#`CMRwVqPWfR@5}_^51}QD6%1{P%C$)5$g&nDW@*&k2Vu#v6m4 zQpNERY^tqIuMf9(*CnqjYr7EY6s-aI??ONB&UJxlaE9OGn$Lcdl=zW(g5mMc)fnJT%TP*3ucZi3yDVs+TJ`aJAeBlbTdu^=L3JB7cyA%FQfa&n#=orVOeS zelM4S*~ZZ*zkAP@E(au4jPmQXSijrAURDProw-+1?WrelWuU9X2&27B&c)pFq^#pc zt!S<;;}XUYv|`%{S~0puXvLt|=y9t28>{OBdA?DDmnnCj<>?t_$vpxE?3mLc$0C>4 ztQmshx%C8O2eu21rczA<2PqJ{m-F;IyfsmHJWYmGQHYWdQ78k4%`g4D4XsS9 zZ}|ofKJq?oS)$$56VQvxiZLpxzA4lnGRR|tad$^{Et%-$3GW{dg90)UgaZ6|28thF znzvBUN)QR?MP;Sx7x~}xsSaVdXj8{DRY8GWWkP`y zZBPIo7Q%4ZM(2K>miAcc!$7sK$)HptwT-z=3v+t|361wh`RHxDlSP%S4H&9Fy=53)kC&9_u7^x7SiN&VFz{eX(mM3K<#sRO9N6m{F^Q2Rgj$A(PcClTLJR6}{u6El8R;+tt zi8Csn_FgjjX-rh?3BtR=Lsp&ksBNRMBi9LFkp1krdzo9BH zQNo~RG?kRE)#r}`oTI55O0yD3DP22y8WvynS=JbmF*iN3bSa;@M^EJmgGafHT5fI0 zptjLqVtzzl-zq2qw)TW)9Hi{UINIbC+Dxt(`(*5XA6XK}59rpjd{Vm|kU9BHy}>Eh zi9K$!js24Genv~nNJl}Lck$o@BSG(S?JgrKw&ZnU1NKj4wIfS~63ImH;OiRY&^xvc zB@AXpHXrhd`rnK`h5e|G;F>~H$xWfZda!LAcyqQ>lt{3nQ*S9p=U?a#8G!dReAuLC zd*6HT{b@}6&@0Q__EL_k#>$`Ek`bS0E}S>&_H;*e8GU}Af3JV(8_1$(Zkan!Ix1$| z@TnyRG2Vzh*JIvO)BNb*X`sod1`%ctIGA#)w~(a<3mJq3?Emrf9`IDY|NppEhluQK zrJN!Z*}HccSw;2=iO6<1Hc675BwI2wv-iyA*n1ruj&;*v#kw%!vEW zM%6eP3)eY)*KB;iR&(j&iB7#o;Y(EGq4z8&)_IcWNx|nH9cr89=}y^~`sXAxx0N5h zGB%G**U6DSNv3k{ZOwO8y?p+D9d(h>H!}^koN^SD|#VtzeSXtLJ^H z%WGfG{I-12EhMoC+)_1~AxSrCz|cqSYPhOMkm&H4=<8z&4f6Mp?kg3*o{;7tVUgEv z89@0M9cWGd%Z^uP$i|T_KEiaDz(;RK#aizEZsVjX^cQqxCNd{tA#yJC|GM^6uCf1^ z<0V=Cs323Uoa-9{eqVqLKzG!e(Xcs6!x-BXI6jcut~DX#WJ$tA_$Y>n@NQ{8ZK3~< zE_y=oNPqB#nEvIij_Fx|Upn!K#tYLkbEa_W%jVLmwG3o$PQwA0gIWJye$virBJ&mv z!Z)u3E@Z2So+uoCxC~U3L!(d3JeME89`(Xf`9i+NJJtuwI$l4CC);M;5tr6GoWliX z7Z;(ptVv1Zx2_Xf{;%epF95nem8?kV$3E%m@o8z{_jf*BioRCX<*JI+vU_0v`iyK< zMP`qiGJwJieP68KzX8RiD<>rcmP6dG9$0dmK-pP<>U>P*)l&6&BG)c!rHC#U=pHyc zMW~i&{d73&jM;?w^Q~M^&3#x3TOHkiXe^E}$p<@v;Zq(hE~Y*>rkM6KQK&$BTS2z! zW6WBH?p)OUr1S#$d5<^iXGxo-4dqF>1)i|e=4}!_x2%H#wi7GkO}e(o>s0(^DffN6 z6jMw*KKTmIaZ3d+NkucAp?GQ9X{&FWAu>}{q4oRJcS!epx$pFP+}S)nz4)(S#Z5BP z=Dz5j)Pk}RWQ?Vjir*(`H>=%_OJi|hZNPnw+vsc*?uq7_;0#ld6Bv=HqWPXH^x>1x zf>C_S-2vv0$s4L*IoRtD)g1Rxbnb@BYD(&;>%oOL^pnfJ_I8mB$eTKX(d?oV@wM_s z(X8D4rkJJ^&lUHuhz(Z5>fyOe2Fg!#wCt`v@OYQ&gX*0AYjRfC!f=P8ZfiuekoE&_ z(bj-_%G@cbM37Tc1nr_!<9`1n%6HT&X}6jZ-?I>pIb!m)=lh~7$r^zcNA6g{uZ>+jVN2J%8EH#O#^mvx;|l%$8^hQCxaX+)3kLqr2fR;VxTqL>E^P zEmrD}-?*tG_erurjzgX5W}K*(-J46Kow!fDnl-ge4mshZwG1|C=PJ9T!)gh(zyL%J|X2neB2)w3hzs~ zD)T{5yv@@2TKCf`RduJ7Rdk4NO40|ZK~rcp9A7yO=KN*&so!4mSYZa0%pgazAa#m& z_Mg_?rr_xz8`HN-MU|DK4ZZEA4$~8#EKxsW$FJ7I;CIgDc3HVRNe&I@^y6~?JGQDq~B-fm$MW<5PKX?+K6 zg!mOuj~~`6j5K&rMKDue+*xPQ*goRY<8JBk8jQe@UX>wu&-HL`{N2^5%Sfa6+uH-~ zg+c*NPrWLnPxfpY>-&3XLQQ7EhmK^n|J(3FV|+~FnGcMykSwn=NFV>;?~zMUm8_nb>Ai4le#<7iHLJM8D^Y zP!(kXOE6w=qCG3#^qZDb4F<*)_h)hWIQzI`0!lnJW#RD?jWUxaQ=~jZHtic9dCjV5 z%DNCux#R=3PLeLW+P*t9N29KEH#A0we9&jqAy6B6BLqF8#s~q;T=nAiM~JQ*V%vY{@Kk+d z+eOt89a>IpPW}Yzr`zhSoQh+*RqIR)MWx{H>*^s=9BcP!-&r;8HeN^!D zf?44U-|cf)M27vJBN?SHeC#QBx)wrt*5Kpo$HBgT`o1gU%oF)6&f45iDyt7;&c0fy z7=?{fnY1({M3}IdFEaQ+@56xYE%c9iRArLTCkI93V*9l>zfcnF)2BOK>s>5?uS?)>O!244R z*p1-(Q**N7bAcHVjga_}XCo#Vs`W_!zN)!E*M?5^GcX?As)4>z`MgDjIm_35AA`eU z;RuThJBPZgLNEE!kErI57nfZJDuVatluv;d}ELXn>HT*z*nI>7(c~`z5|OZ z)0Vm~@dEohHG}_ah3CCzitkMAzNH7238*~xg!FzFy(=79w){f5kE`EO>%N2xcJDVG zTZ&J-&Qz5*Z({v-2g}m&y#kvfFw@j`G-@wJR?Ti&G|wJZ&!qCP^KC2=)d?KOdMfv^ z+qIc~Vd!jl^6*)nh3U&52hZaR>6Qzq(A#p6sjQLrfj+VCG)ol>3(l4sm-Fv?UY?bVXK;sCsEXw3|Kh7u^Mv$$ zM`u~jZuEF)xo#d>HY?XN*?BnG(05?kx0lBLO}0G@RJGo-9W;A6p0Uw~@^z{Y-Wz9h zgx0sdC2we)iTjwnqS71Le;km1V7*&6IF`?0WV@<6Z2TlRVD;@A`k;gdQwHZ?G=8t< zH^l2OMvO!rI~7D(=ZUf?c|s5gw5)ZXf#Wy)m$WN$vHNU)!e@kSIBhqgZARqAMDDC*O4tg0 zJG^w6#5Pm!OntDj_+_P(YlkeMs zbLt__nX`k;*=;18;P{_C7*8M`x?9sTs-X#=Spau6FlHqHMt zIla&B&?NnEM8#8JFk$bmt8p1&)?eeH1Oo2*B5O5b+8X{TxaOk$yS{6UGqr44G|y6C zic%*SCHy+KNCX}1Vd*ZB-A%hXv8jhs@Ry5~uQz1l4~Q}j^FX2>;Z%3OJh2V)*r7X` zbp4L(jk3=Zpakw=Q zJ|8iOYxs*BHTav=!d@|XLL;#fd0ea(B>jHvvxMJtuc8loLACu5OX^sLXzt}N*uI)4 z?MMxgUv(wDy3IMffrM|5h|sDoHOwR`k|1(4k#r)pfX<&2}6; z{d2_^u6yumvBBPWglg`|Yax?B_u*=xYU5Og;PZ>$+>6u54jNX|FRaI28_uZ*{Q0mR z*%M6VtWK$Cn4w3C1#24D!d_nD_GK?LdeHlg=gWR4J(q`Z+nHa- zQ^vLjV0x4sm&>M?mpUgVOdf`|lW_;)*LJ3`o(JIE!&r#kToPzyAJ&xqk9FA(j(V709#U9I%}|5u#d z*Hn=%jYO0c{+ARdBW|C{TB#9HtFc$wh=vMmpNcbU)bAbmnhadGM*>Jz@Y_C>7Z?u! zPtnKcs6Fr`bk-x8)!iBQS|d6@V{e%IbQ1sTz{D&abvZ%DsoWUZPrQ7K*RW0V{VMu% zk%mK9CzS+zj&Fz3vUeNN*{I|5Z|A-y&*NDuRU-_c9R)GxBg(#0cVt2H?3_(1x|}Mz zoUo`WpRcLlhm)3tdGu=-&-JSqgVIbl(gQ(@8_;NubkFm#a4*9OlLmSVNuG}SjDz^G zuq(%}jTJ6!&N&`HF4tr)_b*RR?>)ZA<&W}fte_Z2qEN475&$quUPiD2zh^u=%wy$6? ziK;toq2niE!<=(ge#h_u*XmEgBFjF~%<}|*58vZ=>Z5zP)U`}EbaM9H=y0OUP|;U| zzaW;MCot<6mGZ0g$BYI4fhvE_-MFzZM-LT#+=;*?sc6g6C+F1d<DC7F{lV!Hd}E!%TyN)6jN_`4+Ib;%>3LNVR&s7hLW=(<_BXxX@aaa+ zp`cjtmY%Baa?{}Ap)XP2*vHp+bvyE`?{FX2<^IG%Q* zy4vb}-C)0utgMx5c^@V1CmTt;WA@y&9Ur|t&CHvX#tp$PvS@V?js3Y9iVU^BNf+bu z(@o_L^HjUt-gUzULw`EkPV*45z02nuzTYHkizAnZk0q<54SaDHMiZB*S=_BbxCii0 za1TPfc;Tm$g zmWN$PDxC*<6JN}Zp4i^SVXx-TOP(C%$$Xz0f@S(aF1K$O_zG`oI-71mESIqy8!~%G z8q@9{_2OQlF;#xJNVWTdDOvAN-mn&wM)~QZLKO^cvyb9XbcPIlBhEyF#f9m7+beYIbb#*{=+_F90;_2rPtkb5ER?Sbw=zZ5=e}1(E zo{#W%S>N8%ZTJvkT6W>tP{wR~D0n_%+m+}0;u5!cQ)Thm`AFVzo(TbVqgvsX*mEoh zJ$b63m70MvkNe=qCSEAIb3_1^08|@HLS_=1uIyAx{2J@A_Xu6np2wB;^||Fn}||8or!Xq_6| zZ+Y$NF<}gpKuV*E?csOfyRw|(un$GR8Fa2j_C#2hAoqfhm58@-kNBYua1y89!0k;! zSWYSGqzD4*=6Pg1Qh>>$5T9!vXHXz6@-=egU<6*E;N2TI5_?V)!`quN{!Shu00l5{ z!jNc@XBPM^dbD7woP_Oxm;6O{|G5A7nRq4UtCeh&B7q&B=8%FJ)tu=mLNAYBrZ$KP zxP)4CT-kOw0?izaWgw3S?Yx`zVVHe}#?e%DN@c4V?*qIQdjGC0Klt3J!FuS7;GZY! zgJ>f_u-I!XAsH*D&>OxS6`d32BHD;3LnE%PEl1??)k*?hD|EaR!o#ey9Sgid3xn4G z`BsYKsaI*@X4{9tYk|}9zZL^$MMR<`Q3T2R)TGU^BVtmI;p%cPfbYnTcH8fIx#JsX z*ko_lF65Eax)tFsI;9Kbfl)kyhh$QO$JGV^DU+3K6)K04 z^^f|Mt7n24TS+=nmePR#LezUYl^e30cn0XP*4J<$(Qev5SffHq%4veg~uPWPx#)X39NhNZ%!l(z>oNWREAQ%-gHW5`obH9#cQ zo?#{2Zg)_0BiGNJ1CX`*P|ZHW6c?4Iw)n+QrwRPfX%w(L?$Nbo-c%hM1 zh27;f!kc;D{jKIU4*2;@@_Za!?$mHQe&YF+q#HUThWW+Yx_+?<;E$0?O=rY3Vgz}0 z`Pzr^+9@+#;IOa29^bbTAZW<$e=Cnoj(1$L58<`;tWADJ>cAmd2f|D)6D^WZOc& z8m7T0qGKc$JW&+o9RhZxHM;OAHXm0z6rwS*eFKRKm3lD zTTa$K-Zl~+9FssBVHge205aG04vkg}Rtsjd52Y!tk&X5r#Z5M6#;6l0zKYPDd!6jQ zr$pgUY!IEoi(W8jPO$~l`)^(Vun}1`Py-3wu@xi^=lMa3$hgN~CY3i-^ZM)RJ;ZNr z!qcs*X;FD>2Dfzg5 zJEYFdvZuUOmU{L2?Uq?KAO|IBLIl0%PBu4_2>G1C2`hEn_Q&!_%{Q}Fcu=Pm@Gw94pKI>FIMEN z*ydq}GqrL{Lt9fZB%R-EYQGEIpq}c}YfG%(&R~~~-wx-F6^1ypXF5bYc7*4XhGI(t zdm2&+?rXjq(_8t?qG|7@q4RERGbMj0f0X>`9q?ano9NrYvHPO6(7X|w2Yx*oAhG4> z34uJzo0axGT}vH%H(Q;5Hd1JMK-Cp%mNO?p0-Q$yOJmcxU6vZ4xh-W^7mc2ZOeI z3H0Z9saAz@JzIvj%4u`6Y%gfjcYMG!NjOpak#NRUGU+uBMu6@@O)aaCF2HGW~ur0WU53%~*jUqW&`oN^P-HCMfJ6uBAwS9lmaFMhQuvVogfvvY>=09@UAxp>Qw?;*^H_g{h5UU|2gQ?W?)YY;4Y%UL3FNJ-^DH>c(2A*=e z5VF#C_%A)vD=<40-3upGtu!gy{>>J~wyDmlM8k&^OB8 z<{tiFHT3DW>j0{InhZQRVV^>M{3LQrBzZ zAIV%3_{1OI4&SIklRtef{=o-o4jpJ5nW%JGZ5b1h3@%g~9!}OuRQnH!_(^WKzR6ni zSiT2)>f7@@vT^nI)%xVX*tc(?ap_kplLK#xT%QvEY{z`yMroc9pPY0SJ>+s@x_g(= zCHJKrdNl|W#Fx zQ6lWHj`{9RX*K1sjtoCCK1x24tlFGreBTJL03af2QzZpow8mRv{(A``-Ve`Ukm25O zlu~oPaxKh7N4^0&VZ@0NGj9wvvVA&B3bKbY__kRm%2M?!U~co(hI8UWc5Eqw7Ms_c z19!syOOhCPYZz5&=0y|X%A#-^9IoC0hjtEdsCc;`RCRp`U^l?%9#r^iY<$!JDuq?W zvPo@XsqS&AVY*uaUTD>t^WsJX_fch3kl`x^1W$*`?|$S{zfi|~;UgPy`w$d<^fk+^ z?fbRg<$YbbT>Upg&VuawXkMujk)`EHaML9#lO2W98BQvjl(81of&g3YYVeTlv)|5!|Aa_jOs8 zezE!C!j!f#F*gi>5RcwN_AAiR2y1(J82rChVD@3xeZH`bura^-dbde1X6>}vP8h0Yf*}U#K z<&Vylw~d+E6io+;0=NmKGX0Q({>_fJFB%IJYNeL+@;oNGb_JH1I*^4pU51F}gkXBZ z$y~d@(1OiP+5bKkZ)~+c(%UcCGhCm1blZ_)>n63>Roisfjw;qz;>(*$!k}V{r^vKA zaG0nQa4jBYr*xa@7;U@d+0oQ=%fL{b(a^ER!W&6E$hJffY;Oni*-@Cw7|tPd!@|UQ|{sZ=So3>0-MxugEDlup@%(Sg_iRxMgL8*1U2Ue^TIqCt^pHV zVHVXD&t%1~m{SQx?9RU%Bzu!nU-!N*|E(g#0Y#rA;aCi^HVn(KvbrS*gjh6k|JOD) zw^f6MekVvcpfilAwKx|3-(?;)c_~SbD4g>4lj}qkpWHj#oKK0OZ-2jTtF(aF~Aj>)L&Lb zc#I}3BE#|D3Ig`IPr6S92VzKSCZXTh^AbMO6$UG8lm=Q&={in`vR2ag|FNv zp2Pkta~mHLJ*`kYht)e+WOOc~(@AR^023U73kPsy0nesP6+L*ft$emBJJcO|)Wr5ooy2@?=&-|ElALNH8gDB4ZA~>7rjq><65%?vjVu{2@K$3~;`%&uC+i7#K2bz2VGmU(^b#AN zD_VolmDqH(w4LPHbd{vPHP9Af;2)9&VV2z|RG>c1XIK?wdFt}tSEa-+GBrQmUZ1g7 zmWV&HYjFqQ9vfs)*R)aZX&oKKhTfSrzvGBNkw3Nf<>c^t1gSADZEpHKwT_HGqE9-~ z2F~bv;01%u-DH(Y^XC`uKGiF3OMR+A_CvPSa;U6_gcA%Ca2(-%^Ir=sVXN}J20y50 zzfW?SX`5+I;E#ekN5nK^v+>DvI`quYY@duA6HaF>rV{%euKDJ>`O|ydg6gm*X4#(t zd#K;O`6b()2T4vy3(Qlj2! z$tRBCc`p*HC(N3XB!OiyWjy<(T0joO5~)>T%*k3UpYa*C^qa=ha<(dq_8vX6ojT=z ztWzHJMI!-R;dMRixq2M_5uBntKXndv(W=)A}-qho#EWX|c^5akZBrpq+yJEVcq>6eWg@=#>0VhUwcr<}vQ0+sNVS0D7 z?^&S8_Gz-HbOptz+cQBBN<|3}vPuHupt#rXIBe?VHj(jWbo}=loq~0#L$LWBFCp>R zB@4^%m6C7Me}LS{j&9}1$QkBXQmfN~TyVl9agTgcOG?MPe~-!EEoB&s7uakrS^p?= z!tZrg(3F!R_5#vRlXD}DE1Wgr=5}$MuLKBzB_d1(>Lo#SdIlOIuR2*XOwL~#Oo+YJ zmW%_b0P#L`vR5b7q9K61&wUV^R=>hQ}M0( zL!x&i3xrloHVH2#Hm&=r)+oK>A_NOlAXF($|0zmST=RmO4KG51Wy4u7{IfX*2x-~Y z9%gX0FDL=fhezz$n@MfowEflwlKSapyu?!n6D}YQLtP_-KfCs+ln6v;16rx?NDoK9 z1+nxONS#gDff}p@Mcfns28;{%ZDNR=;vE|@f;OpA%mmD4ydWl-fMCU>R*s|)1-~$; zQ_@cbKW6dsM{`ZNRDMmOXMwaOgB(jc*;?7^@RMPI8xiMnvKNcFgG)it9Y4i;v~&|8 z<3BxY9JhsinV!GS=4vk`aV&JhbK$YG*ln;H0m&!-|0z%pqVy2kBUbiXJ3~^Gq^GmN zuZ-V%F+o|UvFYxVoKZ6bYEnf1PVZ$CJI@^nz%vk`orOQCTjEa|+3JaoCxKe}tC71{ z8lIjj_KHL!#YiUYc*#9!68OqGF>*l4lAc|xh?Y`O@Nj_?E=h?4{MsbgnV)Ndtw#J8vkh&UoLT$q zB=(_{X?_mkK)DXEAnsznE*LG-g`i}n&>p@^+5vUwGnFC<_MR438a=KikN?%r(JB z4b9Xh^ei|h5TbJjs@2$>QrEbr?&mm|)63?O;&?I~)-wlI`xv7FCD|TPDmp&O1@{vt zTB(gyr}qdz=}Fbj0TLg?4z4ezF9c@uN&82L?@$n)es#1)+qEowE-bY3Au#49?ib?k zfK}(F_68m^#YOOomdASRNG*X;W~bxyFZ;!+%`5B=XzD~OhjF;!Pgw-?2p2mJRG63T z5%`X|FvBxo2YGdtFQS%@W6Lhb8o~~JhlI7TxT3hW>f0+J1CtBG&2~Uq8y_n{uZ!=% zIJ4Fz%A!2lGXwK#hB?r*UMGQhp4s?{RcyHC_Gelq9!g zXC6tL5;|yV2F)1j(R+Wj}Yb2W~X;&C8GNYuZasE%5R? zfUU(}h^0<|6|~hzJiv(gMU?JpRAq?t`QZ45Ub`p(&}rxZ9B9!{rxQ{0H1R_%pg#=` z;qiiC#nY$V6ee;OxV_N+k*S|{j$6VGl({2}p?Vz=g(mBW7h;ci&-hst=YK)FRUz6L zSP)9C!?2c};;u zzp2>bTGS*_6Fp45cVUyt;U6zZspv{IDI{cxPZQR-PItCsbQCV7A$IO0^JlxLPNwN^ zboa+K_MlN0u52z8ZnQy~K*;iZHo9GP?QMnV+irO3u5NNTp7KU)soMGEkB74?$9?D2 z+xV7qnE_ceN)GXvEe4xi;xg%GIc1v)d;&!b$_CSEQI@;E-sFEIPy4l6k@Sbe9tC)6 ziWw$uB1WGYzg{Hmo(s@mmA7|o?2wtFN&(SEhPjXI%rCbzd4`UwNwOo(EDva#K&(Lwy0s*64r>-ihsCtpcs<>{o(Er!3yWU6=%B^XoOzN z-^vi##Jgcq8Yi`>(|63K|Mp)HxMr-Yx^nlzlXpZGfJ7pD5)Yd2x-ySLU7tD4OI;(I zV{1;Yif)wn`u+H0rsYbfSBs{nhTd~p*N+M#C(hm3AFrLEU}xdPNM_87*TFn|C=hI ztmow)5@`3eSj4GVjb@CSvO>9~0`lXHO`=$pnqava>kN1N6gTCOa?9D@ei>!G)RRdl z*#@mF$0zLHW}mOcUH78&8YTT2Juus%lQ56QuELn+*T!==eQ02FP2K#P@vhg!S&$!M zKjsuOVS?B4*QCuh-q7g{-!WUdS6mpQT0jdKSI!y>3Ea4&weh!EN?9-EBpymu`ce)y znlMjZD%{BU*4D%SxZ8M3JKGNFdkSk$}OIjF{l$nG$j2(|5Nza3KzOFWt6K>CKcrroy0!fP^j{p|6D*1Mj7#9)?Z~u<%ZzXahe)h5g6cJuLiqd_Vj!C% zDM5>9Wam8F1FD@f6a&kstT?68#SyrgkUg2EU9lj>Qme9h`>jNTV>m!}nNMl6@WTdF zLW()AfLvu(yT5OjN{713{PDw%-sVuDg6y2-J!>Cog~a*X^~1(Pt-Qn>&M{UgNm$|m z`6Y+O|CVP?EeCEUNwT@}`z+yNWdm^LiWf(r#)w++dmWfc9rOUy53Bqr?j?aOonp#7 z5N*ESguyn*8No%5 z2miV75fO|8L#PnGg}lkJv*gQvG#Bbrr@i29K9BED{gpSrecB@^D$;tB@8Z1^f}6Hw zsGcAftWja`lD}Ff<6j^5^C`4v*87n1UWabG4|n)bu~`+*n79A}?i+M2fyfmhq8*hg4XEG2VCVrKjF-v*7%Xnmh% z+pJW);Wwv}J@hGzClk3iy0+g1GXl25^@S>#sU!_9K{2O{QcczjR-3*&zbRyj-8}m>1}-l#B#4 z^4;zMLo0vf@NtG&i=Z=UjP9i9x=S^i2Yw{Ns>&XgBh$wM&|nlQLMavFO(Ou6GVK3h z_!FR~y=SFdv=LA(-o;zH%*ecY83MM)5xt9oMXobo%dIk`o-d*AkF`i-=v$Z8sr%?T}F%+Q4c^7~3L;>yWpgo)pQ;`L)x`Ejc|0Pkx%=%{@kOq>8}p&x#B1HScQm+~5J3cEW4@Zk_(6f3EZ0-TWoHwyPG$@8yrc=#=kCTgt|67DtMl1s6S480ou9YWaw zKtro`W&SA2RarCR+=1A*!Y?Gj-0V`-=fnXoye7XL4B@I+VRU$eBo&r?Rn;g>x(TOo zLd(yPF1gKNE1f9W-XHy+zf=!9Bc#1Q0N@BVIW75^NBTQDZ63^npy1Tt1m=BCmUsc} z0Fp04w}0Kz=`BrkZdWhoid$lg>R_~{a=fIySpU>^Kt<`r#!S542?d{l*C-Ts`l9HM zMglI|7yN|zy><=>Z$|MGJ?Q2TXYB{TgzbQd(3c#9;=%QDp4XB3+- zp5Z+L3MF)umkyN?;VdVG0ALLDYx!Rhkj5)Je~TP-H==t?&^(!t7uN6$-GM<16o@WR zr!Gd+Eya&-1>nn#lXaGk?`47jfv>5$#Lh4|778tE#A7L^pC0Ff~sMkD|4F|0m)t7jI<`!BGY3Id(;CJaz`YItW3C+# zVFHsWgHn3fOQDe0CP|+%BU#%u3h_$jhRnh>DgfdWWRPp;I%6im!C76AO+KPAUQWI|qi^zJQ5F9WLF zmYYj4B1Y6VBA*s=L<+9dMqZ{ttXUGdN@t$v(?jum>v_R{vP6cPh>*n zh%=(CCh>O1#4;5KJ>gxXT=+lrmslQBrt&A5;k;PBa{G;p%YXyZ4@w+7NZNze0gp~P zXQq)%>9uh^iGr8470?7NGv(ZiKbZH_OFq@xA#Y+`+e}hgV z){H;LNeh}ls$*rL^!oq=DV=NSD5yx|=e%GZJ@39wPA>7yy9U<;X; zgj;xBYf1}PI>QA9PeWgzxZ()Z19qRMWP(@J_%rRJ_Go`qwmVH?G*LN$B%Djo6I@WI z0~jpS%LauWI6X;zX?OjfWO2aOsGYUf;$_R5Y3V`m3=tya`4txc04yFM%rw~UuKd&J zt`HY_{UZp13$}J|R57lU@E&~<=P&GJXDL2sKYYb0&s7lfZ(MUsn{N{Zpl|c@Sx=9bC%miN1g94kmIYg4+oeO*9B=`vyv%|yd>gJ>jcKg zr6p9yfgD|FbGCyjQGDJ5aUe~S``HI@Qp|uQ?fmn*I%r={SklAspAePGbcyg5^+d;x z+^W!;=NMY}UrvId*vn_MHIQIw7|&)`B6x3A#Q{tbwxf}?Vnzp!I$akcZUS+8rKT%% zr`PSt)s+UvwltCy@e%$3LnJqIhs^{Ue0J1h2)^d*Z8c|pzVo}z%W)FA9yQF)^r$dw zhLeQT?GKhe{Zs267kS0OU6KU#z)<+a7;YUU?4WfCtwf+(tzL$ z9UFADtN4<(f0iOS;J?&k8t_A=>m)Bs5!xb}b|+nUAy9?7OWU~8@k#LaAaba^@FT8v zdO#SV`piJ9;{r*+Z%3B~-<{k0WOTQn)~j|^j<#xvZA`rLaCKP}qs6gtxi;a*Z zO1;;&{JK2o&V*&eHt$o+S(|2B;E9*FYT@Z-@YFuRxQ|?!U{P5A6^{lmB9D18>RcLLvCu|YN0O+XT)EAa zBr=*6MuG1F#z&pVqvomKwD(8xq zZt{po@Q?x1(qEL-lPVo>Jlv08q7{3kvg~j}t8WV}tj6{4wB=1Ao7-p0qE(l7twh&( zz6AYY#8VOr{;<=)^fm1AV|$mrCf2()4TIJFj4(%rUxk?+L(3d56cm~!)l{BSTRmEL zK0^X;tl*i=Ob1~0nd782D}9TfpOZYlVkK&)%_A22a(F~`hvj3_mD}`5vlK~R?Ajg7 z5ays$@w1^z6xa78u=WXOj%7d?vp=V7pr}lJq`ti3Y>MzjH{|P zqigo$WrJz&aV53o5N~$9OEr^S2$9eSpd|liW~xK6Sw7!=(gB6(^)!xwhW;)QpIVNnl3I) zGWde9-y>c9-=vLgm*wNmm8~S0y2w@aHpw_dB~`)(LDOyLQ>11-m}zuo+Mp2npAz=d zLkA*79j7V>s^vcR>Ic>wruPjtVB-hE^yy^P*alJ)=gyX=BC*k7~z-A9QVmmyZ`AkUQUga$la{hz!)o65z5 zJ;P?@&ZhsI4CexzbF(_F8{C}kBj9)Y!3IEoNviYutJNsc-KY2%Vi&*&JNHfo8;D#+ ztBCmwyxk`0-sBJ?`KGTfIdXQYA*Cg&V%43il^J6_e?h69-1)xVUr>?=p)>O445eu# zYh4W|%8VGM6i!tT=->c_@2cV+B&&kS-&24YhE;l2E%y}B^{M97uj?y;T!g2D->4DX zNxDlkDfonJ_?aLa@Et$>{l`gCvD_E`*_A!57u$rg#`AF;yp|5o?C$zB<@76VBhqTp z3Iub5PU(U_MV>Kg$FFHevs7SoG~0Q5G!*Z=HbTqEmGJ63fJ-aky{kIU2>ZE))Jv!e zbWlbxk^cr0D&vXI-@%fTj#^-*xZ=y!q zv7PThsNTO%pIWhDSXA8UlP6|YI!ZTbF{PU^_&|MkHVF07^YYt|q!}CO^}kRh9T`i2 z=)-G6g#8?wiX;T&fQo1#wQP~Gy#!u1$$QuRg=@JD{ZMp;L+{EQkpG%1yVUzlGg;~p z+YxQV$5^VP28jT{LEVYohfHR7v~u zVH-wtC)2~gg0HAUZXyE%+PCr13fK?T|Me>tF-)7ojTWxe?_ZS6H&UYcy?RKiH>V}# z>>vNMeb9aANXv354V3ab;+q!Kc}Czym*ssK$7^MceJrrsqluPJLBNzp$nQV+lkdkN zfF?ZMZ)(?^Il2DCi6p=CkW0=!zwXe2sG`){8K8yq0uZmo_%pRLz*0Wo1E{5*taMn^ za-mL?c;h7&D7uxJF?lI!|n`o)RQc zu}q5)LS6HtRDZ#A)BQ$olKj*|i6##5pi*EAOE*D6fa{aY)IS0Q;4q|t2kYrPk-$f* zesviDj6!khG1)h#?>z_r%hrE#IX!{hKoklYd!BmHW3OFZG=2F;2trK3{g`z!z^s>q z_Vmc7#&FS-A#oIoY=snGRVhr19dqW7PQ_6X17IQ+i~<8W=a(zIINPI&xOk5)JHe(7BTcWf73 zBLO%h`afOv0=3{NF~i|ZUqmgM`hMI>{+wur`cmXYu<^Lwely~bcgFAJ+62XhRbA#E zCbv006M)#I9-~eU|EUhK@{Vg@l0%%_O7<^GtA4bfbdIF4cOX!F%;vBek}SGJ!=FWr zIJun+aMS44?u4eEU)>)m?ielFQOVO0uV#9vp#~#v>hs6f`3O&~7{w|_nw;8-p)sMC zPC?kcaOu2#e13qCrIdZ8?nSn8_QQtIQs*TG*+bo^S$~gfz5pBtaO&(<9OCq<d9Q5*hI%j+L&PH4CUOa3>?+&h73_LRl*FVYjT#qav98v){#I2<5}r2&2BS z^FOvVPd2QJ7$GtNxN}-JP3ck}uqW;AkL0Hrg#OOH-Y6E-+~<#G1-NWXL+=PhEk3!v z?u9f(vv8WHBzP7yC1^^&%^6E%-VPJ^hvo!^AjP^LGGa@K2VCmet~|jk+!y|6|7;eK z{?WBh0sjc5IC%^rBb9X!}8i>kr)NC95oUbvixrPEB17?dQ4~cKu zQ&cex!;FMpr*t=bSP4S7Ow_$_q)U*I3*hh{I(UtOD0U(EOqiZZ0X(6l&;#Oz56?pn z{w_0u2^eL85aNV?j>!HKjBPyZU*mfT@Olyyv_dt|cuw#w3;L+M`0-L-Rm3Y*NYd^z ziI80sgDgxZs^D>Q%zOH;hCgwpj_W`xjZdG>{9~Z+-@b>=p>e)Vx>43|$ow;82vXLf z2Lu4()mx+c&;F4PK#fvz@NA|!xkg#4&Nz>UAQCDG6k?yD*y!${;QD`Z(6S<2Hvf4@ zd>|l8khwPbs#-5$_pOm4&;QitFWfmaRf6}X&Xkm{rCtz9#~NZKRYgAqGvfW0wELfg z&b{h3jb!?y2W39~7htRSv2;EUR7GemY(oh^GclQF#N5ZBy zbNg}QK7*snZ{Gf>#E2?!u?!C9_0O`euK`UiMQLt25Nz1{lNvu8`@|d!@to7ipM9Ki z6JL}-1CcGqAwi)Y);fuPeKE)69%xEi!3*cmai@Y8dp_saSM)E4+IXSai&B%1L=`*f z{kG%{2~`J$8qCA9ci-$vxR#uB>AnunPT8fml?rzcjZ?3LJL%-)-S+y4lCPtmy1bGd zaYFp|1FzB+sW?4)^y>^YV6-}UmH3mwlijQLPOp9*Xp7c)9i9k~s3_U??In;v=_28v zPx?mt{((D5q=&9JTKr{<5)jP_)FPOnfEPx4d;CQ5A$3;qcu&OkS1HyK2^Sr9&ffSS zuODd05R|VwSydqVK^c#dp>Z4_pgryXN7t9YL-oD?8*2#J!&pj^tc7f2ZBf}OWr++4 z$ueUf84NxmyCh^Uk+M|^na0>9Arx81491eoSTgqi8GU}=*Z22-y?$OVckVss-m^T< z`+4rU=XuU~0ER%cMSzl~A9C6n5Fb@n2h$QUChsdML6A!IDe{{&=mHAusW zNZ}&qS^be5n02Y~gPGRWU3PS(5UR=mW1)gFY4aR$Vg^bP)ofKn(DNJZf4a5q1$a|0 z%h|$cdb_MEueKOHw}i}FgX>)6f1(Z{(lzGnrU4~yjQ^!%9&?%)2+2j%jn!b9aXAdQ zuIhSmA^~Dy2s_xJtfK#BD}>>8fJhruLX(Fm7A&MZaD<1L)~f#=gjZ7C;L3abxbpff zEf##pAp|STdWi?fQJX)qEvf`GS5p?=h-B(8BgpyJLg?JnZUG>Wl8Ao7bA=$*W#LI$ z9e>+ojZwm~QFE!bpyyX&V6hfLdjp$M+1*nbnbyuhOdWCz<54J8M+Xr25C~_7$VBT` z20D_u52*$HlM_dn>A>b>TXTbo(|PtQ)>|Y^6b#EQEO-}$Z$n<1PnaEB7)q$Y2+!a7 zn{T^qdFxh0fM<6NV=JfP9(`YZmeGTm!ht?acuUZaF3$Qgsd@+q$Mh>Y??2gtp@IKq zFX&WC)q!lt!tgNL2@2CNN9U%+l>{8?X5-bO6*|Z#qSoBI=?$lOQJP}demsaoS zM@uAOq`kxJ$~rinK=BX(-J<_A_|ZXQVTct$?xt5d)}ER=PC2YKdU#n12O4&m-UPsI zcga20F+E$cAf_D>&#VV?=XdLS`Y2A#&H#xP**0=`W2rQbabYSa8&7iXs-MHKQ56ZU zrW?W_z=2$>#THGhtKpY&@2H&b$w&Yu=v+8ikpNs84k40qKn=!dmhXv2G$PvP=i5^D zAH_4c!B8B;e!X3;?r=R@YZTK?8mRXiT>45h&S9tVYVL}Z9t-s$b%h68hS%071A#Rk^lfAhuqE0 z1^_{GqS*brFdH%904`f7wvCUY489x*f2o6jC`WJ~JdD6}BC5`HS^StBV|de}zuii2 z1&&@_+3&hH+eC$NnqC;ok0mY8lihwUDITI7-lV`eTF6KKFE1;R|%nVkej)`~WPb}!mc`(n{`M;|^lLX$eZCVo?He%!)I=*vW zi0}|C59R-%%3pMt#lGGp99Q_azzl=-gC~ft3>%>JF~XN|f|}SXzHL4Ai<8ljc0l!^ z>5K@2kd}=;%k}Pc0}=S!YnlHA8@q@})_vy5^qykM%iiJfZcAArStFHLF$k9KRt97n z4;A!lA8epsl?6bn+Cc)nH)t3URMdyax4UU?FBGXCRpeda*C}U{3QST(AL;FkqS1Nw zn9gl?K*M{2-98VO7H#+8JHhC{n9JAlFD0HB<$IJVM8~a5;piBW>6QDpsIX zUdw;1ezcs|B+fld`-ecXdaJV;)y5L`zfQ>|VTy8pJTAR~o5bzhy!LDG(d~vO0$cBk zhAL%RK~~^S$(v%$9|G-9i8hu{k2q&DpE9{5^;VELQ{rbvwYP)~3iPPA{_oSA{LQDg z;g-KV;zG^dRLcEPZ^ac44GL^My`7&w{Peb+rRUbxcV%4t-s*>jgt1*?{t;3x{s9<|&@j$xtP2{o4J;(}ffgX~|dBp38L>o;#T)>w{){VJ*5 z8M5%%dXrE5BUn9M$g*KxaBEA>yZ(>mXy78J*C4LF2sT@A(%e?)-NGi`$i15rO-x(Rqf2;IuBBK#R9zLvKmF&D~;4U9@jpR*6E(_pb9C z9a&akqNlGF&;hlc`~-xLjjb%$$?y;x`jebvH==JB?zJRIH7H?4Z+U)tU(ovp&bf!p~|R1W8SSQE>ptxgz0m<(^ePXWkW^C+}>V$@nQ5&nM7;|Nvg@54LRb`3&& z`hT`r<(PlUvrUAx{(SO;lOAxr6!PyyPI?1=`*B&d#$kC1;tAYtCOl#-e?aiDpl5^R zUSqN@yfNZ_4x&MWe0+lw0pjeP+>pcV$egnalV+lavzwoejzomBU(J|k!q^pByuL<1 z8T?ONy@VH6sEf&AVU)$>2nV}idSqW6jl*j*gjDf<7+}o%_#iPKxB9O|5lo!|r;RTN=3ygAj zg_B2QMAX3>Kf;`g7HAf?$go7ugnBd30;!3Mv%mimSs5}>ar#w^@JSrT7Yrz+2_r66 z(?DN?NH!o04WmogSP zD@^BH7bf_5Y%qJR2m|OL*NphHkMvlM;46}g#&7^@c+1e%@TU5e^MhPRrsS)}zkqm% z-XoFVNifgGRqJJ^vWG8kKTYS!Tj604_^|v6MaXAAT=Bn@5$n<)Qcl78Qb3i*k1LNK z1)5nD`5L`DEmH~ICk}n+|8y6;sKAr6`!zy5LYV}9uC50vNZC2NpqbG%F5l=+hJq(Y z-p+35Y&UwzOH-VC5R`4m%w!vgA(>6@GZVd0G;fW{oDAJ;JRPNT(XLcpjeyu{F%0OejBX7|$f{M0q)HWR!Tn;s=o{mMpD-3wMV=UY6~U*Y*-mfss=@SYbVQr0H^>Ni{G}0fQwRf zr_9jIqX!x2HJo2IYV^Q|X3T)wmA)#|vFJ*a*Bs@gQx*4ppQ=dhWmXrK4_kjY3p&bP zZ`v5Vt>~qBhg0F`nk8C%O>B1etp%RH{YUKx)1z26b9oS@3H9Q0pmkvAr%u4h>`iQV zo9qf}FfDh5vy_N2;6h{#|BQ_!*2JE58O#@OhP1LnLCeT(i;r77&m1%#arq0OOTft zxVJ5}Wa-A@IiABqO0N23`wbADT+Ira0kmgLmR(W(p?Q3S7sIGF8V1)+$1Q~`1_g^l5ceoQmx zdq)|N^n{CjHAao)y~lE6Gcf%a4NW1+ZV^I+C(x@qhDtv4<@MV}FSHPC2biXcaOT?kR>TS zvJ1O4-bu;rj7tpF*DD%;7}PcL|gt@I}JDNqoBFR7}@!txpQ<_R3EYlaT|qZ z$?D9^2-IrbC_qEwIa44X1iFI z-^kcSq^fV;i1&j}ba~0RRr2=T9a3v;l7NrRw}4JQjEr1_}AiB z>FP)N3buF7hq@0QaiLAfAr#YND;IYs_qMAy(_citxM6$KD#i9@dPdfZxi@?oOK_C5 zavumg9|xi-8>RUlrFRz^vJ3em_l-$cl8D>o1oKmRt6&*AsEp2HfTe#~9Q@6PX%}Nm zKsL{O)*Mwqsvf@-(}reX9%qnQnCq%~a>K&u!TYkdeg5bZgfz?j)b%<>CW6Tj$z+QK z>$S&I)nVUjML-0xwqvc?fVBjc@%&W^srtcn9xZc^i;2&Tbh|Xmk9|*<(K>POR`=EZ=6VQvc=&$E0i z7Xzjx6$;=p1W=!u?Z&gXSH>)VU+THB;O+Ggbp#Pkei$qyzOJH2vc^BS*>gP1NA>ij zm_-+Fn5hg)2GbJz=Khyt$pBjbM*7H^=cP*(fmVQ{d^76jgBbrrmbhgta z?ISU}5};q!;XvaUN1p#5cequDw?L1nMA>=SIZbO=SxD5N>?e!U0f7|o z<}f}r<9p5`Jq%10Yab!|WTRx1`;^-R^Izr*goANLF`nbCRxKOFw7*jln&*D==g;|t zIxFIrTdy2KH;C!jhr-K8bISK$>Q?)u?8gdlU_f^wT1$pymyWIS24IX%!~LfoOk3wv94m8KxX3QL{sg zf?aB_BTilnWzu~1nQtQf9cXnEPy@HLy9M(zE%^02T!uJ?s0-9lWxb9x{~UPRPPV zuj{6m&@bX|hErQe9~xGo&^ND3sxVnQerj;#*Qcofa5nb|&WgW%sLXS7L*NG!7EXVA zYFeY|ZY&!DbhU?(Lq@ralbP>(coDxm*quZ$#H-2D8nb=an~#1Kr2Q#;B|N?h;+O5@ z%FbB&+mqw=yRKQ&YZ4Idqp4$6P|6ICzP{F>L`}SSBl#B}JmLD9-gwKp^bOsICTO#3 z-$#P^C*rQJ5f+8@JDg1Mk}52EgnWw}*eJvJ-P4q)W>qs*t8u}J1bN)mh8+<0@aRq~ zh(>s+OIA+rFFcrLWYyQvt(MCSEgcgK zGxvn>gJw5$-)=P`q};6W>`4wK+pXEO11LAv^PAm*6F?)1t??6?xL<#fz}vsygWJbZ z^!iJiX$e3tN2*GEiq)cW+1%)v52JQ79j6PI!dFQtC-rrGx~ec57vw1!SluJA1#{LO z48ikvY$QuG?hu=6{XykV?wRAQz02(xpURhE<#D2P1_1>GIHY5}ZK&id0ho#pw%%h1 zE)Fu-Bfhj4DrXU`|;Af#dZ6DX3H(?;8OL$0VG2+g7PSp?$$X$QHF?A{i=>C;USy-Cf~ z0N$oxJCo6)x%cd;K+G-yG2>}pkie6SN^&e(#b8#HbYDam8(pWdU@tbPralAl$OI-{ zrX7+*+T5}#NXlAz_4?YU zYr-2r`$EZal~b4K;`({Qpm;F84caDih$7aJ$_C@j37P-a3&>p}C~u8+bKNJ~XMKaF z`REOY##dw*9B}q9P`MHcYM`!zNIepwpwk-v`R(4RdkD3@S`NZ7E)~p z*ne**%tubhS5wd@C+$K>LUU;p)SP7-Vo6}_7c&ohtvxIWCBG^+$W@v+V7vC1D2$+M zHQpY&0=eu^v}8|>hCVnCJ7sbB!Y(UtI(|TZmde_H)Vx?;dn}V@S%`9Ur@JeTC_PTs ze9{u_V@x^6s$1j%GGZ)w+b?zccWUNncT7HB5ebiSXCV^=#j92{@LGn=U$ju?o9ihR zG*^3So0ukec{$rt)1lEqq&V4dCY#Wz-e%xM|GJkA7LgNYU0ukTFX$WR@a;hxRyFmQ zI5{qc#G7e+a@N{fsrB->;C@S!W(6L{CussHsptf%Pv0@@8tL=;6RWSB{1bLJ)w#B3 z;xa~_j_nW=u@@(_h{2=(ph$5s>m)A6%}Y!-Y6O5=(#ciUS4c3|?yqtMoGfS8}Z` zHL0nA2uK@vTKq5p-8Ny6UWu{bkK3$?t~A)O^M=LDPgsrfT7c1V8<>ed3dYTmG>r9M zFuAz1uyG-|A?%WB1V`9G6(uI`f=)@wIm)T%`#JFKsnkdQ9l|N!8>Zv@US61q4{Ftz z_ixo`1Lm-~>!di1BqP=2cq4#q#;9oR3{3Y8359_Tpa#oQEM!i9`~arcFHvQAG6f`k zXUiDIp!rsJ{QWqWVoS`8VOj@YPO0GQh^on)OaMFD^0Nv8B?lGl>LN~}@uH_Zh|X~d zzdfhp^n+)MXlYb-Dgl8>i%{~M=2IExRd!mvbE8iopp(#K3ha;2PhG>(C28NI^3`1S^rF)Tov_?fO~M8OwE;`81M z;sTd})c;Hed}Zhkb}8(RLJ2%AO*YDv-!0(g2mt!gG3?|CFBWvm3yq|3Y6U9CGvJ0i zFqrffh8e=D14Y_+>>qL--CD^U)M1@ke=y{G8( zDngO3!ZcT{k$pgh{sPg*Cr1q+37Lf2Cu4i7PesAID`cQ8Ur`_)EoBdR4)CmM$y;xK zJ}cX*GtnN7`RChO8I%zW+G=n1t|>sH&&(HGAsfQzUA4ePb{hG-zfMPoD@Z0O5>BiZ z=y9m$MWO^ZAuIe(bJ+82^cq{_of~S2KJN?aZBeN}Z(Hgg0VYVmJi)?x+dAfi=U42{ zXGMz|7Xsy|4Y+vMRqw=OG|R(i;>pVj1Hgm^__`y66x^g+;HoclJ*7>pt}UCe&3aoY zLF~_iABh&QR!VsGFtEFWSIN~A8^AY(r4qAEJS1gY_`Q~xSc7@q8+=L-8|*VJ7i(N7 zvMGO9VW5tgo1I5!9u{j|sBtBl$41)K3{`ALIocTA&{*pHES}!Bpca$Ex4N4C1DX8?Km5PeC$Q(B6$ywG{Q-H2%TRzlR zmNjoFf-mi~lpF%)t*n4i6EFejnAyYD=;wMI0l4Hn26kbgc|hzF)uPr|DqqWUDf@tE z7Oj0YsB-#pLwp7Pd9CzL3K06_DnymS>J9+K<$w-mnZJyyWk!B4mERy5qS=X@e^vh$ zO#|Jkbaz!ts_Qee+JAz)47px|DnKbN|dh~gUPq+5ZbBDp1m=`}Xs zxB3lacry2LdJb*QS5qfCf)6D}BG&?irT~=;vo=q}qFEB%y6c#cKhFr=!?FPS{M(DVAcs$0 zEO%wBaDoQ%--6ekrP}~s5=co{B#sI!pd6e{)AlwgR6yf`3p+v-M)73-`(U}B@Ly{x*l@b76IyS5vtE6u2qVn4S zDXv~A4vW+D$6F)c($l`Qc7`09h(gKi{sn1~3^Ct=cKlDqrwU9y{4C!t-I8w9odW z+wtGP6Jd)ccB*YaC*_c=I}4br{(Hou*cCu zu?dcRQ|w<6&Dgj2i%n!TcwZU{3!c`E1&s$z8k);JUHltql?ArMsZjY*1s~sI6*@wH zv>pEKkqfbH%j)oRcMwl7q9RwUNV#e|ZG%0{dY8aF!V*q3zBfsC4YrNu z;RVaX={^QUu}cAZ2I zjnT3RUU#j57Xpx)>lDF*$Gt(CJdA!HMPOFiXY4CHyIC7wRf{|Z$Ou43H!12-SD3Cl z)wF!hs-;Tdn5lF>3iNQBhJ#4AVC7Y2otG7WGBUBsboYmEB6BIW_odSsU~D3>?zgt%XzRFB}Q6jkEbkW4Y_l^ z8WOcGs<3s>NcnK1zoUkMLG;*fI&pJSYF9!dp#&jZqJAVyt6OV=9L^lJd6T8VR#=L~<@8uI9KstWu(w{Rp1idySI|0$bv;X znCyO>bcee>VMgFsW+e0e1&eVC;my*4>Xl*Pi9qI!1AO|C%*mhL%?*>YDnj$-hQ-Pj z`*5@~J^qVa=GKqNeK$(pEPd!7vf%Ig`_I-l)Z=_jW#Y0kOCz%a3xBykD?2NHm4In= zUtNBCQNXkG$cI+*5$_t$$@`$2y42QoOwauN7l6Irc^*0NN0qw&9|;=?JDztl#Vvl> z-VkudFcKO=p&a|Zb2jt(I*-$}G5voo1+cj6JwB1XWBKg*pX8cr?0gFQ{VtP2Au}V{ zV&uS&CbuFic0y0nbb+bqaU`>4iX#W#R0dq*ggm9BQtqoo`*O9gn6+82b+6SGs}!rN zu#baj3#4IvNPF=$`6bAN-bek4o#$I1<(YnkR1fcz)*pB()OqhlhZ9^+WFqXmZ^DlC zk@lf*!N{E)9O_WR&PnnCoJSMROnxG4(4ToiU7XCZU4MLeG2_nAKc8{wI=||VlSUBK zB;4L6_5IIBOk1Zv)BY^9KWAYd&bX9_-r6~F=F&|X$1f9U(ui^(Pni1@dh6>-V`WNC z39#$FVg6hL&yxHj&p{elT747st+OB4Q&kN~Ba+h4#!8ta!XV#J@ZD+m9_thZ!{x;u z828qVf5P;AhboVP(F#b$ziF-u>+xgS8|!Y&XB(FajhN0K&EKt)_YA+VO03W4o#uSr z+GKL5VJo*k+r({LcO{@4AJ*rFUS}#t|Fqal$dJlwkg|+y<`bxaQ?z)9K@FHMLx*D; zs=f^F3ZplM_pT|!x8U>BU&2bs_S8hW%0;j9e@6dQUd{obu+(PefWSf~b{?z3|1e=K zH->m9IS7*cia(I-u&L}imak1F_#pk|uS<$92L+C!U>ooVA8=QJFc8W$($+iC^W=Kh}6PIA843`BR8EiE9 zS|m@mMMP&vh1DE9>Svh_Z1;i*@$8}#IU02b zyby<1^grBF^s`8fUb1Taq(zh&%`57b!)>YJou&67#MoPV5eO)CbOR8e43*2kf=lYv;6&1Ra!O%eSW=Bc{ zvmEA|=h_rTaZ?!-cH5>C>i3%Z^Y-4HAIjSDR6y5s==kX;c7FNe7};jo5~Mw<$Mkws zwWi00<}DaHb4J7_agQiNsGr{3z((*vO?$zv9H|8r*StgCQ1~cg4>&hIeTt4q12kGZ zkXn~WM+KtZb}C>cZhIU{TwyoU|IpDfntOXvF?M2+{6&|rrFIf!xIRJdL}&c*u$1v6&uS2t1| z(`It0r+mf#?#O7$cJktzjC!xzccJW}h`7PlZ(hKjTW`khH-WgES;y7X$yTimC5HA>jfDX%SPtIE||tqSRBF^saGy0qIK`1tGS6~{?kqPCQIV6juTVE9$p z#f|))VqBIE6U?!ukrHo#5ORn_{Z1%5OHx8W;T&*SP3=!*bWN1sj4J~t`N_c4=p#|k ziSxO2v*DAYCLtWiZ+lI3Zf62g3aHiOa8G_p*ZFPELDm@!l)_VnkXH-PGb1wjqwD59 zUjv&pl0NkZEI(bmkAP<5tmmZ+Ntdt@mH5Ed?VsO`%_n|2ng<(69|sd*rBTd{ZajCF zL^{5nDbyI3*jdHMQXr_1=0@qsMXsscGY>{zUAw6{-iGnJk@{%mWbWl$KFV2bO_u-` zlo?k;DgN3S^oIufj>^9F>00pjED`@}J@^QJTlgE$ip6C(U$Pr}sJI-V6j=!phpl z(34Bf(9q0;6BG=9Fu*mF=^>nSk>;<-}l7M4wh^6wIP610!u=v-vSQ z{3uc^9A6#lM~U3hAP=+DP@{#B^U%JCfB=@-Y`@dDIQai@11CxL9Q*>015Y8Oh2HEvnIJ{shM`Ce8Geh!v%X3rS#4T<>{#aln@*!zKsnBDljV_&V4n4MC8 zFJ|=@+h>tmY3^4;@%%K5JFnZNn}XZKs+VCH#`tSDQ5g&hCoet}ngcbLe~xC)OV9@} zh4A6Wa{7hKcg^a?9Ph7SviM4rMQVJytybM=GjSGgW8XbWsNZhN_J8UFrZ9!6Y19t> zq#{A6=^>#x`TkLaiYlqdC~6?yjpH#!vM<6CP^2)rUe51F3NnQJ%OJE>*f>%a7-Y@c zqz$#6<;hR4u&NrK*}1oi12dS>z)w95tXV=aX-1CdglzbqHY8o76mj#xM3R}6+j^MW zK`7)mcSZ+=^?D}DmYcIt-%^nwZt;1@jK2r@yiE>dM*j2V80p3%PNUO5qhv>?B_rmz zbe8@Z*=p_d34zlN#YoqFdiTw|_6z;+6CBcmn=qTKFR?0yq+5)_8)=SN1fdhZGT(nV zQ@XZAQDd=76hG{V_bO3^tK87+ZP2&<*M^%ss39ay$PC4r@12RF8Iwe|`S_ zrA1ndMypsQvC1fDqsiQQf?PyDLO<25^344FoJua>UdDE~byN;0Y%>h{ zsIL5)>fmQ}RB(J3GhIE6Z8S|W8!Kct)N*HFjY(s@YjkC3z)8!?2~qMQsI)4k+4rXI z(;doDy_Zoi<>f4LeT|oxVJ{@jhB{q%hbJoT`6_79Wl%WE8*o zQE5ht)66E^uD9?7=h;?i75;NEe~9N{sXff!9=;H}4oTp>+sG$0E>!KPT+wAS=@=QM zeUWlBoUcp+6?)rURC>O#PwnbUHk1CtS9z>nJ~Yy`DQt_qU)98)^~NYse9bZX>6rol zs4U0ZhN36$-8S6hC^#-uRQL=WDjC4Fs1Z1DvA->7gI&GioR>kB)te|V#~Z0X8tD4> zkoz(#`kp~4IA~v0pW92KQnpYIX&OP77yaiT){VtmlkKlW7AvI1<#>w6LeEu@gtZCm?mM(I# zv2QdkY^t^71zGH_O6jZ)-;b#uX)Nw4rLvrTj@K*z0S?c9j#SCaQMfB*D5N&dknU(I zqLC>ZQpI=hcvSk`qu1Ev8>Lij*jdvWFL(2#s9X>Z6tH<~_fzZO!P6Fmu7q5sC%>!$ zxtX4G>ck8_8td$Ka#hdU|P2BX$7!|dT)oidV^)no2n=d;*r`&8F89L&~ zg>aU|sACgz8GcaC1TsE9s%rJD-D=TMM%3cwCT38!_!Xr14a;HCi=uXzLh_A7MFy*t zGandizK#qrI;tZ`SM=5J)$8DR!yP*R=4p-3OG!~T(U|*XGhWd7n6|98BjT^oVlk_! zsZth(VSoWX#p!?nnooGRR1%#xdNunspE!R-kZknSg505B;2vROY}7)ZVxO+In*=Mm z8q7>Qk56>|(~Ixdbam!Ml7?Vv?(WcM?ZyTo&T9x#X6U?0QroGm0b0N2Ip-Tl(xJo~ zjfrQ~tXMRX21?rU+jutyR0cGUI{P3=vFQ0=8;>q{^OXgS%x+)DwtJIp5gYU;RT``~ zHIh0?X0#KfmFcHga9XbF%w@5N1^PmIx%84bsUzf`{XJ+fJpQ3rjD4A$H`Mn{ajnlk zy6Q^r(?mQ__eiX&&PiCCJENLSI0yF$&FVbS4AB^?G9_=Q6ZO?rTb-7butrjqY}*TC zN7pXiKP|N_sjKftYkiLCs_(f&FM$mv+O%k;8y2LCaX~*cl8Q>&UhiVszjUXr4j6pb zLgIc%26dNOBwFW6%YoDV{BBzsBrN&^c0nQ zvpbbK@a{wRXKHE0E5DYO`_tVJzJ+uhwhl`s8!|)ie5xIp;h}JU5O<1A_o$XZo@lo= z$ZNJ$*2sJ)e7wS^%r6I~o;(kgwcv8gB^`O96_gZXGc&5CkSDsR9n)(0{AMM%)6b^m zBiTW>@vMD!;MZ5$G3U!denVO6FVuHmt@5Z|0H?UVrA1a%AMS9Q17lj-0bG8$&G12$ zL9Ojcu9e(oNh`A;?(kdhs>T7^XZ{iR@zvn)`HF#pk!i2)}C>W9$nfqMaL3yH`%k>rXdsXj8KRr1~ zpEG)}OC_n2L{kb+MT~Io?2y(gCrBV!s#u)w?2B|mgd23UZs|;Y~XFIw48qo_Bw zkLJ5Q8+Ek025USEEpdpN|NMsiOB=Y+MmNmYb6)g1{q#l{&|a4{)th~n-`$Ul(abRm z?_K&!btyZcx-TNV|o+NeL$(?iEf8z$Ds8Vi5bQ%SC`^Vbd* zYU_rjKG=WfaY1 zGv+LNAqiEa#BDZyWAyk1t`*)(qSao?ouxP05&u~TV?GqRm2mnGmacOSikj}<#6D;!E+(VlM##T55|6dKXS19M3q+~iqS1P zL)JI8xpckZ9P3cKBgM6{GLe!_jtOFiPgz6@yB*m4$T9IL3?*0z0!=+}dlXjEeWV+< zhsAk^g9fGG!*lYg^vBTc;h2KFz6v^W|DCYNaiV zTR0^67a2*NDX%2Jv0@qB=W}hHR?u;%qxphoddP=jEO)q4VztxaH4vT$eAWQjHqZxj z&e>ywonvC-dtz4o4sN^7$J-r!x4gpo?oQzRX5oRbowXSR17e`w9(!q_>La-iL`rNj8O|1ZzMg(iVL+Xl|dG66wdEF zOFgMWz6-{8H>&E9R-ZrfmlB0H%3vXS?|kby-EQRRN+M@gfbYpA)}dRPpBV!CHaDcg z7PcqNu-Cwj1*dRv>liz|{>=MPa~Idnar+!!mi0S3S2U(iSy=^gB203dUgT7XI%_pF z_9sQ~D=7WRMPA8-{PcBT1h1E6Op2$Wm|P0-r*Ee=$M8@A!8d}VY)|VrlcO2?liVpzirm zzTmQp6L@~DmaLBau5wmteTYKYCgq6UVg2Bbjn3z9$G?t-Xi;ZgMKfn@O*>6Ote%DX zm`uHy{iH8$r9JgnH}&i&AH?lwc-7@4mX8X=DE9K79ZRW=q2KZ)Nv65der zMd|fdO>SJLFSiiz`Qj!T7z(kNIs{ti`!>Wi)iSV0P4DoYW~L;?=J8^lJN8P!L%XA1 z93e*u!X(KB-;ah_cz#Hw0IG~Mu59xz%WwG{o=zFXn^pG&LN1|chm4Q5iIGa80!yxw zy2LFLgnQhy?`oKm=%GmHO-vgUGj)3{hK#j4N-&+`LUIfQyVW~I)`*#!EIuf;f{k$X zqny$FHRqnLDyDu0NkY2s%Kg1J=xx>Z7?Ys|N4d|F($i>1a{mp@8I}$jx|S@UlUu-aPu+0*3Blx)H#5BZRLbYu#$_eZhEG6SXb`& z)7K7Zw%Fl=kAnp$d!OyA62(1F3%e>~ALRo%64aIy+M7*K#y?W`AT`ZJ17EtiV(G9` zVR^5=a52=Nm;FY;uLjdwu)Zs15J7kdBllW~>y$tc_|84Y9-ypOsqegp7(8;t^jz%N zeqn8~{tsFKo_n;aU9iCJG?AdcnIo%HgU<8YvUT;qbD_*%DV79`N1}8bpX|P${?3(2 z{s5MWsS$Y26bc=dp#EMM@N}^&o`=3alj~WZB_e6K?XvyN4Nk_YgynhGmxwHZ^)cv~ z^T;p=y=yPQ3g%jIqyuw_n92r{zJE!u&p~(&s>S%Jw}DM zv&%V^A2*1(O$2FT@d`8uS5J=50E@u0AuGOmJGE)a$^=&W#c4^fz-ee%IpPbCO>4O? zu)|5r+D!LMoYf!o0n@`^1vO#gT@}l{sX5MC2v~d+h$k<#|Hx3x(N1&IV8Alxux@|; zaOD?kv*9;c`DLyvxe|xwz5L;5&)^+K*dwY@1X;OH$U_UzPSS4tZV7eAXT*0@wr7SF zldQM|$8#}&1!Iy`TabL6Mz_SrkV+_)8JGZ|Y?5|*Pe^(Xlc#G5#nQmX;>|RiF++zN zr%#w>D+ z(Zg5PHGW?G`@ar;*T3QLVUI<;^g9ziNdCAm!sGEgJc58(_-YOsI>9she2Hz<=eO9y z|N75;)>7c3)$QL`Flb5%CHoYI61|(Dt(gFxrR9o!aavx4#@2`H1j%a0|E;6ri9Dix8 zV)*F?`@Ep(mCETqtiM;h?Z$nW@`K<%%0lr2yGwgoyV%{OhCRgIn*LnR-n#7G#>C#n z-rjo89yMGde0LjMA@(-ldmH4@jPSOSLH-%*sb&4W%DoM{#y#gvvi|X%c7I1i=~h3L z9RA4qH`CDS*p-&T6@A8?J!<|Qc`0D8Mt^T1e>j|NhhEYiq=I95g~P@Bj-S7;i<1}q z3~~D&hnAFz)c&uEs+#KmB%V6`Zz>EXrF=aY=^~|k0paa`&c)B!*B$AP@a5RQHh1y@ zIakdL!Jqs&mwQfD0hUhQejwA+$<0qnP3rG$Mn(vfl%3)!mD4Ip>M9zVr=%3s&YV`# z)KWQfT3zathL)<5n!1LHs?@2|s;8CIPM^}$kh;a8e97O*)7{z7+s)Gjw1z|3(9hY$ z+h0oSw5HPO)2GjvhnYv9>S=PX_PrB0pMztq$F4@v0xbV9*k9LlCH-fsSP zPl2|pYMche{JjR3f6D+R`n&i7nt3|;yPR`zM%;1X044oPcvQ}!U8twsU_e6#pUcw6 zOGPw!NfW>Gc27R_K6YyT;acK#_u-?vmPMx5d~@$S%%@4Tx}~9ig(N>pce?nL_WV+| z%mRmOsjh~>=rhVk?h|o^X9@Zr)hg1?%+9opME;N;RGf-37lEZ4qKh?~W_X8>%mn6~7J|K$Ngz7(jgPb|!s+*CBud!N}mlRf9x!n0bxqs(WdFQ=cqn|dd- zhL331*EdjI`;+XkYkduieq#Tz%WBDX(3B?o;2m%Js0D5Mi{h7)UuLa^1!r|a4t&~< z|8j1dcJYI6k=BAfZy5D?hxc6au@6vCUJMvZ)Hsfz!zIgt^&y0M$cJ}AcyWoD& zTATN!b_bUh?cbfOdg=15_^WPq{p9Xto6*%Q@->~g73J$;I~|$DCDo?X$4?4*b|(iL zb_U@h=D$pJ^TTMo{e1&tOA>k)-S-?V*{s|0boAFwy!bNw1N)KsL1#8eyhmb)_r*jl z;tRRVB7FYG3-O6h_2ubF3lkfp)e}-Mql+tn2DML9xbaGvB#*CZYF{!MS+KK~xFvW( zmYaEye6IN82&z7=thD(-q*ZM8ktaMTBjb~*o@^Vn57`v6TBQ+K+fz>tZyQpUI z8q3$)J28(+`)*l;v)AtO)#?%SjSnBuUG&;D6Cd-!V~uL89g2p<{I^>9WwXm0m$)1> zgZLDD4ph?=9I?@(T=$u;n}inN1a-4VZW@M)Sy$ZkxHZtne*f0|!p~n@#5XjR7tK`X z_6H>Dq_Z%$xF&M83Bdz>w%QBXqr@A_Bdo%9o6@z`2RI_<>6Fy-`|uXDKhJF$u$pmG?RJm60Ug&_YyLb zAru+QOk}RiV^@+XGTcH!#zYx1WTxm%#^1TT-@d)Q&-ZzL`+3e?XRou@+Uu;d_WtAZ zdCs|O&)4S4#*Aue+?RvH1&^O^Jh#=nz?Nho(&pc zAG7cZB^$RK_kB`)=DtlA@dqtI3v@?~Zdm(&QBidh8Log8zq!>nmq2ppMx~p<(sco-{oUMm9@4byA++z}F+W=l#myYfY`7!}hpUTy z4^dqxML$^bX5G9ml)BUMM5gc30Z)OheDJ}S+*>gN9rA~)Z-Tw#-vu~3&z*ddCHAi7d*tqb~qW7q;JCy@vOs@9xCt0O?dSs3*2g{wS&JL509r^}ftyU)cHHFDT zBxTC952qkQl8cvo7qiwm5;d((S_e=z1r^kpqk`2N&(D{&R4wFUa#-(PVPg*{{I#D=fs4%H_Pxj^@uTkOxU2P?VSxs6+wr$)*1Cx3AAb)UFAeXUzTZ=+w~ zOU7lF1cR{!p8dQxt*Ek9gQDKmkN(^EL%%t<+@`GogY0Z;5IXDos()&**+cT#=3C0M z_jvxg5j|NuMo0SUoD^rUfz0}{|7_~+f%XLw#4SN}y&oyh42!PzJyx0%RdPMJcX$wa zT4vwvV!>GM$h-uX@k?dHLoH95;P{qxu^{WPnC7Gl5!SFFY`V@9j3gP zcy@OS8I<{g5}BGzQWeGWQcZMlvY>I%c%d&)CDWG;!Sw3h?@CAs~0JA zNcL$VM!-~Q(&M%KfX-nDt3)ECnX$Ky^|g6cRc>HgZ3k=ZIm)8Y3u3$g=p{J=nXTky zZogc6@N*3Kz+BJlM1@^%-F@K?JI1|ru}HU|Gm-=IlOZnb35_fUPeTtnS$8{acXAM} zp+D~;TzQ28NgY=BYJcuzdz?0%7Xm)H#51*f-r+^q^{@Latdm1k1&hhO+xAK0-$D?s z5TWPt>oMbf7B3fDPvc&!-poj;9pZZf54kv zq!vg3zg5$u8;QKi!0YelkyJkLbVKWo?wP-IS{kW^KFHR&R(f70FQS^NVg8H3rl#dl zOODXY=*he$yjx0&jVk4cpvL8P5VaN|M2l`h8*fjuu!N6m4<@e@@;i%evdTl^mS&_f zpY)W^uuEKDn71ULp_&=gefR6T7TxB^W>VoyIse^f5E3X;YK{!@E$&dZ_}2Xsj!;7rg56`DiPW3XBcTJ-~ zYxng;*@G;v^$<&TSZ z&aGw}A5S)ol%Y`!|Z{IK#zYZfWfHdI*YsZGAIf|C_;F-jov2U5M#V^-h0oWpxX zv;tDE$(ek>p<+R4_A+W$8%@xnrPQkdaSnpS^b#V5kY3luqHG&3;oP$~lr=vUfQ__Q z2GUS8ev3Cl*uD7TRRgPYp0Vt@)6>eepyEZDa6j?xP#4dR0&R8D z>b1h#UN!Xf@CQCt0{y91B3s#(qdZ7g1zdWZ;87Xl6we{P(jBj`SgZo+fl*_vg0 zAKVtdiMIC=Nl}~-0t9Dtn@4hC^kP)l4hd4Y!pq<_udst54(+j0~*#3z`l3&Y-=}zdIH5 zvFd^B*>K}3C!?mEbf+67ud`UHvNxgT!$aU`a4to3`~9ip+hxcPCw#E<>#80V9L7=*5SRH!qm#)&@V zs?#FdSN!72a!oxkr-hTSqY*Phi!Q0>(Px&XPdU+JrNGcG1d}1lZO4m-@6hDErFXf2 zTZr7%Pd&w*keUz`FzvwL1F? z2wYp2%;vI`e%ft(=D)>eKX5G*`p$Q;p&4=5kayL`;nZO+-$)8j-v71}jyu>PA3h*CA<|7sk;z?pYxj=U01qgE zdqb`6Sg2YexNrJ}^0f!f{Z%KwQ#n|(W_l}YN3pz;6WR2wA$T7=g6AZ-xnGLvDdRM> z?s9N~`Zv+~&h1GT$4U#YTwE6V${Gau>oBMwqdDol&GyR=#pQLH=1c!ObKs+EOts z+PYTsyfP?;YH=}mR3vkN;v;zk*_;^rUU z47#cAYM;HYlK$bV&FO{hp&O2sxQ<9mOozifV!TaSzPxY-VLUl)=6!92b!(l^swjku z(}R&apY!EINOLJ~^NP`a@C9pUM>fvEZ?9(?&gFJ6k1GX&FJa!8T8O3442~G)pkm3+ zW{!nzLZE3UgA*%-Uuqb4N5JBDi|qw1l*)qH<12DEBYXBHCV;;h?I`jX5h-*ZH#Dqx z_+Ym#F&iCTWz!?Pz!dBVd<_A6~vJ!m0^gR4EF=zc{eb8<+r+pM2#^@JWI#79AA?kh+z;`IA`EF2K1xNfC0+43r^6{Up3%~>Bu>%#Ha*ak zQG?e<#v#Rp#NZuBE%c56h6>xwRL=aF!e8;cKEnH)5Oeu79fNEmO|m)LD5=ph$uQk5 zHfK`W&eX8SF1(n?7+$X_6NLCF1h(6rj>W0&3aSoAQIFvmFT-7jU7YP^^LH2L(tFKc z+K|vb)@H)iDPZ?HEP`G`Oo;@;$ImAIpf#zwO%T#YJ}c)smlBv0>limeLa7q1&sTJX z;w99-6vgj&rqIQhpRA(`81p>W{h3tqoI+TZ;vKR`Pd~mzEI{4uKs3T*FyE$>>8#-l6ztb1VWJS`LV}@8jh6ZfFI42A zD#hVL-gSjNG7fRvb%|oBOwl@L^~S{NyhQ$rf+@`dhm#zGW{g z3_&d%h7nWwVP2-g@Ne9u{ny(xsQd1%kXC&b6&6SjdLwy?(!&J#z3X3`nZ&E zS4}^3RFCI_%&D|#mwppulT7DM-lF6?3U^J8`Sxk7^q36&X}+Hu8lIAgcS&O9!R+y` zcavKna0XJY+j?%ce!g&N=mXleLZl9*sWj`_J(}=jna~SS5QSNu4hWY5719rDcpZHv z!VcCfl7TKrPTWW^5qgF=r%R)1a?6%yf*GxT-lJKB>~UD_D{gxBr~&rV>XD%@#AF#b z&V8=aA;1OVifDPoX;eNZVkMboOQtp43~3bN`?2{Fu|kSJ9!9gYtzIS=BxYZPk}g{? zDL!|Gh6*d(hVCt^V)=tF5E#7&qOlFsIAbIu)i-x|20hcFE-y(X-r+H)+*OJ#ff#1r z4Qlr2`t%9XGJ>mmOS6U~C4~VQr4GKKX38>gpS1#edh$l7d7+w49^99*-x)_w!py<1 zH1$U6LR%F@Vt7&9Ti1yYk+_+MI!zDcTND_h7 zT3>0?(!VS*g-8}`QHfk)46P}Y>^>K5^_I=Pl@2tVR#jV95Zcp`Sm!kQ7TTel1uL(s zZ1B!l&y}^vl#c3>juOkKwY35fm`zSrrTZi#&muIsJo{oOPm~(t@nt3V>&%6!;?*k% z`JhVk3Ievp3R6L-w5cGxFpCz$8`=8!w30bN#DPm+tc9vF%u2)qm)}4n`?xeM$~>c` z=ICFRdPWbozWVao*%Ds$xb;Z`*7p4~tS!2Wz-$GtT7;ap*61RzrRwZyY4N7g#?!j! zo*JEQJ%Lewo4n^#dU4)JohXdVL0 z%j6yTChbUU%_YM(_ZzZe>Jwhh?O6k#7LWNirhJ5h60uSg0xtOrKt=unP=Cjisl^O@ zY`u>lW&nC6qRd3R8BQ&xWaEV+cEJG1`(I?42+#U^0YI*%tGyeT2v6(d92@~;)gR2S z2msnR0kD{^Cl=@F=;{FG*Te$gu%ov>pfvDtcfXEvA%ewFFaXB>0cZo*^D%Y{Okn7y zuBmd_#vLs7AA4#sW3V|CY%UIkf)U~{usA~e*MWdZ{yIO>OG$zy00f(Oz$H*%aRAsx zNWsBKaVUUlTM~ibzoB;+>bH#i4Iv}NrHIk}4HXm7XiZl~AP0J$Zgxam+TPI>>xuJn z^YOI9fo*|Xcclg*$2!`1|J^&%*}2%bA4A#y|6E~|4gQ@bFcxPII0BM!jP)O-SmZxx z5dm!eJK~5Z0Ac|wV8_7_fp0H&8#|n*jjIEWS{4d4d9bV+(2kPu$94Fzm970hhyQJL zB#wN*03Se4O&)9xNVv*T7lAky@MH1!7^&kPaaJfma2yNNkw6|BppHi(0dGWS5=R1c zY#?!L00(?p5+B-sA_Gp3GW;tqc0Qh-KxG|0hNEl|v*758`}L&U-Q0-|jt=6_KRFm8 zZE#F2{EwAEOfB?(EeTauAdkdy|5dIc8vkKNur{vV#6{!Vr}LDrR#rJryTOv)S@8;$WKHf5qwksP?@1w^g1;hZ1g+Nqw)u{gu DwEa6~ literal 0 HcmV?d00001 diff --git a/figs/action-angle.fig/Makefile b/figs/action-angle.fig/Makefile new file mode 100644 index 0000000..0f08a68 --- /dev/null +++ b/figs/action-angle.fig/Makefile @@ -0,0 +1,33 @@ +PROJECTNAME=action-angle + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) +DATS=etas.dat + +all: $(PDFS) + +$(PDFS): $(DATS) + gnuplot $(patsubst %.pdf, %.gnuplot, $@) > $(patsubst %.pdf, %.tikz.tex, $@) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +etas.dat: + python etas.py > $@ + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + rm -f $(addsuffix .tikz.tex, $(PROJECTNAME)) + +clean-dat: + rm -f $(DATS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/action-angle.fig/action-angle.gnuplot b/figs/action-angle.fig/action-angle.gnuplot new file mode 100644 index 0000000..2913583 --- /dev/null +++ b/figs/action-angle.fig/action-angle.gnuplot @@ -0,0 +1,25 @@ +set xlabel "$n$" +set ylabel "$\\Delta_2\\gamma$" norotate + +set ytics 0.6660909, 0.0000001 + +# default output canvas size: 12.5cm x 8.75cm +set term lua tikz size 8,6 standalone + +set key off + + +# set linestyle +set style line 1 linetype rgbcolor "#4169E1" linewidth 2 +set style line 2 linetype rgbcolor "#DC143C" linewidth 2 +set style line 3 linetype rgbcolor "#32CD32" linewidth 2 +set style line 4 linetype rgbcolor "#4B0082" linewidth 2 +set style line 5 linetype rgbcolor "#DAA520" linewidth 2 +set style line 6 linetype rgbcolor "#555500" linewidth 2 + +set pointsize 1 + +plot \ + "etas.dat" using ($0*2):1 ls 1 ,\ + "etas.dat" using ($0*2+1):2 ls 2 + diff --git a/figs/action-angle.fig/etas.dat b/figs/action-angle.fig/etas.dat new file mode 100644 index 0000000..911f526 --- /dev/null +++ b/figs/action-angle.fig/etas.dat @@ -0,0 +1,49 @@ +0.6660909569512501 0.6660909569517841 +0.6660909569512068 0.6660909568572682 +0.666090881906312 0.6660908743475806 +0.6660909492953624 0.6660909569530551 +0.6660909569551257 0.666090956950617 +0.6660909569498561 0.6660909569513911 +0.6660909569508098 0.666090957737135 +0.6660908737928897 0.6660908734660769 +0.6660909574099513 0.6660909569513342 +0.6660909569510589 0.6660909569494657 +0.666090956949009 0.6660909569517685 +0.6660909569506912 0.6660909541449044 +0.6660908740414464 0.6660908770567131 +0.6660909571597132 0.6660909569500939 +0.666090956949887 0.666090956948812 +0.6660909569483953 0.666090956957927 +0.6660909270003836 0.6660909609769751 +0.6660909042589079 0.6660908702167302 +0.666090956890927 0.666090956948973 +0.6660909569482991 0.666090956948177 +0.666090956947643 0.6660909569220005 +0.666090854453731 0.6660909700231973 +0.6660909767366049 0.6660908612031291 +0.6660909569565692 0.666090956947536 +0.6660909569466511 0.6660909569475739 +0.6660909569469959 0.6660909569534238 +0.6660908799755063 0.6660908742427215 +0.6660909512195086 0.6660909569483087 +0.666090956948585 0.6660909569461111 +0.6660909569454749 0.6660909569473561 +0.6660909569467641 0.6660909579964863 +0.6660908731876072 0.6660908731982366 +0.666090958006647 0.666090956946789 +0.666090956946648 0.666090956945196 +0.666090956944789 0.6660909569490353 +0.666090956946594 0.6660909496555434 +0.6660908743084164 0.6660908815363281 +0.6660909568839747 0.6660909569457019 +0.6660909569453359 0.6660909569445963 +0.6660909569441398 0.6660909569542799 +0.6660908386690396 0.6660909880746821 +0.6660909925361512 0.6660908431096189 +0.6660909569319037 0.666090956944476 +0.6660909569436897 0.6660909569439311 +0.6660909569433824 0.6660909568669846 +0.666090876905717 0.6660908743223477 +0.6660909542783731 0.666090956945574 +0.6660909569495623 0.6660909569429838 +0.6660909569421549 0.6660909569434601 diff --git a/figs/ellipse.fig/Makefile b/figs/ellipse.fig/Makefile new file mode 100644 index 0000000..a15ccff --- /dev/null +++ b/figs/ellipse.fig/Makefile @@ -0,0 +1,24 @@ +PROJECTNAME=ellipse + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=$(addsuffix .tikz.tex, $(PROJECTNAME)) + +all: $(PDFS) + +$(PDFS): $(TEXS) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/ellipse.fig/ellipse.tikz.tex b/figs/ellipse.fig/ellipse.tikz.tex new file mode 100644 index 0000000..96b6bff --- /dev/null +++ b/figs/ellipse.fig/ellipse.tikz.tex @@ -0,0 +1,86 @@ +\documentclass{standalone} + +\usepackage{xcolor} + +\definecolor{darkgreen}{HTML}{329D32} + +\usepackage{tikz} + +\begin{document} +\begin{tikzpicture}[scale=2.5] + +\draw[color=blue](0,0)ellipse(1.5 and 1.0); + +% angle of L +\def\tt{15} +% foci ae=sqrt(1-b^2/a^2) +\def\ae{1.118} +% collision point +\def\Px{0.5} +% Py=sqrt(1-Px^2/a^2) +\def\Py{0.943} +% slope Px/(a^2*sqrt(1-Px^2/a^2)) +\def\slope{0.236} +% phi: atan(\Py/(\ae+\Px)) +\def\ph{15.2} +% phi': pi-atan(\Py/(\ae-\Px)) +\def\php{108.2} +% lambda: pi-atan(\slope) +\def\lam{151.7} +% psi: pi+phi-lambda +\def\ps{43.5} +% h=sin(\ph)*sqrt((\Px+\ae)^2+Py^2) +\def\h{0.491} + +% foci +\path(-\ae,0)coordinate(F1); +\path(\ae,0)coordinate(F2); +% collision point +\path(\Px,\Py)coordinate(P); +%Q +\path(F1)++(90+\tt:\h)coordinate(Q); + +\draw[densely dotted](F1)++(\tt:3)--++(180+\tt:3.5); +\draw[densely dotted](F2)++(\tt:0.7)--++(180+\tt:3.5); +\draw[line width=1](P)++(\tt:1)--++(180+\tt:3.5); + +\draw[color=gray](F1)--(P)--(F2); + +\draw[color=darkgreen,dashed](P)--++(1,-\slope); +\draw[color=darkgreen,dashed](P)--++(-1,\slope); + +\draw[rotate=\tt](F1)++(0.4,0)arc(0:\ph:0.4); +\draw(F1)++(\tt+\ph/2:0.5)node{$\varphi$}; +\draw[rotate=\tt](F2)++(0.15,0)arc(0:\php:0.15); +\draw(F2)++(\tt+\php/2:0.25)node{$\varphi'$}; + +\draw[rotate=\tt](P)++(0.1,0)arc(0:\lam:0.1); +\draw(P)++(\tt+\lam/2:0.2)node{$\lambda$}; + +\draw[rotate=\tt+\lam](P)++(0.2,0)arc(0:\ps:0.2); +\draw(P)++(\tt+\lam+\ps/2-7:0.3)node{$\psi$}; +\draw[rotate=180+\tt+\lam](P)++(0.2,0)arc(0:-\ps:0.2); +\draw(P)++(\tt+180+\lam-\ps/2:0.3)node{$\psi$}; + + +\fill[color=red](F1)circle(0.03); +\draw(F1)++(0,-0.15)node{$O$}; +\fill(F2)circle(0.03); +\draw(F2)++(0,-0.15)node{$O'$}; + +\fill(P)circle(0.03); +\draw(P)++(\lam+\ps+90-\ps:0.15)node{$P$}; + +\fill(0,0)circle(0.03); +\draw(0,-0.15)node{$C$}; + +\fill(Q)circle(0.03); +\draw(Q)++(0,-0.15)node{$Q$}; + +\draw[dotted,<->](F1)++(\tt:-0.5)--++(90+\tt:\h); +\draw(F1)++(\tt:-0.6)++(90+\tt:\h/2)node{$h$}; + + +\end{tikzpicture} +\end{document} + diff --git a/figs/trajectory.fig/Makefile b/figs/trajectory.fig/Makefile new file mode 100644 index 0000000..3af70f0 --- /dev/null +++ b/figs/trajectory.fig/Makefile @@ -0,0 +1,30 @@ +PROJECTNAME=trajectory + +SIMPLEQ=~/Work/Research/2018+bose_gas/cmp/simpleq + +PDFS=$(addsuffix .pdf, $(PROJECTNAME)) +TEXS=trajectory.tikz.tex + +all: $(PDFS) + +$(PDFS): $(TEXS) + pdflatex -jobname $(basename $@) -file-line-error $(patsubst %.pdf, %.tikz.tex, $@) + +trajectory.tikz.tex: + python trajectory.py > $@ + + +install: $(PDFS) + cp $^ $(INSTALLDIR)/ + +clean-aux: + rm -f $(addsuffix .aux, $(PROJECTNAME)) + rm -f $(addsuffix .log, $(PROJECTNAME)) + +clean-dat: + rm -f $(TEXS) + +clean-tex: + rm -f $(PDFS) + +clean: clean-aux clean-tex diff --git a/figs/trajectory.fig/trajectory.tikz.tex b/figs/trajectory.fig/trajectory.tikz.tex new file mode 100644 index 0000000..8854a5c --- /dev/null +++ b/figs/trajectory.fig/trajectory.tikz.tex @@ -0,0 +1,26 @@ +\documentclass{standalone} +\usepackage{tikz} +\begin{document} +\begin{tikzpicture} +\draw[rotate=263.6206297915572,line width=1,dashed](-4.47213595499958,0) ellipse (5.0 and 2.2360679774997902); +\begin{scope} +\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0); +\draw[rotate=263.6206297915572,line width=2, color=blue](-4.47213595499958,0) ellipse (5.0 and 2.2360679774997902); +\end{scope} +\draw[rotate=231.45599405039505,line width=1,dashed](-3.725076664340498,0) ellipse (5.0 and 3.3352366999638683); +\begin{scope} +\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0); +\draw[rotate=231.45599405039505,line width=2, color=blue](-3.725076664340498,0) ellipse (5.0 and 3.3352366999638683); +\end{scope} +\draw[rotate=172.42531416885245,line width=1,dashed](-1.2502080963147824,0) ellipse (5.0 and 4.841175447751193); +\begin{scope} +\clip(-5,2)--++(15,0)--++(0,8)--++(-15,0); +\draw[rotate=172.42531416885245,line width=2, color=blue](-1.2502080963147824,0) ellipse (5.0 and 4.841175447751193); +\end{scope} +\draw[line width=2](-5,2)--++(15,0); +\fill[color=red](0,0)circle(0.2); +\draw(0.5,0)node{\huge$O$}; +\draw[densely dotted,<->,line width=1](0,0)--++(0,2); +\draw(0.3,1)node{\huge$h$}; +\end{tikzpicture} +\end{document}