c2345678901234567890
      program test2
      integer pgopen
      real x, y, x2, y2, angle, delta
      if (pgopen('/xwin') .le. 0) stop
      call pgpap (6.0, 0.75)
      call pgenv (0.0, 1.0, 0.0, 1.0, 1, 1)
      call pglab ('x-axis', 'y-axis', 'My test PGPLOT graph')
      call pgsfs (2)
      call pgcirc (0.5, 0.5, 0.1)
      
      angle = 0.0
      delta = 3.14159268/36

      x = 0.2 * cos(angle) + 0.5
      y = 0.2 * sin(angle) + 0.5
      call pgcirc (x, y, 0.02)

 1000 continue
      angle = angle + delta
      x2 = 0.2 * cos(angle) + 0.5
      y2 = 0.2 * sin(angle) + 0.5

      call pgbbuf

      call pgsci(0)
      call pgcirc (x, y, 0.02)
      call pgsci(1)
      call pgcirc (x2, y2, 0.02)

      call pgebuf

      x = x2
      y = y2

      do i = 1, 1000000
       idummy = i**3
      enddo

      goto 1000

      call pgclos
c     write (*,*) 'Hello, world!'
      end