      Virtual x(41), y(41)
c
      do 10 i = 1, 41
         x(i) = .05 * (i - 1)
         y(i) = sin(3.141592654 * x(i))
10    continue
c
      call iniplt(7, 10., 7.)
      call scale(0., 2., -1., 1.)
      call axis(.5, .5, 'x', 1, 2, 1, 'sin', 3, 2, 1)
      call Vline(x, y, 41, 0, 0, 0, 0)
      call endplt
c
      end
