program circle
      i_dummy=pgopen('/xwin')
      xmin=-10
      xmax=10
      ymin=-10
      ymax=10
      call pgenv(xmin,xmax,ymin,ymax,1,0)
      call pgsfs(3)
      call pgsci(3)
c     do i=1,16
c       call pgsci(i)
c       write(*,*)'Now color is ',i
        x=0
        y=0
        r=4
        call pgcirc(x,y,r)
c       do j=1,1000000000
c         z=sqrt(36.0)
c       end do
c     end do
      call pgclos
      end