ve/del *
 ve/read e,numu1,numu2,numu3    check/2900-numu-0.vec
 ve/read e,nue91,nue92,nue93    check/2900-nue-90.vec
 ve/read e,numu91,numu92,numu93 check/2900-numu-90.vec
 sigma e=e*1.0e-9
 
 1d 1 'delta=0' 1000 0.1 10
 1d 2 'numu->nue delta=90' 1000 0.1 10
 1d 3 'nue->numu detla=90' 1000 0.1 10
 
 put/cont 1 numu1
 put/cont 2 numu91
 put/cont 3 nue92
 
 opt logx
 title 'Effects of CP phase'
 null 1 10 0 0.2
 set htyp 
 hi/pl 1 s
 hi/pl 2 s
 hi/pl 3 s
 atitle 'E(GeV)' 'Probability'