Add worked example

1 files changed, 1049 insertions, 0 deletions

diff --git a/figures/error_diffusion.eps b/figures/error_diffusion.eps
new file mode 100644
index 0000000..1405259
--- /dev/null
+++ b/figures/error_diffusion.eps
@@ -0,0 +1,1049 @@
+%!PS-Adobe-2.0 EPSF-2.0
+%%Title: error_diffusion.eps
+%%Creator: gnuplot 4.4 patchlevel 3
+%%CreationDate: Tue Mar 24 10:35:40 2015
+%%DocumentFonts: (atend)
+%%BoundingBox: 50 50 410 302
+%%EndComments
+%%BeginProlog
+/gnudict 256 dict def
+gnudict begin
+%
+% The following true/false flags may be edited by hand if desired.
+% The unit line width and grayscale image gamma correction may also be changed.
+%
+/Color true def
+/Blacktext false def
+/Solid false def
+/Dashlength 1 def
+/Landscape false def
+/Level1 false def
+/Rounded false def
+/ClipToBoundingBox false def
+/TransparentPatterns false def
+/gnulinewidth 5.000 def
+/userlinewidth gnulinewidth def
+/Gamma 1.0 def
+%
+/vshift -66 def
+/dl1 {
+  10.0 Dashlength mul mul
+  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
+} def
+/dl2 {
+  10.0 Dashlength mul mul
+  Rounded { currentlinewidth 0.75 mul add } if
+} def
+/hpt_ 31.5 def
+/vpt_ 31.5 def
+/hpt hpt_ def
+/vpt vpt_ def
+Level1 {} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (error_diffusion.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 4.4 patchlevel 3)
+  /Author (stoltz)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Tue Mar 24 10:35:40 2015)
+  /DOCINFO pdfmark
+end
+} ifelse
+/doclip {
+  ClipToBoundingBox {
+    newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
+    clip
+  } if
+} def
+%
+% Gnuplot Prolog Version 4.4 (August 2010)
+%
+%/SuppressPDFMark true def
+%
+/M {moveto} bind def
+/L {lineto} bind def
+/R {rmoveto} bind def
+/V {rlineto} bind def
+/N {newpath moveto} bind def
+/Z {closepath} bind def
+/C {setrgbcolor} bind def
+/f {rlineto fill} bind def
+/g {setgray} bind def
+/Gshow {show} def   % May be redefined later in the file to support UTF-8
+/vpt2 vpt 2 mul def
+/hpt2 hpt 2 mul def
+/Lshow {currentpoint stroke M 0 vshift R 
+	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
+/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
+	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
+/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
+	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
+/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
+  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
+/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
+ {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
+/BL {stroke userlinewidth 2 mul setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/AL {stroke userlinewidth 2 div setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+/UL {dup gnulinewidth mul /userlinewidth exch def
+	dup 1 lt {pop 1} if 10 mul /udl exch def} def
+/PL {stroke userlinewidth setlinewidth
+	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
+% Default Line colors
+/LCw {1 1 1} def
+/LCb {0 0 0} def
+/LCa {0 0 0} def
+/LC0 {1 0 0} def
+/LC1 {0 1 0} def
+/LC2 {0 0 1} def
+/LC3 {1 0 1} def
+/LC4 {0 1 1} def
+/LC5 {1 1 0} def
+/LC6 {0 0 0} def
+/LC7 {1 0.3 0} def
+/LC8 {0.5 0.5 0.5} def
+% Default Line Types
+/LTw {PL [] 1 setgray} def
+/LTb {BL [] LCb DL} def
+/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
+/LT0 {PL [] LC0 DL} def
+/LT1 {PL [4 dl1 2 dl2] LC1 DL} def
+/LT2 {PL [2 dl1 3 dl2] LC2 DL} def
+/LT3 {PL [1 dl1 1.5 dl2] LC3 DL} def
+/LT4 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
+/LT5 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC5 DL} def
+/LT6 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC6 DL} def
+/LT7 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC7 DL} def
+/LT8 {PL [2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 4 dl2] LC8 DL} def
+/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
+/Dia {stroke [] 0 setdash 2 copy vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke
+  Pnt} def
+/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
+  currentpoint stroke M
+  hpt neg vpt neg R hpt2 0 V stroke
+ } def
+/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke
+  Pnt} def
+/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
+  hpt2 vpt2 neg V currentpoint stroke M
+  hpt2 neg 0 R hpt2 vpt2 V stroke} def
+/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke
+  Pnt} def
+/Star {2 copy Pls Crs} def
+/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath fill} def
+/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath fill} def
+/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke
+  Pnt} def
+/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath fill} def
+/DiaF {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath fill} def
+/Pent {stroke [] 0 setdash 2 copy gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore Pnt} def
+/PentF {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath fill grestore} def
+/Circle {stroke [] 0 setdash 2 copy
+  hpt 0 360 arc stroke Pnt} def
+/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
+/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
+/C1 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C2 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C3 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C4 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C5 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc
+	2 copy moveto
+	2 copy vpt 180 270 arc closepath fill
+	vpt 0 360 arc} bind def
+/C6 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C7 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 270 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C8 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C9 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 270 450 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
+	2 copy moveto
+	2 copy vpt 90 180 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C11 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 180 arc closepath fill
+	2 copy moveto
+	2 copy vpt 270 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C12 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C13 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 0 90 arc closepath fill
+	2 copy moveto
+	2 copy vpt 180 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/C14 {BL [] 0 setdash 2 copy moveto
+	2 copy vpt 90 360 arc closepath fill
+	vpt 0 360 arc} bind def
+/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
+	vpt 0 360 arc closepath} bind def
+/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
+	neg 0 rlineto closepath} bind def
+/Square {dup Rec} bind def
+/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
+/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
+/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
+/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
+/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
+	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
+/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
+/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
+/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
+	Bsquare} bind def
+/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
+	Bsquare} bind def
+/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
+/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy vpt Square fill Bsquare} bind def
+/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
+	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
+/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
+/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
+/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
+/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
+/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
+/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
+/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
+/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
+/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
+/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
+/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
+/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
+/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
+/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
+/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
+/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
+/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
+/DiaE {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V closepath stroke} def
+/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V closepath stroke} def
+/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V closepath stroke} def
+/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V closepath stroke} def
+/PentE {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  closepath stroke grestore} def
+/CircE {stroke [] 0 setdash 
+  hpt 0 360 arc stroke} def
+/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
+/DiaW {stroke [] 0 setdash vpt add M
+  hpt neg vpt neg V hpt vpt neg V
+  hpt vpt V hpt neg vpt V Opaque stroke} def
+/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
+  0 vpt2 neg V hpt2 0 V 0 vpt2 V
+  hpt2 neg 0 V Opaque stroke} def
+/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
+  hpt neg vpt -1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt 1.62 mul V Opaque stroke} def
+/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
+  hpt neg vpt 1.62 mul V
+  hpt 2 mul 0 V
+  hpt neg vpt -1.62 mul V Opaque stroke} def
+/PentW {stroke [] 0 setdash gsave
+  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
+  Opaque stroke grestore} def
+/CircW {stroke [] 0 setdash 
+  hpt 0 360 arc Opaque stroke} def
+/BoxFill {gsave Rec 1 setgray fill grestore} def
+/Density {
+  /Fillden exch def
+  currentrgbcolor
+  /ColB exch def /ColG exch def /ColR exch def
+  /ColR ColR Fillden mul Fillden sub 1 add def
+  /ColG ColG Fillden mul Fillden sub 1 add def
+  /ColB ColB Fillden mul Fillden sub 1 add def
+  ColR ColG ColB setrgbcolor} def
+/BoxColFill {gsave Rec PolyFill} def
+/PolyFill {gsave Density fill grestore grestore} def
+/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
+%
+% PostScript Level 1 Pattern Fill routine for rectangles
+% Usage: x y w h s a XX PatternFill
+%	x,y = lower left corner of box to be filled
+%	w,h = width and height of box
+%	  a = angle in degrees between lines and x-axis
+%	 XX = 0/1 for no/yes cross-hatch
+%
+/PatternFill {gsave /PFa [ 9 2 roll ] def
+  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
+  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
+  gsave 1 setgray fill grestore clip
+  currentlinewidth 0.5 mul setlinewidth
+  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
+  0 0 M PFa 5 get rotate PFs -2 div dup translate
+  0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 M 0 PFs V} for
+  0 PFa 6 get ne {
+	0 1 PFs PFa 4 get div 1 add floor cvi
+	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
+ } if
+  stroke grestore} def
+%
+/languagelevel where
+ {pop languagelevel} {1} ifelse
+ 2 lt
+	{/InterpretLevel1 true def}
+	{/InterpretLevel1 Level1 def}
+ ifelse
+%
+% PostScript level 2 pattern fill definitions
+%
+/Level2PatternFill {
+/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
+	bind def
+/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
+>> matrix makepattern
+/Pat1 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
+	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
+>> matrix makepattern
+/Pat2 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
+	8 8 L 8 0 L 0 0 L fill}
+>> matrix makepattern
+/Pat3 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
+	0 12 M 12 0 L stroke}
+>> matrix makepattern
+/Pat4 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
+	0 -4 M 12 8 L stroke}
+>> matrix makepattern
+/Pat5 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
+	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
+>> matrix makepattern
+/Pat6 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
+	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
+>> matrix makepattern
+/Pat7 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
+	12 0 M -4 8 L 12 4 M 0 10 L stroke}
+>> matrix makepattern
+/Pat8 exch def
+<< Tile8x8
+ /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
+	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
+>> matrix makepattern
+/Pat9 exch def
+/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
+/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
+/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
+/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
+/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
+/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
+/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
+} def
+%
+%
+%End of PostScript Level 2 code
+%
+/PatternBgnd {
+  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
+} def
+%
+% Substitute for Level 2 pattern fill codes with
+% grayscale if Level 2 support is not selected.
+%
+/Level1PatternFill {
+/Pattern1 {0.250 Density} bind def
+/Pattern2 {0.500 Density} bind def
+/Pattern3 {0.750 Density} bind def
+/Pattern4 {0.125 Density} bind def
+/Pattern5 {0.375 Density} bind def
+/Pattern6 {0.625 Density} bind def
+/Pattern7 {0.875 Density} bind def
+} def
+%
+% Now test for support of Level 2 code
+%
+Level1 {Level1PatternFill} {Level2PatternFill} ifelse
+%
+/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
+dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
+currentdict end definefont pop
+/MFshow {
+   { dup 5 get 3 ge
+     { 5 get 3 eq {gsave} {grestore} ifelse }
+     {dup dup 0 get findfont exch 1 get scalefont setfont
+     [ currentpoint ] exch dup 2 get 0 exch R dup 5 get 2 ne {dup dup 6
+     get exch 4 get {Gshow} {stringwidth pop 0 R} ifelse }if dup 5 get 0 eq
+     {dup 3 get {2 get neg 0 exch R pop} {pop aload pop M} ifelse} {dup 5
+     get 1 eq {dup 2 get exch dup 3 get exch 6 get stringwidth pop -2 div
+     dup 0 R} {dup 6 get stringwidth pop -2 div 0 R 6 get
+     show 2 index {aload pop M neg 3 -1 roll neg R pop pop} {pop pop pop
+     pop aload pop M} ifelse }ifelse }ifelse }
+     ifelse }
+   forall} def
+/Gswidth {dup type /stringtype eq {stringwidth} {pop (n) stringwidth} ifelse} def
+/MFwidth {0 exch { dup 5 get 3 ge { 5 get 3 eq { 0 } { pop } ifelse }
+ {dup 3 get{dup dup 0 get findfont exch 1 get scalefont setfont
+     6 get Gswidth pop add} {pop} ifelse} ifelse} forall} def
+/MLshow { currentpoint stroke M
+  0 exch R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MRshow { currentpoint stroke M
+  exch dup MFwidth neg 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/MCshow { currentpoint stroke M
+  exch dup MFwidth -2 div 3 -1 roll R
+  Blacktext {gsave 0 setgray MFshow grestore} {MFshow} ifelse } bind def
+/XYsave    { [( ) 1 2 true false 3 ()] } bind def
+/XYrestore { [( ) 1 2 true false 4 ()] } bind def
+end
+%%EndProlog
+gnudict begin
+gsave
+doclip
+50 50 translate
+0.050 0.050 scale
+0 setgray
+newpath
+(Helvetica) findfont 200 scalefont setfont
+1.000 UL
+LTb
+980 640 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 640 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.5)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 1280 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 1280 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.52)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 1920 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 1920 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.54)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 2560 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 2560 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.56)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 3199 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 3199 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.58)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 3839 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 3839 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.6)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 4479 M
+63 0 V
+5796 0 R
+-63 0 V
+stroke
+860 4479 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.62)]
+] -66.7 MRshow
+1.000 UL
+LTb
+980 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+980 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0)]
+] -66.7 MCshow
+1.000 UL
+LTb
+2096 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+2096 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.002)]
+] -66.7 MCshow
+1.000 UL
+LTb
+3212 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+3212 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.004)]
+] -66.7 MCshow
+1.000 UL
+LTb
+4328 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+4328 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.006)]
+] -66.7 MCshow
+1.000 UL
+LTb
+5444 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+5444 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.008)]
+] -66.7 MCshow
+1.000 UL
+LTb
+6560 640 M
+0 63 V
+0 4096 R
+0 -63 V
+stroke
+6560 440 M
+[ [(Helvetica) 200.0 0.0 true true 0 ( 0.01)]
+] -66.7 MCshow
+1.000 UL
+LTb
+1.000 UL
+LTb
+980 4799 N
+980 640 L
+5859 0 V
+0 4159 V
+-5859 0 V
+Z stroke
+LCb setrgbcolor
+160 2719 M
+currentpoint gsave translate 90 rotate 0 0 moveto
+[ [(Helvetica) 300.0 0.0 true true 0 (Diffusion)]
+] -100.0 MCshow
+grestore
+LTb
+LCb setrgbcolor
+4029 140 M
+[ [(Helvetica) 300.0 0.0 true true 0 (Time step)]
+] -100.0 MCshow
+LTb
+1.000 UP
+1.000 UL
+LTb
+% Begin plot #1
+1.500 UP
+1.000 UL
+LT0
+LC0 setrgbcolor
+LCb setrgbcolor
+2708 2460 M
+[ [(Helvetica) 150.0 0.0 true true 0 (GK/Barker/HMC)]
+] -50.0 MRshow
+LT0
+LC0 setrgbcolor
+2828 2460 M
+543 0 V
+1036 4600 M
+56 0 V
+55 6 V
+56 -7 V
+56 3 V
+56 3 V
+56 -4 V
+55 0 V
+56 1 V
+56 4 V
+558 5 V
+558 0 V
+558 5 V
+558 8 V
+558 9 V
+558 3 V
+558 5 V
+558 3 V
+558 3 V
+1036 4600 CircleF
+1092 4600 CircleF
+1147 4606 CircleF
+1203 4599 CircleF
+1259 4602 CircleF
+1315 4605 CircleF
+1371 4601 CircleF
+1426 4601 CircleF
+1482 4602 CircleF
+1538 4606 CircleF
+2096 4611 CircleF
+2654 4611 CircleF
+3212 4616 CircleF
+3770 4624 CircleF
+4328 4633 CircleF
+4886 4636 CircleF
+5444 4641 CircleF
+6002 4644 CircleF
+6560 4647 CircleF
+3099 2460 CircleF
+% End plot #1
+% Begin plot #2
+1.500 UP
+1.000 UL
+LT0
+LC0 setrgbcolor
+LCb setrgbcolor
+2708 2260 M
+[ [(Helvetica) 150.0 0.0 true true 0 (Einstein/Barker/HMC)]
+] -50.0 MRshow
+LT0
+LC0 setrgbcolor
+2828 2260 M
+543 0 V
+1036 4619 M
+56 6 V
+55 -2 V
+56 -16 V
+56 0 V
+56 -1 V
+56 4 V
+55 7 V
+56 -4 V
+56 -4 V
+558 -29 V
+558 -7 V
+558 -19 V
+558 -33 V
+558 -44 V
+558 -54 V
+558 -61 V
+558 -67 V
+558 -84 V
+1036 4619 BoxF
+1092 4625 BoxF
+1147 4623 BoxF
+1203 4607 BoxF
+1259 4607 BoxF
+1315 4606 BoxF
+1371 4610 BoxF
+1426 4617 BoxF
+1482 4613 BoxF
+1538 4609 BoxF
+2096 4580 BoxF
+2654 4573 BoxF
+3212 4554 BoxF
+3770 4521 BoxF
+4328 4477 BoxF
+4886 4423 BoxF
+5444 4362 BoxF
+6002 4295 BoxF
+6560 4211 BoxF
+3099 2260 BoxF
+% End plot #2
+% Begin plot #3
+1.500 UP
+1.000 UL
+LT0
+LC2 setrgbcolor
+LCb setrgbcolor
+2708 2060 M
+[ [(Helvetica) 150.0 0.0 true true 0 (GK/Barker/midpoint)]
+] -50.0 MRshow
+LT0
+LC2 setrgbcolor
+2828 2060 M
+543 0 V
+1036 4606 M
+56 0 V
+55 -3 V
+56 -2 V
+56 1 V
+56 3 V
+56 -2 V
+55 0 V
+56 0 V
+56 0 V
+558 -3 V
+558 -4 V
+558 -5 V
+558 -6 V
+558 -16 V
+558 -10 V
+558 -19 V
+558 -19 V
+558 -22 V
+1036 4606 CircleF
+1092 4606 CircleF
+1147 4603 CircleF
+1203 4601 CircleF
+1259 4602 CircleF
+1315 4605 CircleF
+1371 4603 CircleF
+1426 4603 CircleF
+1482 4603 CircleF
+1538 4603 CircleF
+2096 4600 CircleF
+2654 4596 CircleF
+3212 4591 CircleF
+3770 4585 CircleF
+4328 4569 CircleF
+4886 4559 CircleF
+5444 4540 CircleF
+6002 4521 CircleF
+6560 4499 CircleF
+3099 2060 CircleF
+% End plot #3
+% Begin plot #4
+1.500 UP
+1.000 UL
+LT0
+LC2 setrgbcolor
+LCb setrgbcolor
+2708 1860 M
+[ [(Helvetica) 150.0 0.0 true true 0 (Einstein/Barker/midpoint)]
+] -50.0 MRshow
+LT0
+LC2 setrgbcolor
+2828 1860 M
+543 0 V
+1036 4615 M
+56 -8 V
+55 16 V
+56 -11 V
+56 -6 V
+56 -1 V
+56 -3 V
+55 -2 V
+56 1 V
+56 1 V
+558 -4 V
+558 -35 V
+558 -27 V
+558 -42 V
+558 -58 V
+558 -70 V
+558 -87 V
+558 -106 V
+558 -114 V
+1036 4615 BoxF
+1092 4607 BoxF
+1147 4623 BoxF
+1203 4612 BoxF
+1259 4606 BoxF
+1315 4605 BoxF
+1371 4602 BoxF
+1426 4600 BoxF
+1482 4601 BoxF
+1538 4602 BoxF
+2096 4598 BoxF
+2654 4563 BoxF
+3212 4536 BoxF
+3770 4494 BoxF
+4328 4436 BoxF
+4886 4366 BoxF
+5444 4279 BoxF
+6002 4173 BoxF
+6560 4059 BoxF
+3099 1860 BoxF
+% End plot #4
+% Begin plot #5
+1.500 UP
+1.000 UL
+LT0
+LC8 setrgbcolor
+LCb setrgbcolor
+2708 1660 M
+[ [(Helvetica) 150.0 0.0 true true 0 (Einstein/Metropolis/HMC)]
+] -50.0 MRshow
+LT0
+LC8 setrgbcolor
+2828 1660 M
+543 0 V
+1036 4594 M
+56 -3 V
+55 3 V
+56 -4 V
+56 -5 V
+56 -4 V
+56 -6 V
+55 -10 V
+56 -3 V
+56 -3 V
+558 -63 V
+558 -80 V
+558 -98 V
+558 -133 V
+558 -142 V
+558 -151 V
+558 -164 V
+558 -176 V
+558 -178 V
+1036 4594 BoxF
+1092 4591 BoxF
+1147 4594 BoxF
+1203 4590 BoxF
+1259 4585 BoxF
+1315 4581 BoxF
+1371 4575 BoxF
+1426 4565 BoxF
+1482 4562 BoxF
+1538 4559 BoxF
+2096 4496 BoxF
+2654 4416 BoxF
+3212 4318 BoxF
+3770 4185 BoxF
+4328 4043 BoxF
+4886 3892 BoxF
+5444 3728 BoxF
+6002 3552 BoxF
+6560 3374 BoxF
+3099 1660 BoxF
+% End plot #5
+% Begin plot #6
+1.500 UP
+1.000 UL
+LT0
+LC8 setrgbcolor
+LCb setrgbcolor
+2708 1460 M
+[ [(Helvetica) 150.0 0.0 true true 0 (GK/Metropolis/HMC)]
+] -50.0 MRshow
+LT0
+LC8 setrgbcolor
+2828 1460 M
+543 0 V
+1036 4606 M
+56 -6 V
+55 4 V
+56 -7 V
+56 -4 V
+56 -2 V
+56 -8 V
+55 1 V
+56 -8 V
+56 1 V
+558 -46 V
+558 -66 V
+558 -67 V
+558 -82 V
+558 -81 V
+558 -93 V
+558 -101 V
+558 -107 V
+558 -111 V
+1036 4606 CircleF
+1092 4600 CircleF
+1147 4604 CircleF
+1203 4597 CircleF
+1259 4593 CircleF
+1315 4591 CircleF
+1371 4583 CircleF
+1426 4584 CircleF
+1482 4576 CircleF
+1538 4577 CircleF
+2096 4531 CircleF
+2654 4465 CircleF
+3212 4398 CircleF
+3770 4316 CircleF
+4328 4235 CircleF
+4886 4142 CircleF
+5444 4041 CircleF
+6002 3934 CircleF
+6560 3823 CircleF
+3099 1460 CircleF
+% End plot #6
+% Begin plot #7
+1.500 UP
+1.000 UL
+LT0
+LC6 setrgbcolor
+LCb setrgbcolor
+2708 1260 M
+[ [(Helvetica) 150.0 0.0 true true 0 (GK/reference)]
+] -50.0 MRshow
+LT0
+LC6 setrgbcolor
+2828 1260 M
+543 0 V
+1036 4601 M
+56 -12 V
+55 -10 V
+56 -10 V
+56 -6 V
+56 -11 V
+111 -34 V
+112 -26 V
+558 -157 V
+558 -170 V
+558 -191 V
+4328 3533 L
+5444 3042 L
+6560 2492 L
+1036 4601 CircleF
+1092 4589 CircleF
+1147 4579 CircleF
+1203 4569 CircleF
+1259 4563 CircleF
+1315 4552 CircleF
+1426 4518 CircleF
+1538 4492 CircleF
+2096 4335 CircleF
+2654 4165 CircleF
+3212 3974 CircleF
+4328 3533 CircleF
+5444 3042 CircleF
+6560 2492 CircleF
+3099 1260 CircleF
+% End plot #7
+% Begin plot #8
+1.500 UP
+1.000 UL
+LT0
+LC6 setrgbcolor
+LCb setrgbcolor
+2708 1060 M
+[ [(Helvetica) 150.0 0.0 true true 0 (Einstein/reference)]
+] -50.0 MRshow
+LT0
+LC6 setrgbcolor
+2828 1060 M
+543 0 V
+1036 4576 M
+56 -57 V
+55 -22 V
+56 -33 V
+56 -44 V
+56 -39 V
+111 -71 V
+112 -85 V
+558 -433 V
+558 -423 V
+558 -423 V
+4328 2166 L
+5444 1413 L
+6560 693 L
+1036 4576 BoxF
+1092 4519 BoxF
+1147 4497 BoxF
+1203 4464 BoxF
+1259 4420 BoxF
+1315 4381 BoxF
+1426 4310 BoxF
+1538 4225 BoxF
+2096 3792 BoxF
+2654 3369 BoxF
+3212 2946 BoxF
+4328 2166 BoxF
+5444 1413 BoxF
+6560 693 BoxF
+3099 1060 BoxF
+% End plot #8
+1.000 UL
+LTb
+980 4799 N
+980 640 L
+5859 0 V
+0 4159 V
+-5859 0 V
+Z stroke
+1.000 UP
+1.000 UL
+LTb
+stroke
+grestore
+end
+showpage
+%%Trailer
+%%DocumentFonts: Helvetica