o1 = cos(q)
        o2 = 2.00000000000000d0*q
        o3 = cos(o2)
        o4 = dnnnnngtil
        o5 = sin(q)
        o6 = o5**2
        o7 = o5*o6
        o8 = dnnngtil
        o9 = dngtil
        o10 = dnnnnnnngtil
        o11 = o6**2
        o12 = 4.0000000000000d0*q
        o13 = cos(o12)
        o14 = dnnnngtil
        o15 = dnngtil
        o16 = o1**2
        o17 = dnnnnnngtil
        o18 = o11*o6
        o19 = gtil
        o20 = byJ**2
        o21 = o14**2
        o22 = o3**2
        o23 = o13**2
        o24 = o15**2
        o25 = o4**2
        o26 = o8**2
        o27 = o9**2
        o28 = o16**2
        o29 = o17**2
        o30 = o10**2
        o31 = o11**2
        o32 = o31*o6
        o33 = o19**2
c        exc33(i,j,k) = 0
c        exc32(i,j,k) = 1.33333333333333d0*o1*o10*o11*o5 - 2.083333333333
c     &  33d-2*o1*(3.5200000000000d2 - 9.2800000000000d2*o3)*o4*o7 + 1.66
c     &  666666666667d-1*o1*(1.15000000000000d2 - 4.1500000000000d2*o3)*o
c     &  7*o8 - 1.56250000000000d0*o1*(8.0000000000000d0 - 4.000000000000
c     &  0d1*o3)*o7*o9
c        exc31(i,j,k) = -8.0000000000000d0*o11*o16*o17 + 1.33333333333333
c     &  d0*o17*o18 + 3.00000000000000d0*byJ*o6 - 2.08333333333333d-1*o15*(
c     &  -1.11000000000000d2 - 5.8100000000000d2*o13 - 2.68000000000000d2
c     &  *o3)*o6 - 5.2083333333333d-3*o14*(6.7200000000000d2 + 6.49600000
c     &  00000d3*o13 + 2.04800000000000d3*o3)*o6 + o19*(5.0000000000000d1
c     &  *o11 + 7.5000000000000d2*o11*o16 - 1.25000000000000d2*o18 - 2.00
c     &  000000000000d2*o16*o6)
c        exc22(i,j,k) = 0
c        exc21(i,j,k) = 0
c        exc11(i,j,k) = 0
        t1 = -9.6000000000000d1*byJ*o11*o16*o17 + 1.12000000000000d2*o11*o
     &  14*o16*o17 + 1.08266666666667d3*o11*o13*o14*o16*o17 - 7.40000000
     &  00000d2*o11*o15*o16*o17 - 3.8733333333333d3*o11*o13*o15*o16*o17 
     &  + 1.60000000000000d1*byJ*o17*o18 - 1.86666666666667d1*o14*o17*o18 
     &  - 1.80444444444444d2*o13*o14*o17*o18 + 1.23333333333333d2*o15*o1
     &  7*o18 + 6.4555555555556d2*o13*o15*o17*o18 + 6.0000000000000d2*byJ*
     &  o11*o19 - 7.0000000000000d2*o11*o14*o19 - 6.7666666666667d3*o11*
     &  o13*o14*o19 + 4.6250000000000d3*o11*o15*o19 + 2.42083333333333d4
     &  *o11*o13*o15*o19 + 9.0000000000000d3*byJ*o11*o16*o19 - 1.050000000
     &  00000d4*o11*o14*o16*o19 - 1.01500000000000d5*o11*o13*o14*o16*o19
     &   + 6.9375000000000d4*o11*o15*o16*o19 + 3.6312500000000d5*o11*o13
     &  *o15*o16*o19 - 1.50000000000000d3*byJ*o18*o19 + 1.75000000000000d3
     &  *o14*o18*o19 + 1.69166666666667d4*o13*o14*o18*o19 - 1.1562500000
     &  0000d4*o15*o18*o19 - 6.0520833333333d4*o13*o15*o18*o19 - 2.66666
     &  666666667d3*o16*o17*o18*o19 + 1.07555555555556d2*o11*o16*o25
        t2 = 7.4755555555556d2*o11*o16*o22*o25 + 7.3472222222222d2*o11*o
     &  16*o26 + 9.5680555555556d3*o11*o16*o22*o26 + 3.12500000000000d2*
     &  o11*o16*o27 + 7.8125000000000d3*o11*o16*o22*o27 + 6.400000000000
     &  0d3*o11*o17*o19*o28 - 2.40000000000000d4*o17*o18*o19*o28 + 1.280
     &  00000000000d2*o18*o28*o29 + 3.4133333333333d2*o11*o14*o16*o17*o3
     &   - 1.78666666666667d3*o11*o15*o16*o17*o3 - 5.6888888888889d1*o14
     &  *o17*o18*o3 + 2.97777777777778d2*o15*o17*o18*o3 - 2.133333333333
     &  33d3*o11*o14*o19*o3 + 1.11666666666667d4*o11*o15*o19*o3 - 3.2000
     &  000000000d4*o11*o14*o16*o19*o3 + 1.67500000000000d5*o11*o15*o16*
     &  o19*o3 + 5.3333333333333d3*o14*o18*o19*o3 - 2.79166666666667d4*o
     &  15*o18*o19*o3 - 5.6711111111111d2*o11*o16*o25*o3 - 5.30277777777
     &  78d3*o11*o16*o26*o3 - 3.12500000000000d3*o11*o16*o27*o3 + 2.6666
     &  6666666667d2*o17*o19*o31 + 8.0000000000000d3*o16*o17*o19*o31 - 4
     &  .2666666666667d1*o16*o29*o31 + 3.5555555555556d0*o16*o30*o31 - 6
     &  .6666666666667d2*o17*o19*o32 + 3.5555555555556d0*o29*o32 - 4.000
     &  0000000000d4*o11*o16*o33
        t3 = 5.0000000000000d3*o18*o33 + 2.50000000000000d5*o16*o18*o33 
     &  - 6.0000000000000d5*o11*o28*o33 + 1.12500000000000d6*o18*o28*o33
     &   - 2.50000000000000d4*o31*o33 - 3.7500000000000d5*o16*o31*o33 + 
     &  3.12500000000000d4*o32*o33 - 3.9111111111111d1*o10*o16*o18*o4 + 
     &  1.03111111111111d2*o10*o16*o18*o3*o4 - 4.2000000000000d1*byJ*o14*o
     &  6 - 4.0600000000000d2*byJ*o13*o14*o6 + 2.77500000000000d2*byJ*o15*o6
     &   + 1.45250000000000d3*byJ*o13*o15*o6 - 3.2375000000000d2*o14*o15*o
     &  6 - 4.8241666666667d3*o13*o14*o15*o6 - 2.40000000000000d3*byJ*o16*
     &  o19*o6 + 2.80000000000000d3*o14*o16*o19*o6 + 2.70666666666667d4*
     &  o13*o14*o16*o19*o6 - 1.85000000000000d4*o15*o16*o19*o6 - 9.68333
     &  33333333d4*o13*o15*o16*o19*o6 + 1.80000000000000d1*o20*o6 + 2.45
     &  000000000000d1*o21*o6 + 4.7366666666667d2*o13*o21*o6 - 2.3822222
     &  2222222d3*o14*o15*o22*o6 + 2.27555555555556d2*o21*o22*o6 - 1.638
     &  09722222222d4*o14*o15*o23*o6 + 2.28938888888889d3*o21*o23*o6
        t4 = 1.06953125000000d3*o24*o6 + 1.11963541666667d4*o13*o24*o6 +
     &   6.2347222222222d3*o22*o24*o6 + 2.93021701388889d4*o23*o24*o6 - 
     &  1.28000000000000d2*byJ*o14*o3*o6 + 6.7000000000000d2*byJ*o15*o3*o6 -
     &   1.76833333333333d3*o14*o15*o3*o6 - 1.27205555555556d4*o13*o14*o
     &  15*o3*o6 + 8.5333333333333d3*o14*o16*o19*o3*o6 - 4.4666666666667
     &  d4*o15*o16*o19*o3*o6 + 1.49333333333333d2*o21*o3*o6 + 1.44355555
     &  555556d3*o13*o21*o3*o6 + 5.1645833333333d3*o24*o3*o6 + 2.7032638
     &  8888889d4*o13*o24*o3*o6 + 8.0000000000000d4*o28*o33*o6 + 1.02222
     &  222222222d2*o10*o16*o18*o8 - 3.6888888888889d2*o10*o16*o18*o3*o8
     &   - 5.6222222222222d2*o11*o16*o4*o8 - 5.3488888888889d3*o11*o16*o
     &  22*o4*o8 + 3.5111111111111d3*o11*o16*o3*o4*o8 - 6.6666666666667d
     &  1*o10*o16*o18*o9 + 3.3333333333333d2*o10*o16*o18*o3*o9 + 3.66666
     &  66666667d2*o11*o16*o4*o9 + 4.8333333333333d3*o11*o16*o22*o4*o9 -
     &   2.80000000000000d3*o11*o16*o3*o4*o9 - 9.5833333333333d2*o11*o16
     &  *o8*o9 - 1.72916666666667d4*o11*o16*o22*o8*o9 + 8.2500000000000d
     &  3*o11*o16*o3*o8*o9
        ksq(i,j,k) = t1 + t2 + t3 + t4