1 brash 1.1 vec/create phi(180)
 2           vec/create phi2(180)
 3           vec/read phi phi.dat
 4           vec/read phi2 phi2.dat
 5           1d 5001 ' ' 180 -50 450
 6           1d 5002 ' ' 180 -50 450
 7           hist/put_vec/contents 5001 phi
 8           hist/put_vec/contents 5002 phi2
 9           1d 5003 ' ' 180 -50 450
10           add 5001 5001 5003 .0000115757 0.0
11           hist/set/max 5003 1.05
12           hist/set/min 5003 0.95
13           1d 5004 ' ' 180 -50 450
14           add 5002 5002 5004 .0000183348 0.0
15           hist/set/min 5004 0.95
16           hist/set/max 5004 1.05
17           opt nstat
18           opt fit
19           fort/file 66 phiplots_normalized.ps
20           meta 66
21           zone 1 2
22 brash 1.1 hist/fit 5003(0.:360.) p0
23           hist/fit 5004(0.:360.) p0
24           vec/create x(45)
25           vec/create y(45)
26           vec/create ex(45)
27           vec/create ey(45)
28           1d 5005 ' ' 45 -50 450
29           1d 5006 ' ' 45 -50 450
30           rebin 5003 x y ex ey 45 1 180
31           hist/put_vec/contents 5005 y
32           rebin 5004 x y ex ey 45 1 180
33           hist/put_vec/contents 5006 y
34           hist/set/min 5005 0.95
35           hist/set/max 5005 1.05
36           hist/set/min 5006 0.95
37           hist/set/max 5006 1.05
38           hist/fit 5005(15.:345.) p0
39           hist/fit 5006(15.:345.) p0
40           close 66
41           shell 'ps2pdf phiplots_normalized.ps'