program x2sinx
      i_dummy=pgopen('/xwin')
      ymin=-100.0
      ymax= 100.0
c n is the number of points
      n=1001
      xmin=-10.0
      xmax=10.0
      dx=(xmax-xmin)/(n-1)
      call pgenv(xmin,xmax,ymin,ymax,0,0)
      do i=1,n
        x=xmin+(i-1)*dx
        y=(x**2)*sin(x)
        call pgpt(1,x,y,20)
      end do
      call pgclos
      end