c c AUTOP.FOR (Autocorrelation for the Kundic (main) c event in 1995) QSO 0957+561A,B - April/May 2001 c c implicit double precision (a-h),(o-z) dimension ta(13),a(13),ea(13),b(13),eb(13) dimension am(13),eam(13) dimension siga(13),udcf(13,13),tlag(13,13) open(unit=1,file='autop15.dat',status='new') open(unit=2,file='95p.dat',status='old') read(2,*)(ta(i),a(i),ea(i),b(i),eb(i),i=1,13) do 10 i=1,13 if(i.eq.1)am(i)=a(i) if(i.gt.1)am(i)=a(i)+am(i-1) if(i.eq.1)eam(i)=ea(i) if(i.gt.1)eam(i)=ea(i)+eam(i-1) 10 continue amedio=am(13)/13.D0 errora2=(eam(13)/13.D0)**2.D0 do 30 i=1,13 if(i.eq.1)siga(i)=(a(i)-amedio)**2.D0 if(i.gt.1)siga(i)=(a(i)-amedio)**2.D0+siga(i-1) 30 continue sigmaa2=siga(13)/12.D0 raiz=sigmaa2-errora2 do 50 i=1,13 do 60 j=1,13 udcf(i,j)=(a(i)-amedio)*(a(j)-amedio)/raiz tlag(i,j)=ta(j)-ta(i) 60 continue 50 continue c *********************************************************** c BINS WITH SIZE OF 2*alfa c *********************************************************** alfa=7.5D0 c *********************************************************** do 70 k=1,51 tau=-26.D0+k bdcf=0.D0 bdat=0.D0 do 80 i=1,13 do 90 j=1,13 if(tlag(i,j).ge.(tau-alfa).and.tlag(i,j).lt.(tau+alfa)) *bdcf=bdcf+udcf(i,j) if(tlag(i,j).ge.(tau-alfa).and.tlag(i,j).lt.(tau+alfa)) *bdat=bdat+1.D0 90 continue 80 continue if(bdat.ge.1.D0)write(1,200)tau,bdcf/bdat 70 continue stop 200 format(E10.3,4x,E10.3) end