program prog20
      implicit none
      double precision T(10), myF
      integer i
      external myF
      
      call tab(myF, 0.d0, 1.d0, 10, T)
      
      do i=1,10
        print *, i, T(i) 
      end do
      
      end
      
      double precision function myF(x)
      implicit none
      double precision x
        myF = sin(x) + cos(x)
      end