Instructions: The following paper, which will appear in the Bulletin
of the Amer. Math. Soc., contains some graphics and therefore consists
of two files one right after the other. The first is written in plain
Tex and contains the full text of the paper; it ends with the line
"\bye".  The other file is written in postscript and contains five
figures; it begins with the line "%!PS-Adobe-2.0". 

\magnification 1200
\let\nd\noindent %						
\def\NL{\hfill\break}%						
\def\qed{\hbox{\hskip 6pt\vrule width6pt height7pt depth1pt \hskip1pt}}
\def\natural{{\rm I\kern-.18em N}}
\def\N{{\bf N}}
\def\integer{{\rm Z\kern-.32em Z}}
\def\Z{{\bf Z}}
\def\real{{\rm I\kern-.2em R}}
\def\R{{\bf R}}
\def\complex{\kern.1em{\raise.47ex\hbox{
	    $\scriptscriptstyle |$}}\kern-.40em{\rm C}}
\def\C{{\bf C}}
\def\vs#1 {\vskip#1truein}
\def\hs#1 {\hskip#1truein}
\hskip 2.5truein{\date}\vskip.5truein}
  \hsize=5.8truein	\hoffset=.25truein 
  \vsize=8.8truein	%\voffset=1truein
  \parskip=0 pt		\parindent=20pt
  \overfullrule=0pt	\lineskip=0pt	\lineskiplimit=0pt
  \hbadness=10000 \vbadness=10000 	

\pageno=0
\baselineskip=15.5pt
\font\bit=cmbxti10
\footline={\ifnum\pageno=0
\tenrm\hfil *Research supported in part by Texas ARP Grant 003658-113\hfil
\else\hfil\tenrm \folio\hfil \fi}
\hbox{}
\vskip 1truein\centerline{{\bf SYMMETRY OF TILINGS OF THE PLANE}*}
\vskip .5truein\centerline{by}
\centerline{Charles Radin}
\vskip .2truein\centerline{Mathematics Department}
\centerline{University of Texas}
\centerline{Austin, TX\ \ 78712}
\vs.5
\centerline{{\bf Abstract}}
\vs.1 \nd
We discuss two new results on tilings of the plane.
In the first, we give sufficient conditions for the tilings 
associated with an inflation rule to be uniquely ergodic
under translations, the conditions holding for the pinwheel
inflation rule. In the second result we prove there are
matching rules for the pinwheel inflation rule, making the
system the first known to have complete rotational symmetry.
\vs1
\centerline{1991 Classification: 52C20, 58F11, 47A35}
\vfill \eject
\pageno=1
We consider tilings of the Euclidean plane, $E^2$, by 
(orientation preserving) congruent copies
of a fixed finite set of prototiles. Prototiles are topological disks
in the plane satisfying some mild restrictions on their shapes, as
detailed below. Congruent copies of prototiles are called tiles, and a
tiling is simply an unordered collection of tiles whose union is the
plane, and in which each pair of tiles has disjoint interiors.

We are concerned here with two constructions associated with a fixed
finite set $S=\{P_j\}$ of prototiles, the most important of which is
the set $X(S)$ of all tilings by tiles from $S$. In particular, we are
interested in understanding the purest cases, in which all the
tilings in $X(S)$ are ``essentially the same''; we will define this
precisely below. Two examples are exhibited in Figures 1 and 2, both with
two prototiles; in Figure 1, $S = S_K$ only produces a
checkerboard-like tiling (and all congruences), and in Figure 2,
$S=S_P$ produces the well-known tilings of Penrose [3,4,6].

Tilings like those of Penrose are not usually invariant under any congruence
of the plane (other than the identity), so to analyze their symmetries
we introduce some elementary ergodic theory, and also another basic
construction which can sometimes be associated with a prototile set
$S$, the set $X^F(S)$ of tilings defined by an ``inflation function''
$F$.  An ``inflation rule'' for $S$, if it exists, consists of a
dilation $D_F$ of $E^2$ by some factor $\lambda_F < 1$, and a finite set
$\{C_{jk}\}$ of congruences of $E^2$, such that for each $P_j \in S$ we
have: $$P_j =
\bigcup_kC_{jk}D_FP_{n_k} \eqno 1)$$

\nd where the elements of each union have pairwise disjoint interiors.
The inflation function $F$ associated with such a rule is defined on
tiles (and then sets of tiles), with sets of tiles as values, as
follows. If the tile $P$ is ``of tile-type $j$'', that is, $P = CP_j$
where $P_j\in S$ and $C$ is a congruence, then 

$$F: P\longrightarrow F(P) \equiv\{D_F^{-1}CC_{jk}D_FP_{n_k}\} \eqno
2)$$

\nd (Intuitively, $P$ can be replaced by a set of ``small-size tiles'' 
by 1), which are then expanded in 2) to original size by the inverse
of the dilation. This process can obviously be applied to any collection of
tiles, for example a tiling, and can thus be iterated.)

The tilings $X^F(S)$ associated with the inflation function $F$ are
then defined as those tilings $T$ such that each finite subcollection
of $T$ is congruent to a subcollection of a set of tiles of the form
$F^r(P)$ for some prototile $P$ and integer $r\geq 1$. 

