program plot_sin
      real xx(1000),yy(1000)
      i_dummy=pgopen('/xwin')
      xleft=-14.0
      xright=14.0

      call pgenv(xleft,xright,-2.0,2.0,0,0)

      m=1000
      dx=(xright-xleft)/(m-1)
      do i=1,m
        xx(i)=xleft+(i-1)*dx
        yy(i)=sin(xx(i))
      end do        
      call pgline(m,xx,yy)
      call pgclos
      end