c234567================================================================!
      program campiona
      implicit none
      double precision myfun, a , b , x , xi
      integer n , i

      write(*,*) 'inserire a,b,n'
      read(*,*) a, b, n

      write(*,*) 'campionamenti:'

c andrebbe anche bene
c     open(UNIT=10,FILE='campionamenti.txt',STATUS='NEW')
c le ulteriori opzioni impongono il default
c
      open(UNIT=10,FILE='campionamenti.txt',STATUS='NEW',
     +ACCESS='SEQUENTIAL', FORM='FORMATTED')

      do i=0,n
       xi = dble(i)/dble(n)
       x= a * (1.0d0 - xi) + b * xi
       write(10,*) x , myfun(x)
      enddo
      
      close(UNIT=10)
      
      stop
      end

c234567================================================================!
      double precision function myfun(x)
      implicit none
      double precision x

      if ( x .GT. 0.0d0 ) then
       myfun=dsin(x)/x
      else
       myfun=1.0d0
      endif

      return
      end