subroutine squarf(x, n, result)
      integer n
      double precision x(n), result(n)
      integer i

      do 10 i = 1, n
          result(i) = x(i) ** 2
   10 continue

      end