o1 = cos(phi)
        o2 = cos(q)
        o3 = o2**2
        o4 = -3.00000000000000d0*eta
        o5 = exp(o4)
        o6 = sin(phi)
        o7 = sin(q)
        o8 = o7**2
        o9 = dnnnnngtil
        o10 = -3.6000000000000d1*o9
        o11 = dnnnnnnngtil
        o12 = -4.0000000000000d0*o11*o8
        o13 = 1.16000000000000d2*o8*o9
        o14 = dnnngtil
        o15 = -4.1500000000000d2*o8
        o16 = 1.50000000000000d2 + o15
        o17 = o14*o16
        o18 = dngtil
        o19 = 5.0000000000000d0*o8
        o20 = -2.00000000000000d0 + o19
        o21 = 7.5000000000000d1*o18*o20
        o22 = o10 + o12 + o13 + o17 + o21
        o23 = dnnnnnngtil
        o24 = 3.00000000000000d0*o23
        o25 = 2.00000000000000d0*q
        o26 = cos(o25)
        o27 = 4.0000000000000d0*o23*o26
        o28 = 4.0000000000000d0*q
        o29 = cos(o28)
        o30 = -7.0000000000000d0*o23*o29
        o31 = dnnnngtil
        o32 = 2.10000000000000d1*o31
        o33 = 6.4000000000000d1*o26*o31
        o34 = 2.03000000000000d2*o29*o31
        o35 = o24 + o27 + o30 + o32 + o33 + o34
        o36 = -4.0000000000000d0*o35
        o37 = 2.68000000000000d2*o26
        o38 = 5.8100000000000d2*o29
        o39 = 1.11000000000000d2 + o37 + o38
        o40 = dnngtil
        o41 = 5.0000000000000d0*o39*o40
        o42 = 2.00000000000000d1*o26
        o43 = 3.5000000000000d1*o29
        o44 = 9.0000000000000d0 + o42 + o43
        o45 = gtil
        o46 = -7.5000000000000d1*o44*o45
        o47 = o36 + o41 + o46
        o48 = o1**2
        o49 = o6**2
        o50 = o7*o8
        hxx(i,j,k) = 6.6666666666667d-1*o1*o22*o3*o5*o6*o8 - 8.333333333
     &  3333d-2*o1*o47*o5*o6*o8
        hxy(i,j,k) = -3.3333333333333d-1*o22*o3*o48*o5*o8 + 4.1666666666
     &  667d-2*o47*o48*o5*o8 + 3.3333333333333d-1*o22*o3*o49*o5*o8 - 4.1
     &  666666666667d-2*o47*o49*o5*o8
        hxz(i,j,k) = -3.3333333333333d-1*o2*o22*o5*o50*o6 - 4.1666666666
     &  667d-2*o2*o47*o5*o6*o7
        hyy(i,j,k) = -6.6666666666667d-1*o1*o22*o3*o5*o6*o8 + 8.33333333
     &  33333d-2*o1*o47*o5*o6*o8
        hyz(i,j,k) = 3.3333333333333d-1*o1*o2*o22*o5*o50 + 4.16666666666
     &  67d-2*o1*o2*o47*o5*o7
        hzz(i,j,k) = 0