1 subroutine submaxima_tan3(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 double precision sec 5 external sec 6 f = tan(1.57d+0*sin(x)*sin(y)*sin(z)) 7 df(1) = 1.57d+0*cos(x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin 8 1 (z))**2 9 df(2) = 1.57d+0*sin(x)*cos(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin 10 1 (z))**2 11 df(3) = 1.57d+0*sin(x)*sin(y)*cos(z)*sec(1.57d+0*sin(x)*sin(y)*sin 12 1 (z))**2 13 df2(1) = 4.9298000000000003d+0*cos(x)**2*sin(y)**2*sin(z)**2*sec(1 14 1 .57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z 15 2 ))-1.57d+0*sin(x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z 16 3 ))**2 17 df2(2) = 4.9298000000000003d+0*cos(x)*sin(x)*cos(y)*sin(y)*sin(z)* 18 1 *2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin( 19 2 y)*sin(z))+1.57d+0*cos(x)*cos(y)*sin(z)*sec(1.57d+0*sin(x)*sin( 20 3 y)*sin(z))**2 21 df2(3) = 4.9298000000000003d+0*sin(x)**2*cos(y)**2*sin(z)**2*sec(1 22 1 .57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z 23 2 ))-1.57d+0*sin(x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z 24 3 ))**2 25 df2(4) = 4.9298000000000003d+0*cos(x)*sin(x)*sin(y)**2*cos(z)*sin( 26 1 z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin( 27 2 y)*sin(z))+1.57d+0*cos(x)*sin(y)*cos(z)*sec(1.57d+0*sin(x)*sin( 28 3 y)*sin(z))**2 29 df2(5) = 4.9298000000000003d+0*sin(x)**2*cos(y)*sin(y)*cos(z)*sin( 30 1 z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin( 31 2 y)*sin(z))+1.57d+0*sin(x)*cos(y)*cos(z)*sec(1.57d+0*sin(x)*sin( 32 3 y)*sin(z))**2 33 df2(6) = 4.9298000000000003d+0*sin(x)**2*sin(y)**2*cos(z)**2*sec(1 34 1 .57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z 35 2 ))-1.57d+0*sin(x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z 36 3 ))**2 37 df3(1) = 1.5479572d+1*cos(x)**3*sin(y)**3*sin(z)**3*sec(1.57d+0*si 38 1 n(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))**2-1.4 39 2 7894d+1*cos(x)*sin(x)*sin(y)**2*sin(z)**2*sec(1.57d+0*sin(x)*si 40 3 n(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))+7.73978600000 41 4 0001d+0*cos(x)**3*sin(y)**3*sin(z)**3*sec(1.57d+0*sin(x)*sin(y) 42 5 *sin(z))**4-1.57d+0*cos(x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin 43 6 (y)*sin(z))**2 44 df3(2) = 1.5479572d+1*cos(x)**2*sin(x)*cos(y)*sin(y)**2*sin(z)**3* 45 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 46 2 sin(z))**2-4.9298000000000003d+0*sin(x)**2*cos(y)*sin(y)*sin(z) 47 3 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 48 4 (y)*sin(z))+9.859600000000001d+0*cos(x)**2*cos(y)*sin(y)*sin(z) 49 5 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 50 6 (y)*sin(z))+7.739786000000001d+0*cos(x)**2*sin(x)*cos(y)*sin(y) 51 7 **2*sin(z)**3*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*sin( 52 8 x)*cos(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 53 df3(3) = 1.5479572d+1*cos(x)*sin(x)**2*cos(y)**2*sin(y)*sin(z)**3* 54 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 55 2 sin(z))**2-4.9298000000000003d+0*cos(x)*sin(x)*sin(y)**2*sin(z) 56 3 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 57 4 (y)*sin(z))+9.859600000000001d+0*cos(x)*sin(x)*cos(y)**2*sin(z) 58 5 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 59 6 (y)*sin(z))+7.739786000000001d+0*cos(x)*sin(x)**2*cos(y)**2*sin 60 7 (y)*sin(z)**3*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*cos( 61 8 x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 62 df3(4) = 1.5479572d+1*sin(x)**3*cos(y)**3*sin(z)**3*sec(1.57d+0*si 63 1 n(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))**2-1.4 64 2 7894d+1*sin(x)**2*cos(y)*sin(y)*sin(z)**2*sec(1.57d+0*sin(x)*si 65 3 n(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))+7.73978600000 66 4 0001d+0*sin(x)**3*cos(y)**3*sin(z)**3*sec(1.57d+0*sin(x)*sin(y) 67 5 *sin(z))**4-1.57d+0*sin(x)*cos(y)*sin(z)*sec(1.57d+0*sin(x)*sin 68 6 (y)*sin(z))**2 69 df3(5) = 1.5479572d+1*cos(x)**2*sin(x)*sin(y)**3*cos(z)*sin(z)**2* 70 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 71 2 sin(z))**2-4.9298000000000003d+0*sin(x)**2*sin(y)**2*cos(z)*sin 72 3 (z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 73 4 (y)*sin(z))+9.859600000000001d+0*cos(x)**2*sin(y)**2*cos(z)*sin 74 5 (z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 75 6 (y)*sin(z))+7.739786000000001d+0*cos(x)**2*sin(x)*sin(y)**3*cos 76 7 (z)*sin(z)**2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*sin( 77 8 x)*sin(y)*cos(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 78 df3(6) = 1.5479572d+1*cos(x)*sin(x)**2*cos(y)*sin(y)**2*cos(z)*sin 79 1 (z)**2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)* 80 2 sin(y)*sin(z))**2+1.47894d+1*cos(x)*sin(x)*cos(y)*sin(y)*cos(z) 81 3 *sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x) 82 4 *sin(y)*sin(z))+7.739786000000001d+0*cos(x)*sin(x)**2*cos(y)*si 83 5 n(y)**2*cos(z)*sin(z)**2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4+1 84 6 .57d+0*cos(x)*cos(y)*cos(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))** 85 7 2 86 df3(7) = 1.5479572d+1*sin(x)**3*cos(y)**2*sin(y)*cos(z)*sin(z)**2* 87 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 88 2 sin(z))**2-4.9298000000000003d+0*sin(x)**2*sin(y)**2*cos(z)*sin 89 3 (z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 90 4 (y)*sin(z))+9.859600000000001d+0*sin(x)**2*cos(y)**2*cos(z)*sin 91 5 (z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 92 6 (y)*sin(z))+7.739786000000001d+0*sin(x)**3*cos(y)**2*sin(y)*cos 93 7 (z)*sin(z)**2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*sin( 94 8 x)*sin(y)*cos(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 95 df3(8) = 1.5479572d+1*cos(x)*sin(x)**2*sin(y)**3*cos(z)**2*sin(z)* 96 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 97 2 sin(z))**2-4.9298000000000003d+0*cos(x)*sin(x)*sin(y)**2*sin(z) 98 3 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 99 4 (y)*sin(z))+9.859600000000001d+0*cos(x)*sin(x)*sin(y)**2*cos(z) 100 5 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 101 6 (y)*sin(z))+7.739786000000001d+0*cos(x)*sin(x)**2*sin(y)**3*cos 102 7 (z)**2*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*cos( 103 8 x)*sin(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 104 df3(9) = 1.5479572d+1*sin(x)**3*cos(y)*sin(y)**2*cos(z)**2*sin(z)* 105 1 sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)* 106 2 sin(z))**2-4.9298000000000003d+0*sin(x)**2*cos(y)*sin(y)*sin(z) 107 3 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 108 4 (y)*sin(z))+9.859600000000001d+0*sin(x)**2*cos(y)*sin(y)*cos(z) 109 5 **2*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin 110 6 (y)*sin(z))+7.739786000000001d+0*sin(x)**3*cos(y)*sin(y)**2*cos 111 7 (z)**2*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**4-1.57d+0*sin( 112 8 x)*cos(y)*sin(z)*sec(1.57d+0*sin(x)*sin(y)*sin(z))**2 113 df3(10) = 1.5479572d+1*sin(x)**3*sin(y)**3*cos(z)**3*sec(1.57d+0*s 114 1 in(x)*sin(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))**2-1. 115 2 47894d+1*sin(x)**2*sin(y)**2*cos(z)*sin(z)*sec(1.57d+0*sin(x)*s 116 3 in(y)*sin(z))**2*tan(1.57d+0*sin(x)*sin(y)*sin(z))+7.7397860000 117 4 00001d+0*sin(x)**3*sin(y)**3*cos(z)**3*sec(1.57d+0*sin(x)*sin(y 118 5 )*sin(z))**4-1.57d+0*sin(x)*sin(y)*cos(z)*sec(1.57d+0*sin(x)*si 119 6 n(y)*sin(z))**2 120 end 121 double precision function sec(x) 122 implicit double precision (a-z) 123 sec=1.d0/cos(x) 124 return 125 end 126