program plot_with_auto_range
      integer pgopen
      real x(101),y(101),xold(101),yold(101)

      if( pgopen('/xwin') <= 0 ) stop

      xini = -5.0

      dx = 10.0/100
      do i=1,101
        x(i)= xini + (i-1)*dx
c       y(i)= x(i)**2
        y(i)= sin(x(i))
      end do

      xmin=x(1)
      xmax=x(1)
      ymin=y(1)
      ymax=y(1)
      do i=2,101
        if(x(i)>xmax)xmax=x(i)
        if(x(i)<xmin)xmin=x(i)
        if(y(i)>ymax)ymax=y(i)
        if(y(i)<ymin)ymin=y(i)
      end do

      call pgenv(xmin*1.01,xmax*1.01,ymin*1.2,ymax*1.2,0,0)
      call pgline(101,x,y)

      do 1234, j=1,10000

      do i=1,101
        xold(i)=x(i)
        yold(i)=y(i)
      end do
      x_shift = x_shift+0.001 

      do i=1,101
        x(i)= xini + (i-1)*dx
c       y(i)= (x(i)+x_shift)**2
        y(i)= sin(x(i)+x_shift)
      end do

      call pgbbuf

      call pgsci(0)
      call pgline(101,xold,yold)
      call pgsci(1)
      call pgline(101,x,y)

      call pgebuf

 1234 continue

      call pgclos

      end