program sine_static
c     implicit none
      real x(100), y(100)
      integer pgopen

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

      x_max = 5*(2*3.1415)
      delta_theta = x_max/100

      do i=1,100
        x(i) = i*delta_theta
        y(i) = sin(x(i))
      end do

      call pgenv(0.0,x_max,-1.1,1.1,0,0)
      call pgline(100,x,y)

      call pgclos
      end