!----------------!
 program averages
!----------------!
 implicit none

 integer :: bins
 real(8) :: a(4),av(2,4)

 bins=0
 av=0.d0
 open (10,file='res.dat',status='old')
 do 
   read(10,*,end=10)a(:)
   bins=bins+1
   av(1,:)=av(1,:)+a(:)
   av(2,:)=av(2,:)+a(:)**2
 end do
 10 close(10)
 write(*,*)bins
 av=av/bins
 av(2,:)=sqrt((av(2,:)-av(1,:)**2)/(bins-1))
 
 open (10,file='e.dat')
 write(10,'(4f15.8)')av(:,1:2)
 close(10)

 open (10,file='m.dat')
 write(10,'(4f15.8)')av(:,3:4)
 close(10)

 end program averages
!--------------------!