program tab_gen

      implicit none
      integer np, n, i
      parameter (np=50)
      real x(np), y(np), x_min, x_max, delta_x
      write(*,*) 'How many tabulated point to be generated ?
     &(less than 50)'

      read(*,*) n

      write(*,*) 'Range of x: x1 and x2'
      read(*,*) x_min, x_max 
      delta_x = (x_max - x_min)/(n-1)
      do i=1,n
        x(i) = x_min + (i-1)*delta_x
        y(i) = sin(x(i))
      end do 

      open (unit=10,file='tab_gen.dat')
      do i = 1,n
        write(*,*) x(i), y(i)
        write(10,*) x(i),y(i)  
      end do

      end