program x3
      implicit none
      integer pgopen, i
      real x(100), y(100), delta_x

      delta_x = ((3.0) - (-3.0))/(100-1)
      do i=1, 100
        x(i) = (-3.0) + (i-1)*delta_x
        y(i) = x(i)**3 + 2.0
      end do

      if (pgopen('/xwin').le.0) stop
c     call pgpap(5.0, 0.75)
      call pgenv(-3.0, 3.0, -6.0, 6.0, 0, 0)
      call pgline(100,x,y)
      call pgclos

      end