dimension b(100,3),avg(3)
	open(30, file='repeat')
	do i=1, 3
	  avg(i)=0.
	end do

	do i=1, 100
	  read(30,*,end=5) (b(i,j), j=1, 3)
          k=i
	  do j=1, 3
 	    avg(j)=avg(j)+b(i,j)
	  end do
	end do

5       continue
	do j=1, 3
 	  avg(j)=avg(j)/k
	end do

        write(*,100)
100	format('Data is generated as yt= 3 + 2*x1 + x2'/
     *         'ie.  beta{} = {3   2   1}'/
     *         'ye=yt + eps ;  where eps = RN(mu=0, sigma=0.2)'/
     *         'Several data sets were generated '/
     *         'Each data set gives one set of b'/
     *         'Here are the sets of b and also its average'//
     *         'S. No.         b1              b2            b3'/)
        write(*,101) (i, (b(i,j), j=1,3), i=1, k)
101     format(i3, 3f16.5)
        write(*,*)'Average of above'
        write(*,102) (avg(i), i=1, 3)
102     format(3x, 3f16.5)
        STOP
        END