1      subroutine submaxima_exp3(x,y,z,f,df,df2,df3)
2      implicit none
3      double precision x,y,z,f,df(3),df2(3*4/2),df3(3*4*5/6)
4      f = exp(sin(z)+sin(y)+sin(x))
5      df(1) = cos(x)*exp(sin(z)+sin(y)+sin(x))
6      df(2) = cos(y)*exp(sin(z)+sin(y)+sin(x))
7      df(3) = exp(sin(z)+sin(y)+sin(x))*cos(z)
8      df2(1) = cos(x)**2*exp(sin(z)+sin(y)+sin(x))-sin(x)*exp(sin(z)+sin
9     1   (y)+sin(x))
10      df2(2) = cos(x)*cos(y)*exp(sin(z)+sin(y)+sin(x))
11      df2(3) = cos(y)**2*exp(sin(z)+sin(y)+sin(x))-sin(y)*exp(sin(z)+sin
12     1   (y)+sin(x))
13      df2(4) = cos(x)*exp(sin(z)+sin(y)+sin(x))*cos(z)
14      df2(5) = cos(y)*exp(sin(z)+sin(y)+sin(x))*cos(z)
15      df2(6) = exp(sin(z)+sin(y)+sin(x))*cos(z)**2-exp(sin(z)+sin(y)+sin
16     1   (x))*sin(z)
17      df3(1) = -3*cos(x)*sin(x)*exp(sin(z)+sin(y)+sin(x))+cos(x)**3*exp(
18     1   sin(z)+sin(y)+sin(x))-cos(x)*exp(sin(z)+sin(y)+sin(x))
19      df3(2) = cos(x)**2*cos(y)*exp(sin(z)+sin(y)+sin(x))-sin(x)*cos(y)*
20     1   exp(sin(z)+sin(y)+sin(x))
21      df3(3) = cos(x)*cos(y)**2*exp(sin(z)+sin(y)+sin(x))-cos(x)*sin(y)*
22     1   exp(sin(z)+sin(y)+sin(x))
23      df3(4) = -3*cos(y)*sin(y)*exp(sin(z)+sin(y)+sin(x))+cos(y)**3*exp(
24     1   sin(z)+sin(y)+sin(x))-cos(y)*exp(sin(z)+sin(y)+sin(x))
25      df3(5) = cos(x)**2*exp(sin(z)+sin(y)+sin(x))*cos(z)-sin(x)*exp(sin
26     1   (z)+sin(y)+sin(x))*cos(z)
27      df3(6) = cos(x)*cos(y)*exp(sin(z)+sin(y)+sin(x))*cos(z)
28      df3(7) = cos(y)**2*exp(sin(z)+sin(y)+sin(x))*cos(z)-sin(y)*exp(sin
29     1   (z)+sin(y)+sin(x))*cos(z)
30      df3(8) = cos(x)*exp(sin(z)+sin(y)+sin(x))*cos(z)**2-cos(x)*exp(sin
31     1   (z)+sin(y)+sin(x))*sin(z)
32      df3(9) = cos(y)*exp(sin(z)+sin(y)+sin(x))*cos(z)**2-cos(y)*exp(sin
33     1   (z)+sin(y)+sin(x))*sin(z)
34      df3(10) = -3*exp(sin(z)+sin(y)+sin(x))*cos(z)*sin(z)+exp(sin(z)+si
35     1   n(y)+sin(x))*cos(z)**3-exp(sin(z)+sin(y)+sin(x))*cos(z)
36      end
37