c234567================================================================!
      program integra
      implicit none
      double precision xnew,ynew,xold,yold,accumulo
      integer i


c accesso sequenziale, formato formatted di default
      open(UNIT=10,FILE='campionamenti.txt',STATUS='OLD')

c inizializzo variabile di accumulo a valore nullo
      accumulo=0.0d0
c inizializzo variabile contatore campionamenti letti
      i=0

c ciclo sommando contributi alla variabile di accumulo

   20 continue
      read(UNIT=10, FMT=* , END=100) xnew , ynew
      i=i+1
      if ( i .GE. 2 ) then
        accumulo = accumulo + (xnew-xold) * (ynew+yold)/2.0d0
        write(*,*) i
      endif
      xold=xnew
      yold=ynew
      goto 20

c finiti i contributi, stampo a video i risultati

  100 write(*,*) accumulo

      stop
      end