! device x11 -g 800x800-100+10 -bg RoyalBlue -fg white
 device postencap Pkevolution32.eps
macro prf 5 { lines $1 $2
              read k $4
		        read pk $5
	           ltype 0
	           ltype $3
              connect (LG(k)) (LG(pk)) if (pk > 0)
              ltype 0
	      }
location 8000 30000 5000 25000
lweight 2
expand 1.1
	      limits -1 1.55 -2.30 4.3
	      ticksize -1 10 -1 10
	      box
 	      xlabel k(h/Mpc)
		   ylabel P(k) (Mpc/h)^{3}
data Sp32x128.init.dat
prf 61 171 0 2 4
data Sp32x128.01225.dat
prf 61 171 0 2 4
data Sp32x128.a100.dat
prf 61 171 0 2 4
data model.dat
prf 17 126 2 1 2
prf 129 238 2 1 2
data model.dat.a0.1225
prf 18 127 2 1 2

expand 0.9
relocate 0.3 3.5 label \Lambda CDM model: h=0.7 \Omega_0=0.3
relocate 0.3 3.2 label PM run with 32^3 particles
relocate 0.7 2.9 label 128^3 mesh
relocate 0.7 2.6 label 20h^{-1}Mpc box
		ltype 1
relocate 0.701 -3  draw 0.701 2.2
		ltype 0
		relocate 0.75 1.9 label Nyquist Frequency
relocate -0.6 0.35 label z=15
relocate -0.6 1.93 label z=7.1
relocate -0.6 3.4 label z=0

 hardcopy
 end