c c ACCP5FIG.FOR c real x1(51),y1(51) real x2(51),y2(51) real x3(51) call pgbegin(0,'?',1,1) call pgslw(2) call pgscf(2) open(unit=1,file='crosp5.dat',status='old') open(unit=2,file='autop5.dat',status='old') read(1,*)(x1(i),y1(i),i=1,51) read(2,*)(x2(i),y2(i),i=1,51) do 10 i=1,51 x3(i)=x2(i)+417. 10 continue call pgenv(394.,446.,-1.5,1.5,0,1) call pglabel('Lag (days)', *'DCC vs. DAC(+ 417 days)','Main events (Bin = 5 days)') call pgpoint(51,x1,y1,17) call pgpoint(51,x3,y2,23) call pgend end