! { dg-do compile }
! { dg-options "-O3" }
      subroutine hesfcn(n, x, h, ldh)
      integer n,ldh
      double precision x(n), h(ldh)

      integer i,j,k,kj
      double precision th,u1,u2,v2
 
      kj = 0
      do 770 j = 1, n
         kj = kj - j
         do 760 k = 1, j
            kj = kj + 1
            v2 = 2 * x(k) - 1
            u1 = 0
            u2 = 2
            do 750 i = 1, n
               h(kj) = h(kj) + u2
               th = 4 * v2 + u2 - u1
               u1 = u2
               u2 = th
               th = v2 - 1
  750       continue
  760    continue
  770 continue

      end