program sayilar
integer, parameter::satir=5,sutun=3
real,dimension(satir,sutun)::deger
real,dimension(sutun)::ortalama
integer::i,j
open(unit=1,file="sayilar.txt",status="old",  action="read")
do i=1,satir
   read(unit=1,fmt=*)deger(i,1:sutun)
   print*,deger(i,1:sutun)
end do

do j=1,sutun
   ortalama(j)=(sum(deger(1:satir,j)))/satir
end do
open(unit=2,file="ortalama.ort",status="replace",action="write")
do i=1,sutun
  write(unit=2,fmt="(f9.3)")ortalama(i)
enddo 
end program sayilar