PROGRAM GPE3	! Demonstrate Log graphics routine usage
      DIMENSION X(100),Y(100)
      DATA XORIG,DX,XMAX/0.0,100.,1000./,YORIG,DY,YMAX/0.0,100.0,1000.0/

!  Compute a quadratic curve
      DO I=1,100
        X(I) = I*10.0
        Y(I) = 0.1*I*I
      END DO
      IDX = (XMAX-XORIG)/9.0 + 1
      DX = IDX
!		Draw First Plot
      CALL GSTART
      CALL PAGE(11.0,8.5)          ! Specify page size
      CALL PHYSOR(1.0,0.75)	   ! PHYSOR must be called before AREA2D
      CALL AREA2D(4.0,6.5)         ! Define Subplot area
      CALL XNAME('X LABEL$',100)   ! Define X Axis label
      CALL YNAME('Y LABEL$',100)   ! Define Y Axis label
      CALL HEADIN('Example Plot 3$',100,1.5,2)
      CALL HEADIN('3A:  Use of XLOG$',100,1.5,2)
      XORIG = 1.0
      XCYCLE = 4.0
      IDY = (YMAX-YORIG)/6.5 + 1
      DY = IDY
      CALL YINTAX
      CALL XLOG(XORIG,XCYCLE,YORIG,DY)
      CALL CURVE(X,Y,100,5)
      CALL ENDGR
!		Draw Second Plot
      CALL OREL(5.0,-0.3)
      CALL AREA2D(4.0,3.0)         ! Define Subplot area
      CALL HEADIN('3B:  Use of YLOG$',100,1.5,1)
      YORIG = 0.1
      YCYCLE = 4.0
      XORIG = 0.0
      CALL XINTAX
      CALL YLOG(XORIG,DX,YORIG,YCYCLE)
      CALL CURVE(X,Y,100,5)
!		Draw Third Plot
      CALL OREL(0.0,4.25)
      CALL AREA2D(4.0,3.0)         ! Define Subplot area
      CALL HEADIN('3C:  Use of LOGLOG$',100,1.5,1)
      YORIG = 0.1
      YCYCLE = 4.0
      XORIG = 1.0
      XCYCLE = 3.0
      CALL LOGLOG(XORIG,XCYCLE,YORIG,YCYCLE)
      CALL CURVE(X,Y,100,5)
      CALL ENDPL(0)
      CALL GSTOP
 
      CALL EXIT
      END