The construction $X^F(S)$ is mainly of interest when $F$ is a
homeomorphism on it, since then $F$ defines a natural representation
for the dilation $D_F$ as a map of $X^F(S)$ onto itself. 
(It is ``natural'' in that it extends the representation of the
congruences to a larger subgroup of the conformal group.)
One can then
consider ``symmetry'' with respect to this hierarchical action, as we
do below.  It is noteworthy that dilational (and rotational)
symmetry is manifested not by the invariance of tilings themselves,
but rather of measures on tilings; that is, the action of the dilation
(and rotations) is lifted to the set of (translation invariant Borel
probability) measures on the tilings, and invariance is sought in this
set. This is one reason for using the machinary of ergodic theory in
the analysis of tilings.
\vs.1
We note the following examples of inflation rules.
\vs0 \nd 
1) The square inflation rule $S$, given in Figure 3 -- $S_S$ has
one element;
\vs0 \nd 
2) The Robinson inflation rule $R$, given in Figure 4 -- $S_R$ has
two elements;
\vs0 \nd 
3) Conway's ``pinwheel'' inflation rule $C$, given in Figure 5 --
$S_C$ has two elements.
\vs.1

(The inflation rule is one-to-one -- and in fact a
homeomorphism -- on the associated space of tilings for the Robinson
and pinwheel rules, but not the square rule, for which it is
four-to-one.) Each $X(S)$ (and $X^F(S)$) carries a natural metric
structure (indicated below), and defines a dynamical system with $\real^2$
action, where $\real^2$ acts by simultaneous translation of the tiles in a
tiling.  Using this language we note the nonobvious fact [4] that the
dynamical systems associated with $X^R(S_R)$ and $X(S_P)$ are
topologically conjugate: there is a homeomorphism between $X^R(S_R)$
and $X(S_P)$ which intertwines the translations. It is immediate that
those associated with $X^S(S_S)$ and $X(S_K)$ are topologically
conjugate.

We now outline our results. Our first result gives conditions
sufficient for the dynamical system associated with some $X^F(S)$ to
be uniquely ergodic, the conditions being satisfied for example by
$X^C(S_C)$. Unique ergodicity of a dynamical system, namely the
property that there is one and only one Borel probability measure on
the space invariant under the group action -- the space here being
$X^F(S)$ and the group being $\real^2$ -- is useful because it means
the dynamical system has a natural measure associated with it. (From
Birkhoff's pointwise ergodic theorem [1], this is equivalent to having
all the tilings being statistically identical [6].) The Penrose (and
therefore Robinson) systems are not quite uniquely ergodic. However,
each ergodic measure is invariant under rotation by $2\pi/10$ about
any point of $E^2$, and the different measures are merely rotations of one
another, by an angle in $(0,2\pi/10)$. So modulo this rotation there
is a unique invariant measure for these systems.  Furthermore, the
ergodic measures for the Robinson system are each invariant under the
dilation associated with Figure 4.

Our second result is the construction of a finite prototile set
$\tilde S$ such that the dynamical system associated with $X(\tilde
S)$ is uniquely ergodic, and metrically conjugate with that associated
with $X^C(S_C)$.  (Metric conjugacy -- the existence of a measurable
bijection modulo measure zero, which intertwines the dynamics -- is
somewhat weaker than topological conjugacy.) In other words,
$X^C(S_C)$ plays a role for $X(\tilde S)$ similar to the one
$X^R(S_R)$ plays for $X(S_P)$; in each case the former defines for the
latter a hierarchical symmetry which for example is manifested in a
symmetry in the spectrum (discussed below) of the latter in each
ergodic component. The pinwheel and Penrose systems differ
significantly in their rotational symmetry: the Penrose system has the
ten-fold rotational symmetry of each of its ergodic components, while
the pinwheel system, being uniquely ergodic, has full rotational
symmetry.  (The actions on a space of tilings of congruences of
$E^2$, and of a dilation $D_F$ -- assuming $F$ is a homeomorphism on
$X^F(S)$ of course -- can be lifted to actions on the set of invariant
measures of the dynamical system, so if a dynamical system is uniquely
ergodic its measure must be invariant under such actions.) The
symmetries of the spectrum of tiling systems, in particular the
Penrose tilings, have been of major importance in their connection
with theories of the structure matter, for example quasicrystals [6].

We now state our assumptions and results more fully.
We assume that the prototiles in $S$
satisfy the following conditions, besides being homeomorphs
of the closed unit disk:
\vs0 \hs.25 $\bullet$\ $X(S)$ is nonempty; 
\vs0 \hs.25 $\bullet$\ in the tilings being considered, the boundary of 
each tile can be covered by tiles in only finitely many ways, up to
congruence;
\vs0 \hs.25 $\bullet$\ each prototile $P$ has small surface to volume ratio in 
the sense that:

$${{area\ \{x\in tP\ :\ \Vert x-y\Vert \le 1\ \hbox{ for some }y\in
\partial(tP)\}}\over {area\ \{tP\}}}\longrightarrow 0 \eqno 3)$$

\nd as the expansion factor $t \to \infty$.

A metrizable topology is put on the space $X(S)$ or $X^F(S)$ of
tilings with the following neighborhood basis
$\{N_T(\epsilon)\,:\,\epsilon > 0\}$ of each tiling $T$: $N_T(\epsilon)$
consists of all tilings $T'$ such that within the circle in the plane
centered at the origin and of radius $1/\epsilon$, wherever $T$ has a
tile so does $T'$, within distance $\epsilon$ in the Hausdorff metric
on compact sets. It is known that $X(S)$ and $X^F(S)$ are compact in
this topology, and that the actions of translations $\real^2$ on
$X(S)$ and $X^F(S)$ are continuous [9]. Given the above then, our
results are the following.

