1! { dg-do compile }
2! { dg-options "-O3" }
3      subroutine hesfcn(n, x, h, ldh)
4      integer n,ldh
5      double precision x(n), h(ldh)
6
7      integer i,j,k,kj
8      double precision th,u1,u2,v2
9
10      kj = 0
11      do 770 j = 1, n
12         kj = kj - j
13         do 760 k = 1, j
14            kj = kj + 1
15            v2 = 2 * x(k) - 1
16            u1 = 0
17            u2 = 2
18            do 750 i = 1, n
19               h(kj) = h(kj) + u2
20               th = 4 * v2 + u2 - u1
21               u1 = u2
22               u2 = th
23               th = v2 - 1
24  750       continue
25  760    continue
26  770 continue
27
28      end
29