Program Moduli

implicit none
real :: vx,vy,temp,v,somma,vm
integer :: i,iend

! Legge un file di coppie di numeri (una coppia per riga)
! Supponendo che siano le componenti X e Y di un vettore
! calcola i moduli
! calcola la media dei moduli
!
! nella istruzione read iostat=iend carica nella variabile iend
! 0 se lettura corretta
! un valore negativo se EOF (End Of File)
! un valore positivo se si e' verificato un errore 
!
! ( quando leggo da tastiera EOF e' rappresentato da un Ctrl D )

somma=0
do i=1,100000
   read(5,*,iostat=iend)vx,vy
   if(iend /= 0)exit
   v=sqrt(vx**2+vy**2)
   somma=somma+v
   write(6,'(f10.2)')v
enddo

if(iend > 0)then
   write(6,*)'Errore di lettura :',iend
else
   vm=somma/(i-1)
   write(6,*) 'vmedia=',vm
endif
end