      dimension f(21,11), x(21,11), y(21,11)
1     format(e18.10)
c
      pi = 3.141592654
      do 20 j = 1, 11
         xx = .1 * (j - 1)
         do 10 i = 1, 21
            yy = .1 * (i - 1)
            f(i,j) = cos(pi * xx) * sin(pi * yy)
10       continue
20    continue
100   type *, 'input xv, yv, zv'
      read(5, *, end = 1500) xv, yv, zv
      call iniplt(7, 10., 7.)
      call plot3d(f, x, y, 21, 11, 10., 5., xv, yv, zv, 0, 0.)
      call endplt
      goto 100
c
1500  end
