program write_sin_to_file

      open(unit=10 ,file='sin.dat')
      xmax=10.0
      xmin=0.0
      n=1000
      write(10,*)n
      dx=(xmax-xmin)/n
      do i=1,n
        x=(i-1)*dx
        y=sin(x)
        write(10,*)x,y
      end do
      close(10)
      write(*,*)'File written to sin.dat'
      end