{\bf Theorem 1. [8]}\ \ Assume given some $X^F(S)$ as above, and assume
there is some $r\ge 1$ such that:
\vs0 \hs.5
a)\ \ for each prototile $P$, $F^r(P)$ contains tiles of every tile-type;
\vs0 \hs.5
b)\ \ for some prototile $P,$ $F^r(P)$ contains two tiles $P',\ P''$ of the 
same 
\vs0 \hs.75 tile-type, whose relative rotation is irrational with 
respect to $\pi$.
\vs0 \nd
Then $(X^F(S),\real^2)$ is uniquely ergodic.
\vs.05
{\bf Theorem 2. [7]}\ \ There is a finite set $\tilde S$ of prototiles
such that $(X(\tilde S), \real^2)$ is uniquely ergodic and metrically
conjugate to $(X^C(S_C), \real^2)$.

\vs.1
We note that $X^C(S_C)$ satisfies the hypotheses of Theorem 1 (with
$r=2$), and close with a few comments on these results of a more
technical nature.

Much of the symbolic (substitution) dynamics of hierarchical
structures makes use of a square matrix $A$, of which $A_{jk}$ is the
number of times tile-type $j$ is associated, by the inflation function
$F$, with tile-type $k$. To prove Theorem 1 this matrix is generalized to
a family $A[m]$ of matrices for which $A[m]_{jk}$ is the sum
$\sum_ne^{ima_n(j,k)}$, over all tiles of type $j$ contained in
$F(P_k)$, where $a_n(j,k)$ is the angle of rotation of the tile
compared to the defining prototile of its type.

This generalization allows us to keep track of rotational information,
and we prove unique ergodicity by using Weyl's criterion on uniform
distribution as it is used in elementary treatments of rotations of
the circle [1]. The natural way in which the important matrix $A$
generalizes is evidence to us that much of symbolic substitution
dynamics can be generalized to tilings $X^F(S)$.

The import of Theorem 2 is also, in part, its relation to symbolic
dynamics. The systems $X(S)$ are natural generalizations of symbolic
systems of finite type with $\integer^2$ action; for example, they play
roughly the same role for statistical mechanics in Euclidean spaces as
systems of finite type do for lattice gas statistical mechanics [6]. Of
particular importance, for example in theories of material structure,
are the spectra of such dynamical systems. To see why the symmetries
of the dynamical system are relevant for this, note that if a rotation or
dilation $W$ preserves the ergodic probability measure $m$ of the
dynamical system, then, as shown in [8], the spectral projection
$E_{\Delta}$ of the translations, associated to a set $\Delta\subset \real^2$, is unitarily
equivalent, in the Hilbert space defined by $m$, to $E_{W(\Delta )}$.
(See [2] for explicit connections between dynamics spectra and X-ray
spectra of scatterers.) 

