program fracvar
c
c     determine fractional variance of series2 explained by series1
c
      parameter (nmax=4000)
      real series1(nmax),series2(nmax)
      open (unit=4,file='fracvar.in',status='old')
      do i=1,nmax
         read (4,*,end=99) idum,series1(i),series2(i)
      end do
99    continue
      close (unit=4)
      nx = i-1
c
      sum1 = 0.0
      sum2 = 0.0
      do i=1,nx
         sum1 = sum1+ (series2(i)-series1(i))**2
         sum2 = sum2 + series2(i)**2
      end do
      beta = 1.0-sum1/sum2
      write (6,*) 'explained variance: ',beta
      stop
      end