Another aspect of Theorem 2 is its relation to the work of S.\ Mozes,
who proved in [5] that a rather general class of symbolic substitution
dynamical systems with $\integer^2$ action are of finite type; our example
is, we hope, the first step in a parallel theorem for tilings of the
plane.
\vs.2
\centerline{References}
\vs.1 \nd
[1] I.\ Cornfeld, S.\ Fomin and Ya.\ Sinai, {\it Ergodic Theory}, 
Springer-Verlag, New York, 1982.
\vskip.1truein\noindent
[2] S.\ Dworkin, Spectral theory and X-ray diffraction, {\it J.\ Math.\ Phys.},
to appear.
\vskip.1truein\noindent
[3] M.\ Gardner, Extraordinary nonperiodic tiling that enriches the theory
of tiles, {\it Sci.\ Amer.}\ January 1977, 116-119.
\vskip.1truein\noindent
[4] B.\ Gr\"unbaum and G.C.\ Shephard, {\it Tilings and Patterns},
Freeman, New York, 1986.
\vskip.1truein\noindent
[5] S.\ Mozes, Tilings, substitution systems and dynamical systems
generated by them, {\it J.\ d'Analyse Math.}\ 53 (1989), 139-186.
\vskip.1truein\noindent
[6] C.\ Radin, Global order from local sources,
{\it Bull.\ Amer.\ Math.\ Soc.}\ {\bf 25} (1991), 335--364.
\vskip.1truein\noindent
[7] C.\ Radin, The pinwheel tilings of the plane,
{\it Annals of Math.}, to appear.
\vskip.1truein\noindent
[8] C.\ Radin, Space tilings and substitutions, University of
Texas preprint.
\vskip.1truein\noindent
[9] C.\ Radin and M.\ Wolff, Space tilings and local isomorphism,
{\it Geometriae Dedicata}\ {\bf 42} (1992), 355-360.





\bye
%!PS-Adobe-2.0
%%Creator: dvips 5.47 Copyright 1986-91 Radical Eye Software
%%Title: bullfigs4.dvi
%%Pages: 1 1
%%BoundingBox: 0 0 612 792
%%EndComments
%%BeginProcSet: tex.pro


/TeXDict 200 dict def   % define a working dictionary
TeXDict begin           % start using it.
/N /def load def
/B { bind def } N
/S /exch load def
/X { S N } B
/TR /translate load N



/isls false N
/vsize 10 N

/@rigin                 % -xps -yps @rigin -   establishes dvips conventions
  { isls { [ 0 1 -1 0 0 0 ] concat } if
    72 Resolution div 72 VResolution div neg scale
    Resolution VResolution vsize neg mul TR
    matrix currentmatrix
    dup dup 4 get round 4 exch put
    dup dup 5 get round 5 exch put
    setmatrix } N


/@letter { /vsize 10 N } B

/@landscape { /isls true N /vsize -1 N } B

/@a4 { /vsize 10.6929133858 N } B

/@a3 { /vsize 15.5531 N } B

/@ledger { /vsize 16 N } B

/@legal { /vsize 13 N } B

/@manualfeed
   { statusdict /manualfeed true put
   } B

        % n @copies -   set number of copies
/@copies
   { /#copies X
   } B


/FMat [1 0 0 -1 0 0] N
/FBB [0 0 0 0] N

/nn 0 N /IE 0 N /ctr 0 N
/df-tail       % id numcc maxcc df-tail -- initialize a new font dictionary
  {
    /nn 8 dict N              % allocate new font dictionary
    nn begin
        /FontType 3 N
	/FontMatrix fntrx N
	/FontBBox FBB N
        string /base X
        array /BitMaps X
        /BuildChar {CharBuilder} N
        /Encoding IE N
        end
    dup { /foo setfont }          %  dummy macro to be filled in
       2 array copy cvx N         %  have to allocate a new one
    load                          %  now we change it
       0 nn put
    /ctr 0 N                      %  go, count, and etc.
    [                               %  start next char definition
  } B
/df {
   /sf 1 N
   /fntrx FMat N
   df-tail
} B
/dfs { div /sf X
   /fntrx [ sf 0 0 sf neg 0 0 ] N
   df-tail
} B

/E { pop nn dup definefont setfont } B


/ch-width {ch-data dup length 5 sub get} B % the number of pixels across
/ch-height {ch-data dup length 4 sub get} B % the number of pixels tall
/ch-xoff  {128 ch-data dup length 3 sub get sub} B % num pixels right of origin
/ch-yoff  {ch-data dup length 2 sub get 127 sub} B % number of pixels below origin
/ch-dx  {ch-data dup length 1 sub get} B     % number of pixels to next character
/ch-image {ch-data dup type /stringtype ne
      { ctr get /ctr ctr 1 add N } if
   } B                        % the hex string image, or array of same
/id 0 N /rw 0 N /rc 0 N /gp 0 N /cp 0 N /G 0 N /sf 0 N

/CharBuilder    % fontdict ch Charbuilder -     -- image one character
     {save 3 1 roll S dup /base get 2 index get S /BitMaps get S get
      /ch-data X pop
      /ctr 0 N
      ch-dx 0 ch-xoff ch-yoff ch-height sub
      ch-xoff ch-width add ch-yoff
      setcachedevice
      ch-width ch-height true
      [1 0 0 -1 -.1 ch-xoff sub ch-yoff .1 add]
      {ch-image} imagemask
      restore
  } B

               % in the following, the font-cacheing mechanism requires that
                % a name unique in the particular font be generated

/D            % char-data ch D -    -- define character bitmap in current font
  { /cc X
    dup type /stringtype ne {]} if
    nn /base get cc ctr put
    nn /BitMaps get S ctr S
    sf 1 ne {
       dup dup length 1 sub dup 2 index S get sf div put
    } if
    put
    /ctr ctr 1 add N
  } B

/I            % a faster D for when the next char follows immediately
  { cc 1 add D } B

/bop           % bop -              -- begin a brand new page
  {
    userdict /bop-hook known { bop-hook } if
    /SI save N
    @rigin
    userdict /offset-hook known { offset-hook } if
    0 0 moveto
  } N

/eop           % - eop -              -- end a page
  { % eop-aux  % -- to observe VM usage
    clear SI restore
    showpage
    userdict /eop-hook known { eop-hook } if
  } N

/@start         % - @start -            -- start everything
  {
    userdict /start-hook known { start-hook } if
    /VResolution X
    /Resolution X
    1000 div /DVImag X
    /IE 256 array N
    0 1 255 {IE S 1 string dup 0 3 index put cvn put} for
  } N

/p /show load N        %  the main character setting routine

/RMat [ 1 0 0 -1 0 0 ] N % things we need for rules
/BDot 260 string N
/rulex 0 N /ruley 0 N
/v {                   % can't use ...fill; it makes rules too big
   /ruley X /rulex X
   V
} B
/V
statusdict begin /product where
   { pop product dup length 7 ge
      { 0 7 getinterval (Display) eq } { pop false } ifelse }
   { false } ifelse end
{ {
   gsave
      TR -.1 -.1 TR 1 1 scale rulex ruley
      false RMat { BDot } imagemask
   grestore
} }
{ {
   gsave
      TR -.1 -.1 TR rulex ruley scale 1 1
      false RMat { BDot } imagemask
   grestore
} } ifelse B
/a { moveto } B    % absolute positioning
/delta 0 N         % we need a variable to hold space moves
/tail { dup /delta X 0 rmoveto } B
/M { S p delta add tail } B
/b { S p tail } B      % show and tail!
/c { -4 M } B
/d { -3 M } B
/e { -2 M } B
/f { -1 M } B
/g { 0 M } B
/h { 1 M } B
/i { 2 M } B
/j { 3 M } B
/k { 4 M } B
/w { 0 rmoveto } B
/l { p -4 w } B
/m { p -3 w } B
/n { p -2 w } B
/o { p -1 w } B
/q { p 1 w } B
/r { p 2 w } B
/s { p 3 w } B
/t { p 4 w } B
/x { 0 S rmoveto } B
/y { 3 2 roll p a } B
/bos { /SS save N } B
/eos { clear SS restore } B

end  % revert to previous dictionary
%%EndProcSet
%%BeginProcSet: special.pro
TeXDict begin /SDict 200 dict N SDict begin /@SpecialDefaults{/hs 612 N /vs
792 N /ho 0 N /vo 0 N /hsc 1 N /vsc 1 N /ang 0 N /CLIP false N /BBcalc false N
/p 3 def}B /@scaleunit 100 N /@hscale{@scaleunit div /hsc X}B /@vscale{
@scaleunit div /vsc X}B /@hsize{/hs X /CLIP true N}B /@vsize{/vs X /CLIP true
N}B /@hoffset{/ho X}B /@voffset{/vo X}B /@angle{/ang X}B /@rwi{10 div /rwi X}
B /@llx{/llx X}B /@lly{/lly X}B /@urx{/urx X}B /@ury{/ury X /BBcalc true N}B
/magscale true def end /@MacSetUp{userdict /md known{userdict /md get type
/dicttype eq{md begin /letter{}N /note{}N /legal{}N /od{txpose 1 0 mtx
defaultmatrix dtransform S atan/pa X newpath clippath mark{transform{
itransform moveto}}{transform{itransform lineto}}{6 -2 roll transform 6 -2
roll transform 6 -2 roll transform{itransform 6 2 roll itransform 6 2 roll
itransform 6 2 roll curveto}}{{closepath}}pathforall newpath counttomark array
astore /gc xdf pop ct 39 0 put 10 fz 0 fs 2 F/|______Courier fnt invertflag{
PaintBlack}if}N /txpose{pxs pys scale ppr aload pop por{noflips{pop S neg S TR
pop 1 -1 scale}if xflip yflip and{pop S neg S TR 180 rotate 1 -1 scale ppr 3
get ppr 1 get neg sub neg ppr 2 get ppr 0 get neg sub neg TR}if xflip yflip
not and{pop S neg S TR pop 180 rotate ppr 3 get ppr 1 get neg sub neg 0 TR}if
yflip xflip not and{ppr 1 get neg ppr 0 get neg TR}if}{noflips{TR pop pop 270
rotate 1 -1 scale}if xflip yflip and{TR pop pop 90 rotate 1 -1 scale ppr 3 get
ppr 1 get neg sub neg ppr 2 get ppr 0 get neg sub neg TR}if xflip yflip not
and{TR pop pop 90 rotate ppr 3 get ppr 1 get neg sub neg 0 TR}if yflip xflip
not and{TR pop pop 270 rotate ppr 2 get ppr 0 get neg sub neg 0 S TR}if}
ifelse scaleby96{ppr aload pop 4 -1 roll add 2 div 3 1 roll add 2 div 2 copy
TR .96 dup scale neg S neg S TR}if}N /cp{pop pop showpage pm restore}N end}if}
if}N /normalscale{Resolution 72 div VResolution 72 div neg scale magscale{
DVImag dup scale}if}N /psfts{S 65536 div N}N /startTexFig{/psf$SavedState save
N userdict maxlength dict begin /magscale false def normalscale currentpoint
TR /psf$ury psfts /psf$urx psfts /psf$lly psfts /psf$llx psfts /psf$y psfts
/psf$x psfts currentpoint /psf$cy X /psf$cx X /psf$sx psf$x psf$urx psf$llx
sub div N /psf$sy psf$y psf$ury psf$lly sub div N psf$sx psf$sy scale psf$cx
psf$sx div psf$llx sub psf$cy psf$sy div psf$ury sub TR /showpage{}N
/erasepage{}N /copypage{}N /p 3 def @MacSetUp}N /doclip{psf$llx psf$lly
psf$urx psf$ury currentpoint 6 2 roll newpath 4 copy 4 2 roll moveto 6 -1 roll
S lineto S lineto S lineto closepath clip newpath moveto}N /endTexFig{end
psf$SavedState restore}N /@beginspecial{SDict begin /SpecialSave save N gsave
normalscale currentpoint TR @SpecialDefaults}N /@setspecial{CLIP{newpath 0 0
moveto hs 0 rlineto 0 vs rlineto hs neg 0 rlineto closepath clip}if ho vo TR
hsc vsc scale ang rotate BBcalc{rwi urx llx sub div dup scale llx neg lly neg
TR}if /showpage{}N /erasepage{}N /copypage{}N newpath}N /@endspecial{grestore
clear SpecialSave restore end}N /@defspecial{SDict begin}N /@fedspecial{end}B
/li{lineto}B /rl{rlineto}B /rc{rcurveto}B /np{/SaveX currentpoint /SaveY X N 1
setlinecap newpath}N /st{stroke SaveX SaveY moveto}N /fil{fill SaveX SaveY
moveto}N /ellipse{/endangle X /startangle X /yrad X /xrad X /savematrix matrix
currentmatrix N TR xrad yrad scale 0 0 1 startangle endangle arc savematrix
setmatrix}N end
%%EndProcSet
TeXDict begin 1200 300 300 @start /Fa 11 118 df<004001C007C0FFC0F9C001C001C001
C001C001C001C001C001C001C001C001C001C001C001C001C001C001C001C001C001C001C001C0
01C001C001C001C07FFF7FFF10217CA019>49 D<03F8000FFE00181F802007C04003C04001E0F8
01E0FC01F0FC00F0FC00F07800F03001F00001E00001E00003E00003C0000780000700000E0000
1C0000380000700000E0000080000100000200100400100800101000302000207FFFE0FFFFE0FF
FFE014217EA019>I<01F8000FFE001C0F803003C03803C07C03E07C01E07C01E03C03E01803E0
0003C0000380000780000E00001C0003F800000E000007800003C00003E00001E00001F00001F0
3001F07801F0FC01F0FC01F0FC01E0F803E04003C02007801C0F000FFE0003F80014227EA019>
I<000300000300000700000F00000F0000170000370000270000470000C7000087000107000307
000207000407000C0700080700100700300700200700400700C00700FFFFFCFFFFFC0007000007
0000070000070000070000070000070000FFF800FFF816217FA019>I<1800401E03801FFF801F
FF001FFC0013F00010000010000010000010000010000010000010FC0013060014038018018010
01C00000E00000E00000F00000F00000F07000F0F800F0F800F0F800F0F800E0C001E04001C060
03C03003801C0F000FFE0003F00014227EA019>I<FFFFFFE0FFFFFFE0078003E0078000E00780
006007800020078000300780003007800010078000100780201007802010078020000780200007
8060000780E00007FFE00007FFE0000780E0000780600007802000078020000780200007802000
0780000007800000078000000780000007800000078000000780000007800000FFFE0000FFFE00
001C227EA121>70 D<01FC000707000C03801C01C03801C07800E07000E0F000E0FFFFE0F00000
F00000F00000F00000F000007000007800203800201C00400E008007030000FC0013157F9416>
101 D<00007801F08C071D1C0E0E1C1C07001C07003C07803C07803C07803C07801C07001C0700
0E0E000F1C0019F0001000001000001800001C00001FFF000FFFE00FFFF03800F860001840001C
C0000CC0000CC0000C6000186000183800700E01C001FE0016217F9519>103
D<1C003E003E003E001C00000000000000000000000000000000000E00FE00FE001E000E000E00
0E000E000E000E000E000E000E000E000E000E000E000E000E00FFC0FFC00A227FA10E>105
D<0E1E00FE6300FE87801E87800F03000F00000E00000E00000E00000E00000E00000E00000E00
000E00000E00000E00000E00000E00000E0000FFF000FFF00011157F9414>114
D<0E003800FE03F800FE03F8001E0078000E0038000E0038000E0038000E0038000E0038000E00
38000E0038000E0038000E0038000E0038000E0038000E0038000E0078000E0078000700BC0003
833F8000FC3F8019157F941C>117 D E end
%%EndProlog
%%BeginSetup
%%Feature: *Resolution 300
TeXDict begin 
%%EndSetup
%%Page: 1 1
bop 249 403 a @beginspecial 79 @llx -207 @lly 653 @urx -86
@ury 2784 @rwi @setspecial
%%BeginDocument: figbul5.eps
/$F2psDict 32 dict def $F2psDict begin $F2psDict /mtrx matrix put
 end /$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def /$F2psEnd {$F2psEnteredState restore end} def
$F2psBegin
0.900000 -0.900000 scale
1.000 .6 mul setlinewidth
newpath 476 108 moveto 449 146 lineto 458 151 lineto 434 186 lineto 425 180 lineto 401 213 lineto stroke
newpath 553 214 moveto 526 176 lineto 535 171 lineto 511 136 lineto 502 142 lineto 478 109 lineto stroke
newpath 401 213 moveto 423 205 lineto 430 218 lineto 455 210 lineto 449 197 lineto 476 188 lineto stroke
newpath 479 189 moveto 501 197 lineto 508 184 lineto 533 192 lineto 527 205 lineto 554 214 lineto stroke
newpath 724 123 moveto 702 115 lineto 695 128 lineto 670 120 lineto 676 107 lineto 649 98 lineto stroke
newpath 648 101 moveto 626 109 lineto 619 96 lineto 594 104 lineto 600 117 lineto 573 126 lineto stroke
newpath 725 125 moveto 698 163 lineto 707 168 lineto 683 203 lineto 674 197 lineto 650 230 lineto stroke
newpath 649 231 moveto 622 193 lineto 631 188 lineto 607 153 lineto 598 159 lineto 574 126 lineto stroke
newpath 175 122 moveto 175 134 lineto stroke
newpath 175 134 moveto 180 147 lineto stroke
newpath 179 146 moveto 174 146 lineto stroke
newpath 174 146 moveto 180 166 lineto stroke
newpath 180 166 moveto 174 165 lineto stroke
newpath 174 165 moveto 179 184 lineto stroke
newpath 179 184 moveto 174 184 lineto stroke
newpath 174 184 moveto 175 201 lineto stroke
newpath 162 112 moveto 162 123 lineto stroke
newpath 162 123 moveto 176 123 lineto stroke
newpath 98 121 moveto 111 121 lineto stroke
newpath 111 121 moveto 129 113 lineto stroke
newpath 129 113 moveto 130 121 lineto stroke
newpath 130 121 moveto 162 113 lineto stroke
newpath 96 199 moveto 96 187 lineto stroke
newpath 96 187 moveto 87 131 lineto stroke
newpath 87 131 moveto 96 131 lineto stroke
newpath 97 130 moveto 97 120 lineto stroke
newpath 175 201 moveto 165 201 lineto stroke
newpath 165 201 moveto 153 210 lineto stroke
newpath 154 210 moveto 155 201 lineto stroke
newpath 155 201 moveto 147 201 lineto stroke
newpath 147 201 moveto 139 209 lineto stroke
newpath 139 209 moveto 139 201 lineto stroke
newpath 138 202 moveto 130 202 lineto stroke
newpath 121 202 moveto 110 211 lineto stroke
newpath 110 211 moveto 110 202 lineto stroke
newpath 110 202 moveto 96 202 lineto stroke
newpath 133 202 moveto 122 210 lineto stroke
newpath 122 210 moveto 122 203 lineto stroke
newpath 313 119 moveto 303 119 lineto stroke
newpath 303 119 moveto 291 128 lineto stroke
newpath 292 128 moveto 293 119 lineto stroke
newpath 293 119 moveto 285 119 lineto stroke
newpath 285 119 moveto 277 127 lineto stroke
newpath 277 127 moveto 277 119 lineto stroke
newpath 276 120 moveto 268 120 lineto stroke
newpath 259 120 moveto 248 129 lineto stroke
newpath 248 129 moveto 248 120 lineto stroke
newpath 248 120 moveto 234 120 lineto stroke
newpath 271 120 moveto 260 128 lineto stroke
newpath 260 128 moveto 260 121 lineto stroke
newpath 312 198 moveto 312 186 lineto stroke
newpath 312 186 moveto 303 130 lineto stroke
newpath 303 130 moveto 312 130 lineto stroke
newpath 313 129 moveto 313 119 lineto stroke
newpath 265 200 moveto 297 192 lineto stroke
newpath 264 192 moveto 265 200 lineto stroke
newpath 246 200 moveto 264 192 lineto stroke
newpath 233 200 moveto 246 200 lineto stroke
newpath 297 202 moveto 311 202 lineto stroke
newpath 297 191 moveto 297 202 lineto stroke
newpath 231 122 moveto 231 134 lineto stroke
newpath 231 134 moveto 236 147 lineto stroke
newpath 235 146 moveto 230 146 lineto stroke
newpath 230 146 moveto 236 166 lineto stroke
newpath 236 166 moveto 230 165 lineto stroke
newpath 230 165 moveto 235 184 lineto stroke
newpath 235 184 moveto 230 184 lineto stroke
newpath 230 184 moveto 231 201 lineto stroke
$F2psEnd
%%EndDocument
 @endspecial 405 522 a Fa(Figure)15 b(1)596 b(Figure)15 b(2)771
1030 y @beginspecial 265 @llx -188 @lly 409 @urx -44 @ury 695
@rwi @setspecial
%%BeginDocument: figbul8.eps
/$F2psDict 32 dict def $F2psDict begin $F2psDict /mtrx matrix put
 end /$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def /$F2psEnd {$F2psEnteredState restore end} def
$F2psBegin
0.900000 -0.900000 scale
1.000 .6 mul setlinewidth
newpath 294 49 moveto 374 49 lineto 374 129 lineto 294 129 lineto 294 49 lineto stroke
newpath 374 49 moveto 454 49 lineto 454 129 lineto 374 129 lineto 374 49 lineto stroke
newpath 374 129 moveto 454 129 lineto 454 209 lineto 374 209 lineto 374 129 lineto stroke
newpath 294 129 moveto 374 129 lineto 374 209 lineto 294 209 lineto 294 129 lineto stroke
$F2psEnd
%%EndDocument
 @endspecial 854 1149 a(Figure)g(3)75 2379 y @beginspecial
50 @llx 50 @lly 554 @urx 554 @ury 3000 @rwi @setspecial
%%BeginDocument: gnuplot6.ps
/gnudict 40 dict def
gnudict begin
/Color false def
/gnulinewidth 5.000 def
/vshift -46 def
/dl {10 mul} def
/hpt 31.5 def
/vpt 31.5 def
/vpt2 vpt 2 mul def
/hpt2 hpt 2 mul def
/Lshow { currentpoint stroke moveto
  0 vshift rmoveto show } def
/Rshow { currentpoint stroke moveto
  dup stringwidth pop neg vshift rmoveto show } def
/Cshow { currentpoint stroke moveto
  dup stringwidth pop -2 div vshift rmoveto show } def
/DL { Color {setrgbcolor [] 0 setdash pop}
 {pop pop pop 0 setdash} ifelse } def
/BL { stroke gnulinewidth 2 mul setlinewidth } def
/AL { stroke gnulinewidth 2 div setlinewidth } def
/PL { stroke gnulinewidth setlinewidth } def
/LTb { BL [] 0 0 0 DL } def
/LTa { AL [1 dl 2 dl] 0 setdash 0 0 0 setrgbcolor } def
/LT0 { PL [] 0 1 0 DL } def
/LT1 { PL [4 dl 2 dl] 0 0 1 DL } def
/LT2 { PL [2 dl 3 dl] 1 0 0 DL } def
/LT3 { PL [1 dl 1.5 dl] 1 0 1 DL } def
/LT4 { PL [5 dl 2 dl 1 dl 2 dl] 0 1 1 DL } def
/LT5 { PL [4 dl 3 dl 1 dl 3 dl] 1 1 0 DL } def
/LT6 { PL [2 dl 2 dl 2 dl 4 dl] 0 0 0 DL } def
/LT7 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 1 0.3 0 DL } def
/LT8 { PL [2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 2 dl 4 dl] 0.5 0.5 0.5 DL } def
/M {moveto} def
/L {lineto} def
/P { stroke [] 0 setdash
  currentlinewidth 2 div sub moveto
  0 currentlinewidth rlineto  stroke } def
/D { stroke [] 0 setdash  2 copy  vpt add moveto
  hpt neg vpt neg rlineto  hpt vpt neg rlineto
  hpt vpt rlineto  hpt neg vpt rlineto  closepath  stroke
  P  } def
/A { stroke [] 0 setdash  vpt sub moveto  0 vpt2 rlineto
  currentpoint stroke moveto
  hpt neg vpt neg rmoveto  hpt2 0 rlineto stroke
  } def
/B { stroke [] 0 setdash  2 copy  exch hpt sub exch vpt add moveto
  0 vpt2 neg rlineto  hpt2 0 rlineto  0 vpt2 rlineto
  hpt2 neg 0 rlineto  closepath  stroke
  P  } def
/C { stroke [] 0 setdash  exch hpt sub exch vpt add moveto
  hpt2 vpt2 neg rlineto  currentpoint  stroke  moveto
  hpt2 neg 0 rmoveto  hpt2 vpt2 rlineto stroke  } def
/T { stroke [] 0 setdash  2 copy  vpt 1.12 mul add moveto
  hpt neg vpt -1.62 mul rlineto
  hpt 2 mul 0 rlineto
  hpt neg vpt 1.62 mul rlineto  closepath  stroke
  P  } def
/S { 2 copy A C} def
end
gnudict begin
gsave
50 50 translate
0.070 0.100 scale
0 setgray
/Courier findfont 140 scalefont setfont
newpath
LTa
LTb
LTb
LT0
LT0
3458 1986 M
3458 1986 L
4779 1986 L
4779 1986 M
4779 1986 L
6100 1986 L
6100 1986 M
6100 1986 L
6917 1986 L
3458 1986 M
3458 1986 L
4527 2533 L
4527 2533 M
4527 2533 L
5187 2871 L
4779 1986 M
4779 1986 L
5187 2871 L
4779 1986 M
4779 1986 L
5848 2533 L
6100 1986 M
6100 1986 L
5848 2533 L
6917 1986 M
6917 1986 L
5848 2533 L
5848 2533 M
5848 2533 L
5187 2871 L
4779 1986 M
4779 1986 L
4527 2533 L
1008 1411 M
1008 1411 L
1825 1411 L
1825 1411 M
1825 1411 L
3146 1411 L
1008 1411 M
1008 1411 L
1260 1958 L
1260 1958 M
1260 1958 L
2077 1958 L
2077 1958 M
2077 1958 L
3146 1411 L
3146 1411 M
3146 1411 L
2738 2296 L
2077 1958 M
2077 1958 L
2738 2296 L
1260 1958 M
1260 1958 L
1669 2843 L
2077 1958 M
2077 1958 L
1669 2843 L
2738 2296 M
2738 2296 L
1669 2843 L
2738 2296 M
2738 2296 L
2485 2843 L
1669 2843 M
1669 2843 L
2077 3728 L
2485 2843 M
2485 2843 L
2077 3728 L
1825 1411 M
1825 1411 L
2077 1958 L
1008 1411 M
1008 1411 L
2077 1958 L
1669 2843 M
1669 2843 L
2485 2843 L
stroke
grestore
end
showpage
%%EndDocument
 @endspecial 854 2049 a(Figure)g(4)249 2506 y @beginspecial
24 @llx -528 @lly 703 @urx -383 @ury 2784 @rwi @setspecial
%%BeginDocument: figbul7.eps
/$F2psDict 32 dict def $F2psDict begin $F2psDict /mtrx matrix put
 end /$F2psBegin {$F2psDict begin /$F2psEnteredState save def} def /$F2psEnd {$F2psEnteredState restore end} def
$F2psBegin
0.900000 -0.900000 scale
1.000 .6 mul setlinewidth
newpath 701 587 moveto 541 507 lineto 381 427 lineto 541 427 lineto stroke
newpath 541 507 moveto 541 427 lineto 701 427 lineto 701 427 lineto 781 427 lineto 701 587 lineto 701 587 lineto 701 587 lineto 701 587 lineto 701 507 lineto closepath  stroke
newpath 701 427 moveto 541 507 lineto stroke
newpath 701 507 moveto 701 507 lineto 701 427 lineto closepath  stroke
newpath 541 507 moveto 701 507 lineto stroke
newpath 701 587 moveto 541 507 lineto 381 427 lineto 541 427 lineto stroke
newpath 541 507 moveto 541 427 lineto 701 427 lineto 701 427 lineto 781 427 lineto 701 587 lineto 701 587 lineto 701 587 lineto 701 587 lineto 701 507 lineto closepath  stroke
newpath 701 427 moveto 541 507 lineto stroke
newpath 701 507 moveto 701 507 lineto 701 427 lineto closepath  stroke
newpath 541 507 moveto 701 507 lineto stroke
newpath 186 506 moveto 346 506 lineto stroke
newpath 346 506 moveto 346 506 lineto 346 586 lineto closepath  stroke
newpath 346 586 moveto 186 506 lineto stroke
newpath 186 506 moveto 186 586 lineto 346 586 lineto 346 586 lineto 426 586 lineto 346 426 lineto 346 426 lineto 346 426 lineto 346 426 lineto 346 506 lineto closepath  stroke
newpath 346 426 moveto 186 506 lineto 26 586 lineto 186 586 lineto stroke
newpath 186 506 moveto 346 506 lineto stroke
newpath 346 506 moveto 346 506 lineto 346 586 lineto closepath  stroke
newpath 346 586 moveto 186 506 lineto stroke
newpath 186 506 moveto 186 586 lineto 346 586 lineto 346 586 lineto 426 586 lineto 346 426 lineto 346 426 lineto 346 426 lineto 346 426 lineto 346 506 lineto closepath  stroke
newpath 346 426 moveto 186 506 lineto 26 586 lineto 186 586 lineto stroke
$F2psEnd
%%EndDocument
 @endspecial 854 2625 a(Figure)g(5)p eop
%%Trailer
end
userdict /end-hook known{end-hook}if
%%EOF

