1 /* Automatically generated functional code: wl90c
2 Maxima input:
3 >> PI: 3.14159265358979323846;
4 >> rho: rhoa + rhob;
5 >> grad: sqrt(grada*grada + gradb*gradb + 2*gradab);
6 >> rho13: rho^(1/3);
7 >>
8 >> a: -0.74860;
9 >> b: 0.06001;
10 >> c: 3.60073;
11 >> d: 0.90000;
12 >>
13 >> rs: (3/(4*PI*rho))^(1/3);
14 >> zeta: (rhoa-rhob)/(rhoa+rhob);
15 >> xa: grada/rhoa^(4/3);
16 >> xb: gradb/rhob^(4/3);
17 >>
18 >> t1: a*rho+b*grad/rho^(1/3);
19 >> t2: (1-zeta^2)^(1/2);
20 >> t3: c+d*(xa+xb)+rs;
21 >>
22 >> K(rhoa,rhob,grada,gradb,gradab):=t1*t2/t3;
23 >>
24 */
25
26 // add "extern Functional wl90cFunctional;" to 'functionals.h'
27 // add "&wl90cFunctional," to 'functionals.c'
28 // add "fun-wl90c.c" to 'Makefile.in'
29
30 #include <math.h>
31 #include <stddef.h>
32 #include "general.h"
33
34 #define __CVERSION__
35
36 #include "functionals.h"
37 #define LOG log
38 #define ABS fabs
39 #define ASINH asinh
40 #define SQRT sqrt
41
42 static const real WL90C_THRESHOLD = 1e-14;
43
44
45 /* INTERFACE PART */
wl90c_isgga(void)46 static integer wl90c_isgga(void) {return 1;}
47 static integer wl90c_read(const char* conf_line);
48 static real wl90c_energy(const FunDensProp* dp);
49 static void wl90c_first(FunFirstFuncDrv *ds, real factor,
50 const FunDensProp* dp);
51 static void wl90c_second(FunSecondFuncDrv *ds, real factor,
52 const FunDensProp* dp);
53 static void wl90c_third(FunThirdFuncDrv *ds, real factor,
54 const FunDensProp* dp);
55
56 static void wl90c_fourth(FunFourthFuncDrv *ds, real factor,
57 const FunDensProp* dp);
58
59 //static integer fun_true(void) { return 1; }
60 Functional WL90cFunctional = {
61 "WL90c",
62 wl90c_isgga,
63 3,
64 wl90c_read,
65 NULL,
66 wl90c_energy,
67 wl90c_first,
68 wl90c_second,
69 wl90c_third,
70 wl90c_fourth
71 };
72
73 /* IMPLEMENTATION PART */
74 static integer
wl90c_read(const char * conf_line)75 wl90c_read(const char* conf_line)
76 {
77 fun_set_hf_weight(0);
78 return 1;
79 }
80
81
82 static real
wl90c_energy(const FunDensProp * dp)83 wl90c_energy(const FunDensProp* dp)
84 {
85 real t[3],zk;
86 real rhoa = dp->rhoa;
87 real rhob = dp->rhob;
88 real grada = dp->grada;
89 real gradb = dp->gradb;
90 real gradab = dp->gradab;
91
92 if (rhoa + rhob > WL90C_THRESHOLD) {
93
94 t[1] = rhob+rhoa;
95 t[2] = 1/pow(t[1],.3333333333333333);
96 zk = (0.06001*t[2]*sqrt(pow(gradb,2.0)+2.0*gradab+pow(grada,2.0))-0.7486*t[1])*sqrt(1.0-1.0*pow(rhoa-1.0*rhob,2.0)/pow(t[1],2.0))/(0.9*(gradb/pow(rhob,1.333333333333333)+grada/pow(rhoa,1.333333333333333))+0.6203504908994*t[2]+3.60073);
97 return zk;
98
99 } else {
100 return 0;
101 }
102 }
103
104 static void
wl90c_first(FunFirstFuncDrv * ds,real factor,const FunDensProp * dp)105 wl90c_first(FunFirstFuncDrv *ds, real factor, const FunDensProp* dp)
106 {
107 real t[20];
108 real dfdra, dfdrb, dfdga, dfdgb, dfdab;
109 real rhoa = dp->rhoa;
110 real rhob = dp->rhob;
111 real grada = dp->grada;
112 real gradb = dp->gradb;
113 real gradab = dp->gradab;
114
115 if (rhoa + rhob > WL90C_THRESHOLD) {
116
117 t[1] = rhoa-1.0*rhob;
118 t[2] = pow(t[1],2.0);
119 t[3] = rhob+rhoa;
120 t[4] = 1/pow(t[3],2.0);
121 t[5] = sqrt(1.0-1.0*t[2]*t[4]);
122 t[6] = 1/pow(t[3],1.333333333333333);
123 t[7] = -.2067834969664666*t[6];
124 t[8] = 1/pow(rhoa,1.333333333333333);
125 t[9] = 1/pow(rhob,1.333333333333333);
126 t[10] = 1/pow(t[3],.3333333333333333);
127 t[11] = 0.6203504908994*t[10]+0.9*(gradb*t[9]+grada*t[8])+3.60073;
128 t[12] = 1/pow(t[11],2.0);
129 t[13] = sqrt(pow(gradb,2.0)+2.0*gradab+pow(grada,2.0));
130 t[14] = 0.06001*t[13]*t[10]-0.7486*t[3];
131 t[15] = 1/t[5];
132 t[16] = 2.0*t[2]/pow(t[3],3.0);
133 t[17] = 1/t[11];
134 t[18] = t[5]*(-.02000333333333333*t[13]*t[6]-0.7486)*t[17];
135 t[19] = 1/t[13];
136 dfdra = -1.0*t[12]*t[14]*t[5]*(t[7]-1.2*grada/pow(rhoa,2.333333333333333))+0.5*t[14]*t[15]*t[17]*(t[16]-2.0*t[1]*t[4])+t[18];
137 dfdrb = -1.0*t[12]*t[14]*t[5]*(t[7]-1.2*gradb/pow(rhob,2.333333333333333))+0.5*t[14]*t[15]*t[17]*(2.0*t[1]*t[4]+t[16])+t[18];
138 dfdga = 0.06001*grada*t[19]*t[10]*t[5]*t[17]-0.9*t[8]*t[5]*t[12]*t[14];
139 dfdgb = 0.06001*gradb*t[19]*t[10]*t[5]*t[17]-0.9*t[9]*t[5]*t[12]*t[14];
140 dfdab = 0.06001*t[19]*t[10]*t[5]*t[17];
141 ds->df1000 += factor*dfdra;
142 ds->df0100 += factor*dfdrb;
143 ds->df0010 += factor*dfdga;
144 ds->df0001 += factor*dfdgb;
145 ds->df00001 += factor*dfdab;
146 };
147 }
148
149 static void
wl90c_second(FunSecondFuncDrv * ds,real factor,const FunDensProp * dp)150 wl90c_second(FunSecondFuncDrv *ds, real factor, const FunDensProp* dp)
151 {
152 real t[44];
153 real dfdra, dfdrb, dfdga, dfdgb, dfdab;
154 real d2fdraga, d2fdrara, d2fdrarb, d2fdragb, d2fdrbrb;
155 real d2fdrbgb, d2fdgaga, d2fdgbgb, d2fdrbga;
156 real d2fdraab, d2fdrbab;
157 real d2fdgaab, d2fdgbab, d2fdabab, d2fdgagb;
158 real rhoa = dp->rhoa;
159 real rhob = dp->rhob;
160 real grada = dp->grada;
161 real gradb = dp->gradb;
162 real gradab = dp->gradab;
163
164 if (rhoa + rhob > WL90C_THRESHOLD) {
165
166 t[1] = rhoa-1.0*rhob;
167 t[2] = pow(t[1],2.0);
168 t[3] = rhob+rhoa;
169 t[4] = 1/pow(t[3],2.0);
170 t[5] = sqrt(1.0-1.0*t[2]*t[4]);
171 t[6] = 1/pow(rhoa,2.333333333333333);
172 t[7] = 1/pow(t[3],1.333333333333333);
173 t[8] = -.2067834969664666*t[7];
174 t[9] = t[8]-1.2*grada*t[6];
175 t[10] = 1/pow(rhoa,1.333333333333333);
176 t[11] = 1/pow(rhob,1.333333333333333);
177 t[12] = 1/pow(t[3],.3333333333333333);
178 t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;
179 t[14] = 1/pow(t[13],2.0);
180 t[15] = pow(grada,2.0);
181 t[16] = pow(gradb,2.0);
182 t[17] = sqrt(2.0*gradab+t[16]+t[15]);
183 t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];
184 t[19] = 1/t[5];
185 t[20] = 1/pow(t[3],3.0);
186 t[21] = 2.0*t[2]*t[20];
187 t[22] = t[21]-2.0*t[1]*t[4];
188 t[23] = 1/t[13];
189 t[24] = -.02000333333333333*t[17]*t[7]-0.7486;
190 t[25] = t[5]*t[24]*t[23];
191 t[26] = 1/pow(rhob,2.333333333333333);
192 t[27] = t[8]-1.2*gradb*t[26];
193 t[28] = 2.0*t[1]*t[4]+t[21];
194 t[29] = 1/t[17];
195 t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];
196 t[31] = 1/pow(t[13],3.0);
197 t[32] = 1/pow(t[3],2.333333333333333);
198 t[33] = .2757113292886222*t[32];
199 t[34] = -6.0*t[2]/pow(t[3],4.0);
200 t[35] = -2.0*t[4];
201 t[36] = 1/pow(t[5],3.0);
202 t[37] = .02667111111111111*t[17]*t[32]*t[5]*t[23];
203 t[38] = -0.9*t[10]*t[5]*t[24]*t[14];
204 t[39] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];
205 t[40] = -0.9*t[11]*t[5]*t[24]*t[14];
206 t[41] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];
207 t[42] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];
208 t[43] = 1/pow(t[17],3.0);
209 dfdra =-1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];
210 dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];
211 dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];
212 dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];
213 dfdab = t[30];
214 d2fdrara = -1.0*t[14]*t[18]*t[5]*(2.8*grada/pow(rhoa,3.333333333333333)+t[33])+2.0*t[18]*t[31]*t[5]*pow(t[9],2.0)-2.0*t[14]*t[24]*t[5]*t[9]-1.0*t[14]*t[18]*t[19]*t[22]*t[9]+t[37]-0.25*t[18]*pow(t[22],2.0)*t[23]*t[36]+0.5*t[18]*t[19]*t[23]*(t[35]+t[34]+8.0*t[1]*t[20])+t[19]*t[22]*t[24]*t[23];
215 d2fdrarb = 2.0*t[18]*t[27]*t[31]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t[24]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*(2.0*t[4]+t[34])+t[37]-0.25*t[18]*t[22]*t[23]*t[28]*t[36]+0.5*t[19]*t[23]*t[24]*t[28]-0.5*t[14]*t[18]*t[19]*t[22]*t[27]+0.5*t[19]*t[22]*t[23]*t[24]-.2757113292886222*t[32]*t[5]*t[14]*t[18];
216 d2fdraga= t[39]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[38]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]+1.2*t[6]*t[5]*t[14]*t[18]-0.45*t[10]*t[19]*t[22]*t[14]*t[18]+1.8*t[10]*t[5]*t[9]*t[31]*t[18];
217 d2fdragb = t[41]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[40]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]-0.45*t[11]*t[19]*t[22]*t[14]*t[18]+1.8*t[11]*t[5]*t[9]*t[31]*t[18];
218 d2fdraab = t[42]+0.030005*t[29]*t[12]*t[19]*t[22]*t[23]-0.06001*t[29]*t[12]*t[5]*t[9]*t[14];
219 d2fdrbrb = -1.0*t[14]*t[18]*t[5]*(2.8*gradb/pow(rhob,3.333333333333333)+t[33])+2.0*t[18]*pow(t[27],2.0)*t[31]*t[5]-2.0*t[14]*t[24]*t[27]*t[5]+t[37]-0.25*t[18]*t[23]*pow(t[28],2.0)*t[36]+0.5*t[18]*t[19]*t[23]*(t[35]+t[34]-8.0*t[1]*t[20])-1.0*t[14]*t[18]*t[19]*t[27]*t[28]+t[19]*t[28]*t[24]*t[23];
220 d2fdrbga = t[39]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[38]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14]-0.45*t[10]*t[19]*t[28]*t[14]*t[18]+1.8*t[10]*t[5]*t[27]*t[31]*t[18];
221 d2fdrbgb = t[41]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[40]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14]+1.2*t[26]*t[5]*t[14]*t[18]-0.45*t[11]*t[19]*t[28]*t[14]*t[18]+1.8*t[11]*t[5]*t[27]*t[31]*t[18];
222 d2fdrbab = t[42]+0.030005*t[29]*t[12]*t[19]*t[28]*t[23]-0.06001*t[29]*t[12]*t[5]*t[27]*t[14];
223 d2fdgaga = 1.62*t[18]*t[31]*t[5]/pow(rhoa,2.666666666666667)+t[30]-0.06001*t[15]*t[43]*t[12]*t[5]*t[23]-0.108018*grada*t[29]*t[10]*t[12]*t[5]*t[14];
224 d2fdgagb = -0.06001*grada*gradb*t[43]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10]*t[12]*t[5]*t[14]+1.62*t[10]*t[11]*t[5]*t[31]*t[18];
225 d2fdgaab = -0.06001*grada*t[43]*t[12]*t[5]*t[23]-0.054009*t[29]*t[10]*t[12]*t[5]*t[14];
226 d2fdgbgb = 1.62*t[18]*t[31]*t[5]/pow(rhob,2.666666666666667)+t[30]-0.06001*t[16]*t[43]*t[12]*t[5]*t[23]-0.108018*gradb*t[29]*t[11]*t[12]*t[5]*t[14];
227 d2fdgbab = -0.06001*gradb*t[43]*t[12]*t[5]*t[23]-0.054009*t[29]*t[11]*t[12]*t[5]*t[14];
228 d2fdabab = -0.06001*t[43]*t[12]*t[5]*t[23];
229 ds->df1000 += factor*dfdra;
230 ds->df0100 += factor*dfdrb;
231 ds->df0010 += factor*dfdga;
232 ds->df0001 += factor*dfdgb;
233 ds->df00001 += factor*dfdab;
234 ds->df2000 += factor*d2fdrara;
235 ds->df1100 += factor*d2fdrarb;
236 ds->df1010 += factor*d2fdraga;
237 ds->df1001 += factor*d2fdragb;
238 ds->df10001 += factor*d2fdraab;
239 ds->df0200 += factor*d2fdrbrb;
240 ds->df0110 += factor*d2fdrbga;
241 ds->df0101 += factor*d2fdrbgb;
242 ds->df01001 += factor*d2fdrbab;
243 ds->df0020 += factor*d2fdgaga;
244 ds->df0011 += factor*d2fdgagb;
245 ds->df00101+= factor*d2fdgaab;
246 ds->df0002 += factor*d2fdgbgb;
247 ds->df00011+= factor*d2fdgbab;
248 ds->df00002+= factor*d2fdabab;
249 };
250 }
251
252 static void
wl90c_third(FunThirdFuncDrv * ds,real factor,const FunDensProp * dp)253 wl90c_third(FunThirdFuncDrv *ds, real factor, const FunDensProp* dp)
254 {
255 real t[104];
256 real dfdra, dfdrb, dfdga, dfdgb, dfdab;
257 real d2fdrara , d2fdrarb , d2fdraga , d2fdragb ;
258 real d2fdraab , d2fdrbrb , d2fdrbga , d2fdrbgb ;
259 real d2fdrbab , d2fdgaga , d2fdgagb , d2fdgaab ;
260 real d2fdgbgb , d2fdgbab , d2fdabab ;
261 real d3fdrarara , d3fdrararb , d3fdraraga , d3fdraragb ;
262 real d3fdraraab , d3fdrarbrb , d3fdrarbga , d3fdrarbgb ;
263 real d3fdrarbab , d3fdragaga , d3fdragagb , d3fdragaab ;
264 real d3fdragbgb , d3fdragbab , d3fdraabab , d3fdrbrbrb ;
265 real d3fdrbrbga , d3fdrbrbgb , d3fdrbrbab , d3fdrbgaga ;
266 real d3fdrbgagb , d3fdrbgaab , d3fdrbgbgb , d3fdrbgbab ;
267 real d3fdrbabab , d3fdgagaga , d3fdgagagb , d3fdgagaab ;
268 real d3fdgagbgb , d3fdgagbab , d3fdgaabab , d3fdgbgbgb ;
269 real d3fdgbgbab , d3fdgbabab , d3fdababab ;
270 real rhoa = dp->rhoa;
271 real rhob = dp->rhob;
272 real grada = dp->grada;
273 real gradb = dp->gradb;
274 real gradab = dp->gradab;
275
276 if (rhoa + rhob > WL90C_THRESHOLD) {
277
278 t[1] = rhoa-1.0*rhob;
279 t[2] = pow(t[1],2.0);
280 t[3] = rhob+rhoa;
281 t[4] = 1/pow(t[3],2.0);
282 t[5] = sqrt(1.0-1.0*t[2]*t[4]);
283 t[6] = 1/pow(rhoa,2.333333333333333);
284 t[7] = 1/pow(t[3],1.333333333333333);
285 t[8] = -.2067834969664666*t[7];
286 t[9] = t[8]-1.2*grada*t[6];
287 t[10] = 1/pow(rhoa,1.333333333333333);
288 t[11] = 1/pow(rhob,1.333333333333333);
289 t[12] = 1/pow(t[3],.3333333333333333);
290 t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;
291 t[14] = 1/pow(t[13],2.0);
292 t[15] = pow(grada,2.0);
293 t[16] = pow(gradb,2.0);
294 t[17] = sqrt(2.0*gradab+t[16]+t[15]);
295 t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];
296 t[19] = 1/t[5];
297 t[20] = 1/pow(t[3],3.0);
298 t[21] = 2.0*t[2]*t[20];
299 t[22] = t[21]-2.0*t[1]*t[4];
300 t[23] = 1/t[13];
301 t[24] = -.02000333333333333*t[17]*t[7]-0.7486;
302 t[25] = t[5]*t[24]*t[23];
303 t[26] = 1/pow(rhob,2.333333333333333);
304 t[27] = t[8]-1.2*gradb*t[26];
305 t[28] = 2.0*t[1]*t[4]+t[21];
306 t[29] = 1/t[17];
307 t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];
308 t[31] = pow(t[9],2.0);
309 t[32] = 1/pow(t[13],3.0);
310 t[33] = 1/pow(rhoa,3.333333333333333);
311 t[34] = 1/pow(t[3],2.333333333333333);
312 t[35] = .2757113292886222*t[34];
313 t[36] = t[35]+2.8*grada*t[33];
314 t[37] = 1/pow(t[3],4.0);
315 t[38] = -6.0*t[2]*t[37];
316 t[39] = -2.0*t[4];
317 t[40] = t[39]+t[38]+8.0*t[1]*t[20];
318 t[41] = 1/pow(t[5],3.0);
319 t[42] = pow(t[22],2.0);
320 t[43] = .02667111111111111*t[17]*t[34]*t[5]*t[23];
321 t[44] = 2.0*t[4]+t[38];
322 t[45] = -0.9*t[10]*t[5]*t[24]*t[14];
323 t[46] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];
324 t[47] = -0.9*t[11]*t[5]*t[24]*t[14];
325 t[48] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];
326 t[49] = -0.06001*t[29]*t[12]*t[5]*t[9]*t[14];
327 t[50] = 0.030005*t[29]*t[12]*t[19]*t[22]*t[23];
328 t[51] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];
329 t[52] = pow(t[27],2.0);
330 t[53] = 1/pow(rhob,3.333333333333333);
331 t[54] = t[35]+2.8*gradb*t[53];
332 t[55] = t[39]+t[38]-8.0*t[1]*t[20];
333 t[56] = pow(t[28],2.0);
334 t[57] = -0.06001*t[29]*t[12]*t[5]*t[27]*t[14];
335 t[58] = 0.030005*t[29]*t[12]*t[19]*t[28]*t[23];
336 t[59] = 1/pow(rhoa,2.666666666666667);
337 t[60] = 1/pow(t[17],3.0);
338 t[61] = -0.054009*t[29]*t[10]*t[12]*t[5]*t[14];
339 t[62] = -0.06001*grada*t[60]*t[12]*t[5]*t[23];
340 t[63] = 1/pow(rhob,2.666666666666667);
341 t[64] = -0.054009*t[29]*t[11]*t[12]*t[5]*t[14];
342 t[65] = -0.06001*gradb*t[60]*t[12]*t[5]*t[23];
343 t[66] = -0.06001*t[60]*t[12]*t[5]*t[23];
344 t[67] = 1/pow(t[13],4.0);
345 t[68] = 1/pow(t[3],3.333333333333333);
346 t[69] = -.6433264350067851*t[68];
347 t[70] = 24.0*t[2]/pow(t[3],5.0);
348 t[71] = 12.0*t[20];
349 t[72] = 1/pow(t[5],5.0);
350 t[73] = -.06223259259259259*t[17]*t[68]*t[5]*t[23];
351 t[74] = .6433264350067851*t[68]*t[5]*t[14]*t[18];
352 t[75] = -4.0*t[20];
353 t[76] = 4.0*t[24]*t[27]*t[32]*t[5]*t[9];
354 t[77] = -1.0*t[14]*t[19]*t[24]*t[28]*t[9];
355 t[78] = -1.0*t[14]*t[19]*t[22]*t[24]*t[27];
356 t[79] = -.5514226585772444*t[34]*t[5]*t[24]*t[14];
357 t[80] = -0.5*t[22]*t[23]*t[24]*t[28]*t[41];
358 t[81] = t[44]*t[19]*t[24]*t[23];
359 t[82] = -0.024004*t[17]*t[10]*t[34]*t[5]*t[14];
360 t[83] = .02667111111111111*grada*t[29]*t[34]*t[5]*t[23];
361 t[84] = -0.024004*t[17]*t[11]*t[34]*t[5]*t[14];
362 t[85] = .02667111111111111*gradb*t[29]*t[34]*t[5]*t[23];
363 t[86] = .02667111111111111*t[29]*t[34]*t[5]*t[23];
364 t[87] = 1/pow(t[3],2.666666666666667);
365 t[88] = 1.62*t[59]*t[5]*t[24]*t[32];
366 t[89] = 0.036006*grada*t[29]*t[10]*t[7]*t[5]*t[14];
367 t[90] = .02000333333333333*t[15]*t[60]*t[7]*t[5]*t[23];
368 t[91] = 1.62*t[10]*t[11]*t[5]*t[24]*t[32];
369 t[92] = 0.018003*gradb*t[29]*t[10]*t[7]*t[5]*t[14];
370 t[93] = 0.018003*grada*t[29]*t[11]*t[7]*t[5]*t[14];
371 t[94] = .02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[23];
372 t[95] = 0.018003*t[29]*t[10]*t[7]*t[5]*t[14];
373 t[96] = .02000333333333333*grada*t[60]*t[7]*t[5]*t[23];
374 t[97] = 1.62*t[63]*t[5]*t[24]*t[32];
375 t[98] = 0.036006*gradb*t[29]*t[11]*t[7]*t[5]*t[14];
376 t[99] = .02000333333333333*t[16]*t[60]*t[7]*t[5]*t[23];
377 t[100] = 0.018003*t[29]*t[11]*t[7]*t[5]*t[14];
378 t[101] = .02000333333333333*gradb*t[60]*t[7]*t[5]*t[23];
379 t[102] = .02000333333333333*t[60]*t[7]*t[5]*t[23];
380 t[103] = 1/pow(t[17],5.0);
381 dfdra = -1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];
382 dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];
383 dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];
384 dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];
385 dfdab = t[30];
386 d2fdrara =-2.0*t[14]*t[24]*t[5]*t[9]-1.0*t[14]*t[18]*t[19]*t[22]*t[9]-1.0*t[14]*t[18]*t[36]*t[5]+2.0*t[18]*t[31]*t[32]*t[5]+t[43]-0.25*t[18]*t[23]*t[41]*t[42]+0.5*t[18]*t[19]*t[23]*t[40]+t[19]*t[22]*t[24]*t[23];
387 d2fdrarb = 2.0*t[18]*t[27]*t[32]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t[24]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[44]+t[43]-0.25*t[18]*t[22]*t[23]*t[28]*t[41]+0.5*t[19]*t[23]*t[24]*t[28]-0.5*t[14]*t[18]*t[19]*t[22]*t[27]+0.5*t[19]*t[22]*t[23]*t[24]-.2757113292886222*t[34]*t[5]*t[14]*t[18];
388 d2fdraga = t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]+1.2*t[6]*t[5]*t[14]*t[18]-0.45*t[10]*t[19]*t[22]*t[14]*t[18]+1.8*t[10]*t[5]*t[9]*t[32]*t[18];
389 d2fdragb = t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]-0.45*t[11]*t[19]*t[22]*t[14]*t[18]+1.8*t[11]*t[5]*t[9]*t[32]*t[18];
390 d2fdraab = t[51]+t[50]+t[49];
391 d2fdrbrb = -0.25*t[18]*t[23]*t[41]*t[56]+0.5*t[18]*t[19]*t[23]*t[55]-1.0*t[14]*t[18]*t[5]*t[54]+2.0*t[18]*t[32]*t[5]*t[52]-2.0*t[14]*t[24]*t[27]*t[5]+t[43]-1.0*t[14]*t[18]*t[19]*t[27]*t[28]+t[19]*t[28]*t[24]*t[23];
392 d2fdrbga = t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14]-0.45*t[10]*t[19]*t[28]*t[14]*t[18]+1.8*t[10]*t[5]*t[27]*t[32]*t[18];
393 d2fdrbgb = t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14]+1.2*t[26]*t[5]*t[14]*t[18]-0.45*t[11]*t[19]*t[28]*t[14]*t[18]+1.8*t[11]*t[5]*t[27]*t[32]*t[18];
394 d2fdrbab = t[51]+t[58]+t[57];
395 d2fdgaga = t[30]-0.06001*t[15]*t[60]*t[12]*t[5]*t[23]-0.108018*grada*t[29]*t[10]*t[12]*t[5]*t[14]+1.62*t[59]*t[5]*t[32]*t[18];
396 d2fdgagb = -0.06001*grada*gradb*t[60]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10]*t[12]*t[5]*t[14]+1.62*t[10]*t[11]*t[5]*t[32]*t[18];
397 d2fdgaab = t[62]+t[61];
398 d2fdgbgb = t[30]-0.06001*t[16]*t[60]*t[12]*t[5]*t[23]-0.108018*gradb*t[29]*t[11]*t[12]*t[5]*t[14]+1.62*t[63]*t[5]*t[32]*t[18];
399 d2fdgbab = t[65]+t[64];
400 d2fdabab = t[66];
401 d3fdrarara = -1.0*t[14]*t[18]*t[5]*(t[69]-9.333333333333334*grada/pow(rhoa,4.333333333333333))-6.0*t[18]*t[5]*t[67]*pow(t[9],3.0)+6.0*t[18]*t[32]*t[36]*t[5]*t[9]+0.75*t[14]*t[18]*t[41]*t[42]*t[9]-1.5*t[14]*t[18]*t[19]*t[40]*t[9]-3.0*t[14]*t[19]*t[22]*t[24]*t[9]+t[73]+0.375*t[18]*pow(t[22],3.0)*t[23]*t[72]+0.5*t[18]*t[19]*t[23]*(t[71]+t[70]-36.0*t[1]*t[37])-3.0*t[14]*t[24]*t[36]*t[5]+6.0*t[24]*t[31]*t[32]*t[5]-0.75*t[23]*t[24]*t[41]*t[42]-0.75*t[18]*t[22]*t[23]*t[40]*t[41]+1.5*t[19]*t[23]*t[24]*t[40]-1.5*t[14]*t[18]*t[19]*t[22]*t[36]+3.0*t[18]*t[19]*t[22]*t[31]*t[32]+.04000666666666666*t[17]*t[34]*t[19]*t[22]*t[23]-.08001333333333333*t[17]*t[34]*t[5]*t[9]*t[14];
402 d3fdrararb = -1.0*t[14]*t[18]*t[19]*t[44]*t[9]+0.5*t[14]*t[18]*t[22]*t[28]*t[41]*t[9]+2.0*t[18]*t[19]*t[22]*t[27]*t[32]*t[9]-1.0*t[14]*t[19]*t[22]*t[24]*t[9]+t[81]+t[80]+t[79]+t[78]+t[77]+t[76]+0.5*t[18]*t[19]*t[23]*(t[75]+t[70]-12.0*t[1]*t[37])+t[74]+t[73]+0.375*t[18]*t[23]*t[28]*t[42]*t[72]-6.0*t[18]*t[27]*t[31]*t[5]*t[67]+2.0*t[18]*t[27]*t[32]*t[36]*t[5]-1.0*t[14]*t[24]*t[36]*t[5]+2.0*t[24]*t[31]*t[32]*t[5]-0.5*t[18]*t[22]*t[23]*t[41]*t[44]+0.25*t[14]*t[18]*t[27]*t[41]*t[42]-0.25*t[23]*t[24]*t[41]*t[42]-0.25*t[18]*t[23]*t[28]*t[40]*t[41]-0.5*t[14]*t[18]*t[19]*t[27]*t[40]+0.5*t[19]*t[23]*t[24]*t[40]-0.5*t[14]*t[18]*t[19]*t[28]*t[36]+.01333555555555555*t[17]*t[34]*t[19]*t[28]*t[23]+.02667111111111111*t[17]*t[34]*t[19]*t[22]*t[23]+1.102845317154489*t[34]*t[5]*t[9]*t[32]*t[18]+t[19]*t[28]*t[31]*t[32]*t[18]-.2757113292886222*t[34]*t[19]*t[22]*t[14]*t[18]-.05334222222222221*t[17]*t[34]*t[5]*t[9]*t[14]-.02667111111111111*t[17]*t[34]*t[5]*t[27]*t[14];
403 d3fdraraga = -0.0150025*grada*t[29]*t[12]*t[41]*t[42]*t[23]+t[83]+0.030005*grada*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*grada*t[29]*t[7]*t[19]*t[22]*t[23]+t[82]-0.06001*grada*t[29]*t[12]*t[36]*t[5]*t[14]+2.4*t[6]*t[5]*t[24]*t[14]-0.9*t[10]*t[19]*t[22]*t[24]*t[14]+.04000666666666666*grada*t[29]*t[7]*t[5]*t[9]*t[14]-0.06001*grada*t[29]*t[12]*t[19]*t[22]*t[9]*t[14]+0.12002*grada*t[29]*t[12]*t[5]*t[31]*t[32]+3.6*t[10]*t[5]*t[9]*t[24]*t[32]+0.225*t[10]*t[41]*t[42]*t[14]*t[18]-2.8*t[33]*t[5]*t[14]*t[18]-0.45*t[10]*t[40]*t[19]*t[14]*t[18]+1.2*t[6]*t[19]*t[22]*t[14]*t[18]+1.8*t[10]*t[36]*t[5]*t[32]*t[18]-4.8*t[6]*t[5]*t[9]*t[32]*t[18]+1.8*t[10]*t[19]*t[22]*t[9]*t[32]*t[18]-5.4*t[10]*t[5]*t[31]*t[67]*t[18];
404 d3fdraragb = -0.0150025*gradb*t[29]*t[12]*t[41]*t[42]*t[23]+t[85]+0.030005*gradb*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*gradb*t[29]*t[7]*t[19]*t[22]*t[23]+t[84]-0.06001*gradb*t[29]*t[12]*t[36]*t[5]*t[14]-0.9*t[11]*t[19]*t[22]*t[24]*t[14]+.04000666666666666*gradb*t[29]*t[7]*t[5]*t[9]*t[14]-0.06001*gradb*t[29]*t[12]*t[19]*t[22]*t[9]*t[14]+0.12002*gradb*t[29]*t[12]*t[5]*t[31]*t[32]+3.6*t[11]*t[5]*t[9]*t[24]*t[32]+0.225*t[11]*t[41]*t[42]*t[14]*t[18]-0.45*t[11]*t[40]*t[19]*t[14]*t[18]+1.8*t[11]*t[36]*t[5]*t[32]*t[18]+1.8*t[11]*t[19]*t[22]*t[9]*t[32]*t[18]-5.4*t[11]*t[5]*t[31]*t[67]*t[18];
405 d3fdraraab= -0.0150025*t[29]*t[12]*t[41]*t[42]*t[23]+t[86]+0.030005*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*t[29]*t[7]*t[19]*t[22]*t[23]-0.06001*t[29]*t[12]*t[36]*t[5]*t[14]+.04000666666666666*t[29]*t[7]*t[5]*t[9]*t[14]-0.06001*t[29]*t[12]*t[19]*t[22]*t[9]*t[14]+0.12002*t[29]*t[12]*t[5]*t[31]*t[32];
406 d3fdrarbrb = -6.0*t[18]*t[5]*t[52]*t[67]*t[9]+0.25*t[14]*t[18]*t[41]*t[56]*t[9]-0.5*t[14]*t[18]*t[19]*t[55]*t[9]+2.0*t[18]*t[32]*t[5]*t[54]*t[9]+2.0*t[18]*t[19]*t[27]*t[28]*t[32]*t[9]+t[81]+t[80]+t[79]+t[78]+t[77]+t[76]+0.5*t[18]*t[19]*t[23]*(t[75]+t[70]+12.0*t[1]*t[37])+t[74]+t[73]+0.375*t[18]*t[22]*t[23]*t[56]*t[72]-0.25*t[23]*t[24]*t[41]*t[56]-0.25*t[18]*t[22]*t[23]*t[41]*t[55]+0.5*t[19]*t[23]*t[24]*t[55]-1.0*t[14]*t[24]*t[5]*t[54]-0.5*t[14]*t[18]*t[19]*t[22]*t[54]+2.0*t[24]*t[32]*t[5]*t[52]-0.5*t[18]*t[23]*t[28]*t[41]*t[44]-1.0*t[14]*t[18]*t[19]*t[27]*t[44]+0.5*t[14]*t[18]*t[22]*t[27]*t[28]*t[41]-1.0*t[14]*t[19]*t[24]*t[27]*t[28]+.02667111111111111*t[17]*t[34]*t[19]*t[28]*t[23]+.01333555555555555*t[17]*t[34]*t[19]*t[22]*t[23]+t[19]*t[22]*t[52]*t[32]*t[18]+1.102845317154489*t[34]*t[5]*t[27]*t[32]*t[18]-.2757113292886222*t[34]*t[19]*t[28]*t[14]*t[18]-.02667111111111111*t[17]*t[34]*t[5]*t[9]*t[14]-.05334222222222221*t[17]*t[34]*t[5]*t[27]*t[14];
407 d3fdrarbga = t[83]+0.030005*grada*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[28]*t[23]-0.0150025*grada*t[29]*t[12]*t[41]*t[22]*t[28]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[22]*t[23]+t[82]-.01654543687061022*grada*t[29]*t[87]*t[5]*t[14]+1.2*t[6]*t[5]*t[24]*t[14]-0.45*t[10]*t[19]*t[28]*t[24]*t[14]-0.45*t[10]*t[19]*t[22]*t[24]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[5]*t[27]*t[14]-0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[27]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[5]*t[9]*t[14]-0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[9]*t[14]+1.8*t[10]*t[5]*t[27]*t[24]*t[32]+1.8*t[10]*t[5]*t[9]*t[24]*t[32]+0.12002*grada*t[29]*t[12]*t[5]*t[9]*t[27]*t[32]-0.45*t[10]*t[44]*t[19]*t[14]*t[18]+0.6*t[6]*t[19]*t[28]*t[14]*t[18]+0.225*t[10]*t[41]*t[22]*t[28]*t[14]*t[18]+0.49628039271952*t[10]*t[34]*t[5]*t[32]*t[18]-2.4*t[6]*t[5]*t[27]*t[32]*t[18]+0.9*t[10]*t[19]*t[22]*t[27]*t[32]*t[18]+0.9*t[10]*t[19]*t[28]*t[9]*t[32]*t[18]-5.4*t[10]*t[5]*t[9]*t[27]*t[67]*t[18];
408 d3fdrarbgb = t[85]+0.030005*gradb*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[28]*t[23]-0.0150025*gradb*t[29]*t[12]*t[41]*t[22]*t[28]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[22]*t[23]+t[84]-.01654543687061022*gradb*t[29]*t[87]*t[5]*t[14]+1.2*t[26]*t[5]*t[24]*t[14]-0.45*t[11]*t[19]*t[28]*t[24]*t[14]-0.45*t[11]*t[19]*t[22]*t[24]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[5]*t[27]*t[14]-0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[27]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[5]*t[9]*t[14]-0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[9]*t[14]+1.8*t[11]*t[5]*t[27]*t[24]*t[32]+1.8*t[11]*t[5]*t[9]*t[24]*t[32]+0.12002*gradb*t[29]*t[12]*t[5]*t[9]*t[27]*t[32]-0.45*t[11]*t[44]*t[19]*t[14]*t[18]+0.225*t[11]*t[41]*t[22]*t[28]*t[14]*t[18]+0.6*t[26]*t[19]*t[22]*t[14]*t[18]+0.49628039271952*t[11]*t[34]*t[5]*t[32]*t[18]+0.9*t[11]*t[19]*t[22]*t[27]*t[32]*t[18]-2.4*t[26]*t[5]*t[9]*t[32]*t[18]+0.9*t[11]*t[19]*t[28]*t[9]*t[32]*t[18]-5.4*t[11]*t[5]*t[9]*t[27]*t[67]*t[18];
409 d3fdrarbab = t[86]+0.030005*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*t[29]*t[7]*t[19]*t[28]*t[23]-0.0150025*t[29]*t[12]*t[41]*t[22]*t[28]*t[23]-.01000166666666667*t[29]*t[7]*t[19]*t[22]*t[23]-.01654543687061022*t[29]*t[87]*t[5]*t[14]+.02000333333333333*t[29]*t[7]*t[5]*t[27]*t[14]-0.030005*t[29]*t[12]*t[19]*t[22]*t[27]*t[14]+.02000333333333333*t[29]*t[7]*t[5]*t[9]*t[14]-0.030005*t[29]*t[12]*t[19]*t[28]*t[9]*t[14]+0.12002*t[29]*t[12]*t[5]*t[9]*t[27]*t[32];
410 d3fdragaga = -4.32*t[18]*t[32]*t[5]/pow(rhoa,3.666666666666667)+t[90]+t[89]+t[88]+t[51]+t[50]+t[49]+0.216036*grada*t[29]*t[10]*t[12]*t[5]*t[9]*t[32]-0.030005*t[15]*t[60]*t[12]*t[19]*t[22]*t[23]-4.86*t[59]*t[5]*t[9]*t[67]*t[18]+0.81*t[59]*t[19]*t[22]*t[32]*t[18]+0.06001*t[15]*t[60]*t[12]*t[5]*t[9]*t[14]+0.144024*grada*t[29]*t[6]*t[12]*t[5]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[19]*t[22]*t[14];
411 d3fdragagb = t[94]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*gradb*t[29]*t[6]*t[12]*t[5]*t[14]+t[93]+t[92]+0.06001*grada*gradb*t[60]*t[12]*t[5]*t[9]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[19]*t[22]*t[14]+t[91]+0.108018*grada*t[29]*t[11]*t[12]*t[5]*t[9]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[5]*t[9]*t[32]-2.16*t[6]*t[11]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[19]*t[22]*t[32]*t[18]-4.86*t[10]*t[11]*t[5]*t[9]*t[67]*t[18];
412 d3fdragaab = t[96]-0.030005*grada*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*t[29]*t[6]*t[12]*t[5]*t[14]+t[95]+0.06001*grada*t[60]*t[12]*t[5]*t[9]*t[14]-0.0270045*t[29]*t[10]*t[12]*t[19]*t[22]*t[14]+0.108018*t[29]*t[10]*t[12]*t[5]*t[9]*t[32];
413 d3fdragbgb = t[51]+t[99]+t[50]-0.030005*t[16]*t[60]*t[12]*t[19]*t[22]*t[23]+t[98]+t[49]+0.06001*t[16]*t[60]*t[12]*t[5]*t[9]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]+t[97]+0.216036*gradb*t[29]*t[11]*t[12]*t[5]*t[9]*t[32]+0.81*t[63]*t[19]*t[22]*t[32]*t[18]-4.86*t[63]*t[5]*t[9]*t[67]*t[18];
414 d3fdragbab = t[101]-0.030005*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+t[100]+0.06001*gradb*t[60]*t[12]*t[5]*t[9]*t[14]-0.0270045*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]+0.108018*t[29]*t[11]*t[12]*t[5]*t[9]*t[32];
415 d3fdraabab = t[102]-0.030005*t[60]*t[12]*t[19]*t[22]*t[23]+0.06001*t[60]*t[12]*t[5]*t[9]*t[14];
416 d3fdrbrbrb =-1.0*t[14]*t[18]*t[5]*(t[69]-9.333333333333334*gradb/pow(rhob,4.333333333333333))+t[73]+0.375*t[18]*t[23]*pow(t[28],3.0)*t[72]+0.5*t[18]*t[19]*t[23]*(t[71]+t[70]+36.0*t[1]*t[37])-6.0*t[18]*pow(t[27],3.0)*t[5]*t[67]+0.75*t[14]*t[18]*t[27]*t[41]*t[56]-0.75*t[23]*t[24]*t[41]*t[56]-0.75*t[18]*t[23]*t[28]*t[41]*t[55]-1.5*t[14]*t[18]*t[19]*t[27]*t[55]+1.5*t[19]*t[23]*t[24]*t[55]+6.0*t[18]*t[27]*t[32]*t[5]*t[54]-3.0*t[14]*t[24]*t[5]*t[54]-1.5*t[14]*t[18]*t[19]*t[28]*t[54]+6.0*t[24]*t[32]*t[5]*t[52]+3.0*t[18]*t[19]*t[28]*t[32]*t[52]-3.0*t[14]*t[19]*t[24]*t[27]*t[28]+.04000666666666666*t[17]*t[34]*t[19]*t[28]*t[23]-.08001333333333333*t[17]*t[34]*t[5]*t[27]*t[14];
417 d3fdrbrbga = -0.0150025*grada*t[29]*t[12]*t[41]*t[56]*t[23]+t[83]+0.030005*grada*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*grada*t[29]*t[7]*t[19]*t[28]*t[23]+t[82]-0.06001*grada*t[29]*t[12]*t[54]*t[5]*t[14]-0.9*t[10]*t[19]*t[28]*t[24]*t[14]+.04000666666666666*grada*t[29]*t[7]*t[5]*t[27]*t[14]-0.06001*grada*t[29]*t[12]*t[19]*t[28]*t[27]*t[14]+0.12002*grada*t[29]*t[12]*t[5]*t[52]*t[32]+3.6*t[10]*t[5]*t[27]*t[24]*t[32]+0.225*t[10]*t[41]*t[56]*t[14]*t[18]-0.45*t[10]*t[55]*t[19]*t[14]*t[18]+1.8*t[10]*t[54]*t[5]*t[32]*t[18]+1.8*t[10]*t[19]*t[28]*t[27]*t[32]*t[18]-5.4*t[10]*t[5]*t[52]*t[67]*t[18];
418 d3fdrbrbgb = -0.0150025*gradb*t[29]*t[12]*t[41]*t[56]*t[23]+t[85]+0.030005*gradb*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*gradb*t[29]*t[7]*t[19]*t[28]*t[23]+t[84]-0.06001*gradb*t[29]*t[12]*t[54]*t[5]*t[14]+2.4*t[26]*t[5]*t[24]*t[14]-0.9*t[11]*t[19]*t[28]*t[24]*t[14]+.04000666666666666*gradb*t[29]*t[7]*t[5]*t[27]*t[14]-0.06001*gradb*t[29]*t[12]*t[19]*t[28]*t[27]*t[14]+0.12002*gradb*t[29]*t[12]*t[5]*t[52]*t[32]+3.6*t[11]*t[5]*t[27]*t[24]*t[32]+0.225*t[11]*t[41]*t[56]*t[14]*t[18]-2.8*t[53]*t[5]*t[14]*t[18]-0.45*t[11]*t[55]*t[19]*t[14]*t[18]+1.2*t[26]*t[19]*t[28]*t[14]*t[18]+1.8*t[11]*t[54]*t[5]*t[32]*t[18]-4.8*t[26]*t[5]*t[27]*t[32]*t[18]+1.8*t[11]*t[19]*t[28]*t[27]*t[32]*t[18]-5.4*t[11]*t[5]*t[52]*t[67]*t[18];
419 d3fdrbrbab = -0.0150025*t[29]*t[12]*t[41]*t[56]*t[23]+t[86]+0.030005*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*t[29]*t[7]*t[19]*t[28]*t[23]-0.06001*t[29]*t[12]*t[54]*t[5]*t[14]+.04000666666666666*t[29]*t[7]*t[5]*t[27]*t[14]-0.06001*t[29]*t[12]*t[19]*t[28]*t[27]*t[14]+0.12002*t[29]*t[12]*t[5]*t[52]*t[32];
420 d3fdrbgaga = t[51]+t[90]+t[58]-0.030005*t[15]*t[60]*t[12]*t[19]*t[28]*t[23]+t[89]+t[57]+0.06001*t[15]*t[60]*t[12]*t[5]*t[27]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+t[88]+0.216036*grada*t[29]*t[10]*t[12]*t[5]*t[27]*t[32]+0.81*t[59]*t[19]*t[28]*t[32]*t[18]-4.86*t[59]*t[5]*t[27]*t[67]*t[18];
421 d3fdrbgagb = t[94]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*grada*t[29]*t[26]*t[12]*t[5]*t[14]+t[93]+t[92]+0.06001*grada*gradb*t[60]*t[12]*t[5]*t[27]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+t[91]+0.108018*grada*t[29]*t[11]*t[12]*t[5]*t[27]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[5]*t[27]*t[32]-2.16*t[10]*t[26]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[19]*t[28]*t[32]*t[18]-4.86*t[10]*t[11]*t[5]*t[27]*t[67]*t[18];
422 d3fdrbgaab = t[96]-0.030005*grada*t[60]*t[12]*t[19]*t[28]*t[23]+t[95]+0.06001*grada*t[60]*t[12]*t[5]*t[27]*t[14]-0.0270045*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+0.108018*t[29]*t[10]*t[12]*t[5]*t[27]*t[32];
423 d3fdrbgbgb = -4.32*t[18]*t[32]*t[5]/pow(rhob,3.666666666666667)+t[99]+t[98]+t[97]+t[58]+t[57]+t[51]+0.216036*gradb*t[29]*t[11]*t[12]*t[5]*t[27]*t[32]-0.030005*t[16]*t[60]*t[12]*t[19]*t[28]*t[23]-4.86*t[63]*t[5]*t[27]*t[67]*t[18]+0.81*t[63]*t[19]*t[28]*t[32]*t[18]+0.144024*gradb*t[29]*t[26]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]+0.06001*t[16]*t[60]*t[12]*t[5]*t[27]*t[14];
424 d3fdrbgbab = t[101]-0.030005*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*t[29]*t[26]*t[12]*t[5]*t[14]+t[100]+0.06001*gradb*t[60]*t[12]*t[5]*t[27]*t[14]-0.0270045*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]+0.108018*t[29]*t[11]*t[12]*t[5]*t[27]*t[32];
425 d3fdrbabab = t[102]-0.030005*t[60]*t[12]*t[19]*t[28]*t[23]+0.06001*t[60]*t[12]*t[5]*t[27]*t[14];
426 d3fdgagaga = -4.374*t[18]*t[5]*t[67]/pow(rhoa,4.0)+0.18003*t[103]*t[12]*t[23]*t[5]*pow(grada,3.0)+0.2916486*grada*t[29]*t[59]*t[12]*t[5]*t[32]-0.18003*grada*t[60]*t[12]*t[5]*t[23]+0.162027*t[15]*t[60]*t[10]*t[12]*t[5]*t[14]-0.162027*t[29]*t[10]*t[12]*t[5]*t[14];
427 d3fdgagagb = t[65]+0.18003*t[15]*gradb*t[103]*t[12]*t[5]*t[23]+t[64]+0.054009*t[15]*t[60]*t[11]*t[12]*t[5]*t[14]+0.108018*grada*gradb*t[60]*t[10]*t[12]*t[5]*t[14]+0.1944324*grada*t[29]*t[10]*t[11]*t[12]*t[5]*t[32]+0.0972162*gradb*t[29]*t[59]*t[12]*t[5]*t[32]-4.374*t[59]*t[11]*t[5]*t[67]*t[18];
428 d3fdgagaab = t[66]+0.18003*t[15]*t[103]*t[12]*t[5]*t[23]+0.108018*grada*t[60]*t[10]*t[12]*t[5]*t[14]+0.0972162*t[29]*t[59]*t[12]*t[5]*t[32];
429 d3fdgagbgb = t[62]+0.18003*grada*t[16]*t[103]*t[12]*t[5]*t[23]+0.108018*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[14]+t[61]+0.054009*t[16]*t[60]*t[10]*t[12]*t[5]*t[14]+0.1944324*gradb*t[29]*t[10]*t[11]*t[12]*t[5]*t[32]+0.0972162*grada*t[29]*t[63]*t[12]*t[5]*t[32]-4.374*t[10]*t[63]*t[5]*t[67]*t[18];
430 d3fdgagbab = 0.18003*grada*gradb*t[103]*t[12]*t[5]*t[23]+0.054009*grada*t[60]*t[11]*t[12]*t[5]*t[14]+0.054009*gradb*t[60]*t[10]*t[12]*t[5]*t[14]+0.0972162*t[29]*t[10]*t[11]*t[12]*t[5]*t[32];
431 d3fdgaabab = 0.18003*grada*t[103]*t[12]*t[5]*t[23]+0.054009*t[60]*t[10]*t[12]*t[5]*t[14];
432 d3fdgbgbgb = -4.374*t[18]*t[5]*t[67]/pow(rhob,4.0)+0.18003*t[103]*t[12]*t[23]*t[5]*pow(gradb,3.0)+0.2916486*gradb*t[29]*t[63]*t[12]*t[5]*t[32]-0.18003*gradb*t[60]*t[12]*t[5]*t[23]+0.162027*t[16]*t[60]*t[11]*t[12]*t[5]*t[14]-0.162027*t[29]*t[11]*t[12]*t[5]*t[14];
433 d3fdgbgbab= t[66]+0.18003*t[16]*t[103]*t[12]*t[5]*t[23]+0.108018*gradb*t[60]*t[11]*t[12]*t[5]*t[14]+0.0972162*t[29]*t[63]*t[12]*t[5]*t[32];
434 d3fdgbabab = 0.18003*gradb*t[103]*t[12]*t[5]*t[23]+0.054009*t[60]*t[11]*t[12]*t[5]*t[14];
435 d3fdababab = 0.18003*t[103]*t[12]*t[5]*t[23];
436 ds->df1000 += factor*dfdra;
437 ds->df0100 += factor*dfdrb;
438 ds->df0010 += factor*dfdga;
439 ds->df0001 += factor*dfdgb;
440 ds->df00001 += factor*dfdab;
441 ds->df2000 += factor*d2fdrara;
442 ds->df1100 += factor*d2fdrarb;
443 ds->df1010 += factor*d2fdraga;
444 ds->df1001 += factor*d2fdragb;
445 ds->df10001 += factor*d2fdraab;
446 ds->df0200 += factor*d2fdrbrb;
447 ds->df0110 += factor*d2fdrbga;
448 ds->df0101 += factor*d2fdrbgb;
449 ds->df01001 += factor*d2fdrbab;
450 ds->df0020 += factor*d2fdgaga;
451 ds->df0011 += factor*d2fdgagb;
452 ds->df00101+= factor*d2fdgaab;
453 ds->df0002 += factor*d2fdgbgb;
454 ds->df00011+= factor*d2fdgbab;
455 ds->df00002+= factor*d2fdabab;
456 ds->df3000 += factor*d3fdrarara;
457 ds->df2100 += factor*d3fdrararb;
458 ds->df2010 += factor*d3fdraraga;
459 ds->df2001 += factor*d3fdraragb;
460 ds->df20001 += factor*d3fdraraab;
461 ds->df1200 += factor*d3fdrarbrb;
462 ds->df1110 += factor*d3fdrarbga;
463 ds->df1101 += factor*d3fdrarbgb;
464 ds->df11001 += factor*d3fdrarbab;
465 ds->df1020 += factor*d3fdragaga;
466 ds->df1011 += factor*d3fdragagb;
467 ds->df10101+= factor*d3fdragaab;
468 ds->df1002 += factor*d3fdragbgb;
469 ds->df10011+= factor*d3fdragbab;
470 ds->df10002+= factor*d3fdraabab;
471 ds->df0300 += factor*d3fdrbrbrb;
472 ds->df0210 += factor*d3fdrbrbga;
473 ds->df0201 += factor*d3fdrbrbgb;
474 ds->df02001 += factor*d3fdrbrbab;
475 ds->df0120 += factor*d3fdrbgaga;
476 ds->df0111 += factor*d3fdrbgagb;
477 ds->df01101+= factor*d3fdrbgaab;
478 ds->df0102 += factor*d3fdrbgbgb;
479 ds->df01011+= factor*d3fdrbgbab;
480 ds->df01002+= factor*d3fdrbgbab;
481 ds->df0030 += factor*d3fdgagaga;
482 ds->df0021 += factor*d3fdgagagb;
483 ds->df00201+= factor*d3fdgagaab;
484 ds->df0012 += factor*d3fdgagbgb;
485 ds->df00111+= factor*d3fdgagbab;
486 ds->df00102+= factor*d3fdgaabab;
487 ds->df0003 += factor*d3fdgbgbgb;
488 ds->df00021+= factor*d3fdgbgbab;
489 ds->df00012+= factor*d3fdgbabab;
490 ds->df00003+= factor*d3fdababab;
491 };
492 }
493
494 static void
wl90c_fourth(FunFourthFuncDrv * ds,real factor,const FunDensProp * dp)495 wl90c_fourth(FunFourthFuncDrv *ds, real factor, const FunDensProp* dp)
496 {
497 real t[300];
498 real dfdra, dfdrb, dfdga, dfdgb, dfdab;
499 real d2fdrara , d2fdrarb , d2fdraga , d2fdragb ;
500 real d2fdraab , d2fdrbrb , d2fdrbga , d2fdrbgb ;
501 real d2fdrbab , d2fdgaga , d2fdgagb , d2fdgaab ;
502 real d2fdgbgb , d2fdgbab , d2fdabab ;
503 real d3fdrarara , d3fdrararb , d3fdraraga , d3fdraragb ;
504 real d3fdraraab , d3fdrarbrb , d3fdrarbga , d3fdrarbgb ;
505 real d3fdrarbab , d3fdragaga , d3fdragagb , d3fdragaab ;
506 real d3fdragbgb , d3fdragbab , d3fdraabab , d3fdrbrbrb ;
507 real d3fdrbrbga , d3fdrbrbgb , d3fdrbrbab , d3fdrbgaga ;
508 real d3fdrbgagb , d3fdrbgaab , d3fdrbgbgb , d3fdrbgbab ;
509 real d3fdrbabab , d3fdgagaga , d3fdgagagb , d3fdgagaab ;
510 real d3fdgagbgb , d3fdgagbab , d3fdgaabab , d3fdgbgbgb ;
511 real d3fdgbgbab , d3fdgbabab , d3fdababab ;
512 real d4fdrararara , d4fdrarararb , d4fdrararaga , d4fdrararagb ;
513 real d4fdrararaab , d4fdrararbrb , d4fdrararbga , d4fdrararbgb ;
514 real d4fdrararbab , d4fdraragaga , d4fdraragagb , d4fdraragaab ;
515 real d4fdraragbgb , d4fdraragbab , d4fdraraabab , d4fdrarbrbrb ;
516 real d4fdrarbrbga , d4fdrarbrbgb , d4fdrarbrbab , d4fdrarbgaga ;
517 real d4fdrarbgagb , d4fdrarbgaab , d4fdrarbgbgb , d4fdrarbgbab ;
518 real d4fdrarbabab , d4fdragagaga , d4fdragagagb , d4fdragagaab ;
519 real d4fdragagbgb , d4fdragagbab , d4fdragaabab , d4fdragbgbgb ;
520 real d4fdragbgbab , d4fdragbabab , d4fdraababab , d4fdrbrbrbrb ;
521 real d4fdrbrbrbga , d4fdrbrbrbgb , d4fdrbrbrbab , d4fdrbrbgaga ;
522 real d4fdrbrbgagb , d4fdrbrbgaab , d4fdrbrbgbgb , d4fdrbrbgbab ;
523 real d4fdrbrbabab , d4fdrbgagaga , d4fdrbgagagb , d4fdrbgagaab ;
524 real d4fdrbgagbgb , d4fdrbgagbab , d4fdrbgaabab , d4fdrbgbgbgb ;
525 real d4fdrbgbgbab , d4fdrbgbabab , d4fdrbababab , d4fdgagagaga ;
526 real d4fdgagagagb , d4fdgagagaab , d4fdgagagbgb , d4fdgagagbab ;
527 real d4fdgagaabab , d4fdgagbgbgb , d4fdgagbgbab , d4fdgagbabab ;
528 real d4fdgaababab , d4fdgbgbgbgb , d4fdgbgbgbab , d4fdgbgbabab ;
529 real d4fdgbababab , d4fdabababab ;
530 real rhoa = dp->rhoa;
531 real rhob = dp->rhob;
532 real grada = dp->grada;
533 real gradb = dp->gradb;
534 real gradab = dp->gradab;
535
536 if (rhoa + rhob > WL90C_THRESHOLD) {
537
538 t[1] = rhoa-1.0*rhob;
539 t[2] = pow(t[1],2.0);
540 t[3] = rhob+rhoa;
541 t[4] = 1/pow(t[3],2.0);
542 t[5] = sqrt(1.0-1.0*t[2]*t[4]);
543 t[6] = 1/pow(rhoa,2.333333333333333);
544 t[7] = 1/pow(t[3],1.333333333333333);
545 t[8] = -.2067834969664666*t[7];
546 t[9] = t[8]-1.2*grada*t[6];
547 t[10] = 1/pow(rhoa,1.333333333333333);
548 t[11] = 1/pow(rhob,1.333333333333333);
549 t[12] = 1/pow(t[3],.3333333333333333);
550 t[13] = 0.6203504908994*t[12]+0.9*(gradb*t[11]+grada*t[10])+3.60073;
551 t[14] = 1/pow(t[13],2.0);
552 t[15] = pow(grada,2.0);
553 t[16] = pow(gradb,2.0);
554 t[17] = sqrt(2.0*gradab+t[16]+t[15]);
555 t[18] = 0.06001*t[17]*t[12]-0.7486*t[3];
556 t[19] = 1/t[5];
557 t[20] = 1/pow(t[3],3.0);
558 t[21] = 2.0*t[2]*t[20];
559 t[22] = t[21]-2.0*t[1]*t[4];
560 t[23] = 1/t[13];
561 t[24] = -.02000333333333333*t[17]*t[7]-0.7486;
562 t[25] = t[5]*t[24]*t[23];
563 t[26] = 1/pow(rhob,2.333333333333333);
564 t[27] = t[8]-1.2*gradb*t[26];
565 t[28] = 2.0*t[1]*t[4]+t[21];
566 t[29] = 1/t[17];
567 t[30] = 0.06001*t[29]*t[12]*t[5]*t[23];
568 t[31] = pow(t[9],2.0);
569 t[32] = 1/pow(t[13],3.0);
570 t[33] = 1/pow(rhoa,3.333333333333333);
571 t[34] = 1/pow(t[3],2.333333333333333);
572 t[35] = .2757113292886222*t[34];
573 t[36] = t[35]+2.8*grada*t[33];
574 t[37] = 1/pow(t[3],4.0);
575 t[38] = -6.0*t[2]*t[37];
576 t[39] = -2.0*t[4];
577 t[40] = t[39]+t[38]+8.0*t[1]*t[20];
578 t[41] = 1/pow(t[5],3.0);
579 t[42] = pow(t[22],2.0);
580 t[43] = .02667111111111111*t[17]*t[34]*t[5]*t[23];
581 t[44] = 2.0*t[4]+t[38];
582 t[45] = -0.9*t[10]*t[5]*t[24]*t[14];
583 t[46] = -.02000333333333333*grada*t[29]*t[7]*t[5]*t[23];
584 t[47] = -0.9*t[11]*t[5]*t[24]*t[14];
585 t[48] = -.02000333333333333*gradb*t[29]*t[7]*t[5]*t[23];
586 t[49] = -0.06001*t[29]*t[12]*t[5]*t[9]*t[14];
587 t[50] = 0.030005*t[29]*t[12]*t[19]*t[22]*t[23];
588 t[51] = -.02000333333333333*t[29]*t[7]*t[5]*t[23];
589 t[52] = pow(t[27],2.0);
590 t[53] = 1/pow(rhob,3.333333333333333);
591 t[54] = t[35]+2.8*gradb*t[53];
592 t[55] = t[39]+t[38]-8.0*t[1]*t[20];
593 t[56] = pow(t[28],2.0);
594 t[57] = -0.06001*t[29]*t[12]*t[5]*t[27]*t[14];
595 t[58] = 0.030005*t[29]*t[12]*t[19]*t[28]*t[23];
596 t[59] = 1/pow(rhoa,2.666666666666667);
597 t[60] = 1/pow(t[17],3.0);
598 t[61] = -0.054009*t[29]*t[10]*t[12]*t[5]*t[14];
599 t[62] = -0.06001*grada*t[60]*t[12]*t[5]*t[23];
600 t[63] = 1/pow(rhob,2.666666666666667);
601 t[64] = -0.054009*t[29]*t[11]*t[12]*t[5]*t[14];
602 t[65] = -0.06001*gradb*t[60]*t[12]*t[5]*t[23];
603 t[66] = -0.06001*t[60]*t[12]*t[5]*t[23];
604 t[67] = pow(t[9],3.0);
605 t[68] = 1/pow(t[13],4.0);
606 t[69] = 1/pow(rhoa,4.333333333333333);
607 t[70] = 1/pow(t[3],3.333333333333333);
608 t[71] = -.6433264350067851*t[70];
609 t[72] = t[71]-9.333333333333334*grada*t[69];
610 t[73] = 1/pow(t[3],5.0);
611 t[74] = 24.0*t[2]*t[73];
612 t[75] = 12.0*t[20];
613 t[76] = t[75]+t[74]-36.0*t[1]*t[37];
614 t[77] = 1/pow(t[5],5.0);
615 t[78] = pow(t[22],3.0);
616 t[79] = -.06223259259259259*t[17]*t[70]*t[5]*t[23];
617 t[80] = .6433264350067851*t[70]*t[5]*t[14]*t[18];
618 t[81] = -4.0*t[20];
619 t[82] = t[81]+t[74]-12.0*t[1]*t[37];
620 t[83] = 4.0*t[24]*t[27]*t[32]*t[5]*t[9];
621 t[84] = -1.0*t[14]*t[19]*t[24]*t[28]*t[9];
622 t[85] = -1.0*t[14]*t[19]*t[22]*t[24]*t[27];
623 t[86] = -.5514226585772444*t[34]*t[5]*t[24]*t[14];
624 t[87] = -0.5*t[22]*t[23]*t[24]*t[28]*t[41];
625 t[88] = t[44]*t[19]*t[24]*t[23];
626 t[89] = -0.024004*t[17]*t[10]*t[34]*t[5]*t[14];
627 t[90] = .02667111111111111*grada*t[29]*t[34]*t[5]*t[23];
628 t[91] = -0.024004*t[17]*t[11]*t[34]*t[5]*t[14];
629 t[92] = .02667111111111111*gradb*t[29]*t[34]*t[5]*t[23];
630 t[93] = 0.12002*t[29]*t[12]*t[5]*t[31]*t[32];
631 t[94] = -0.06001*t[29]*t[12]*t[19]*t[22]*t[9]*t[14];
632 t[95] = .04000666666666666*t[29]*t[7]*t[5]*t[9]*t[14];
633 t[96] = -0.06001*t[29]*t[12]*t[36]*t[5]*t[14];
634 t[97] = -.02000333333333333*t[29]*t[7]*t[19]*t[22]*t[23];
635 t[98] = 0.030005*t[29]*t[12]*t[40]*t[19]*t[23];
636 t[99] = .02667111111111111*t[29]*t[34]*t[5]*t[23];
637 t[100] = -0.0150025*t[29]*t[12]*t[41]*t[42]*t[23];
638 t[101] = t[81]+t[74]+12.0*t[1]*t[37];
639 t[102] = 1/pow(t[3],2.666666666666667);
640 t[103] = 0.12002*t[29]*t[12]*t[5]*t[9]*t[27]*t[32];
641 t[104] = -0.030005*t[29]*t[12]*t[19]*t[28]*t[9]*t[14];
642 t[105] = .02000333333333333*t[29]*t[7]*t[5]*t[9]*t[14];
643 t[106] = -0.030005*t[29]*t[12]*t[19]*t[22]*t[27]*t[14];
644 t[107] = .02000333333333333*t[29]*t[7]*t[5]*t[27]*t[14];
645 t[108] = -.01654543687061022*t[29]*t[102]*t[5]*t[14];
646 t[109] = -.01000166666666667*t[29]*t[7]*t[19]*t[22]*t[23];
647 t[110] = -0.0150025*t[29]*t[12]*t[41]*t[22]*t[28]*t[23];
648 t[111] = -.01000166666666667*t[29]*t[7]*t[19]*t[28]*t[23];
649 t[112] = 0.030005*t[29]*t[12]*t[44]*t[19]*t[23];
650 t[113] = 1/pow(rhoa,3.666666666666667);
651 t[114] = 1.62*t[59]*t[5]*t[24]*t[32];
652 t[115] = 0.036006*grada*t[29]*t[10]*t[7]*t[5]*t[14];
653 t[116] = .02000333333333333*t[15]*t[60]*t[7]*t[5]*t[23];
654 t[117] = 1.62*t[10]*t[11]*t[5]*t[24]*t[32];
655 t[118] = 0.018003*gradb*t[29]*t[10]*t[7]*t[5]*t[14];
656 t[119] = 0.018003*grada*t[29]*t[11]*t[7]*t[5]*t[14];
657 t[120] = .02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[23];
658 t[121] = 0.108018*t[29]*t[10]*t[12]*t[5]*t[9]*t[32];
659 t[122] = -0.0270045*t[29]*t[10]*t[12]*t[19]*t[22]*t[14];
660 t[123] = 0.06001*grada*t[60]*t[12]*t[5]*t[9]*t[14];
661 t[124] = 0.018003*t[29]*t[10]*t[7]*t[5]*t[14];
662 t[125] = .07201199999999999*t[29]*t[6]*t[12]*t[5]*t[14];
663 t[126] = -0.030005*grada*t[60]*t[12]*t[19]*t[22]*t[23];
664 t[127] = .02000333333333333*grada*t[60]*t[7]*t[5]*t[23];
665 t[128] = 1.62*t[63]*t[5]*t[24]*t[32];
666 t[129] = 0.036006*gradb*t[29]*t[11]*t[7]*t[5]*t[14];
667 t[130] = .02000333333333333*t[16]*t[60]*t[7]*t[5]*t[23];
668 t[131] = 0.108018*t[29]*t[11]*t[12]*t[5]*t[9]*t[32];
669 t[132] = -0.0270045*t[29]*t[11]*t[12]*t[19]*t[22]*t[14];
670 t[133] = 0.06001*gradb*t[60]*t[12]*t[5]*t[9]*t[14];
671 t[134] = 0.018003*t[29]*t[11]*t[7]*t[5]*t[14];
672 t[135] = -0.030005*gradb*t[60]*t[12]*t[19]*t[22]*t[23];
673 t[136] = .02000333333333333*gradb*t[60]*t[7]*t[5]*t[23];
674 t[137] = 0.06001*t[60]*t[12]*t[5]*t[9]*t[14];
675 t[138] = -0.030005*t[60]*t[12]*t[19]*t[22]*t[23];
676 t[139] = .02000333333333333*t[60]*t[7]*t[5]*t[23];
677 t[140] = pow(t[27],3.0);
678 t[141] = 1/pow(rhob,4.333333333333333);
679 t[142] = t[71]-9.333333333333334*gradb*t[141];
680 t[143] = t[75]+t[74]+36.0*t[1]*t[37];
681 t[144] = pow(t[28],3.0);
682 t[145] = 0.12002*t[29]*t[12]*t[5]*t[52]*t[32];
683 t[146] = -0.06001*t[29]*t[12]*t[19]*t[28]*t[27]*t[14];
684 t[147] = .04000666666666666*t[29]*t[7]*t[5]*t[27]*t[14];
685 t[148] = -0.06001*t[29]*t[12]*t[54]*t[5]*t[14];
686 t[149] = -.02000333333333333*t[29]*t[7]*t[19]*t[28]*t[23];
687 t[150] = 0.030005*t[29]*t[12]*t[55]*t[19]*t[23];
688 t[151] = -0.0150025*t[29]*t[12]*t[41]*t[56]*t[23];
689 t[152] = 0.108018*t[29]*t[10]*t[12]*t[5]*t[27]*t[32];
690 t[153] = -0.0270045*t[29]*t[10]*t[12]*t[19]*t[28]*t[14];
691 t[154] = 0.06001*grada*t[60]*t[12]*t[5]*t[27]*t[14];
692 t[155] = -0.030005*grada*t[60]*t[12]*t[19]*t[28]*t[23];
693 t[156] = 1/pow(rhob,3.666666666666667);
694 t[157] = 0.108018*t[29]*t[11]*t[12]*t[5]*t[27]*t[32];
695 t[158] = -0.0270045*t[29]*t[11]*t[12]*t[19]*t[28]*t[14];
696 t[159] = 0.06001*gradb*t[60]*t[12]*t[5]*t[27]*t[14];
697 t[160] = .07201199999999999*t[29]*t[26]*t[12]*t[5]*t[14];
698 t[161] = -0.030005*gradb*t[60]*t[12]*t[19]*t[28]*t[23];
699 t[162] = 0.06001*t[60]*t[12]*t[5]*t[27]*t[14];
700 t[163] = -0.030005*t[60]*t[12]*t[19]*t[28]*t[23];
701 t[164] = 1/pow(rhoa,4.0);
702 t[165] = pow(grada,3.0);
703 t[166] = 1/pow(t[17],5.0);
704 t[167] = 0.0972162*t[29]*t[59]*t[12]*t[5]*t[32];
705 t[168] = 0.108018*grada*t[60]*t[10]*t[12]*t[5]*t[14];
706 t[169] = 0.18003*t[15]*t[166]*t[12]*t[5]*t[23];
707 t[170] = 0.054009*t[60]*t[10]*t[12]*t[5]*t[14];
708 t[171] = 0.18003*grada*t[166]*t[12]*t[5]*t[23];
709 t[172] = 1/pow(rhob,4.0);
710 t[173] = pow(gradb,3.0);
711 t[174] = 0.0972162*t[29]*t[63]*t[12]*t[5]*t[32];
712 t[175] = 0.108018*gradb*t[60]*t[11]*t[12]*t[5]*t[14];
713 t[176] = 0.18003*t[16]*t[166]*t[12]*t[5]*t[23];
714 t[177] = 0.054009*t[60]*t[11]*t[12]*t[5]*t[14];
715 t[178] = 0.18003*gradb*t[166]*t[12]*t[5]*t[23];
716 t[179] = 0.18003*t[166]*t[12]*t[5]*t[23];
717 t[180] = 1/pow(t[13],5.0);
718 t[181] = 1/pow(rhoa,5.333333333333333);
719 t[182] = 1/pow(t[3],4.333333333333333);
720 t[183] = 2.144421450022617*t[182];
721 t[184] = -120.0*t[2]/pow(t[3],6.0);
722 t[185] = -72.0*t[37];
723 t[186] = 1/pow(t[5],7.0);
724 t[187] = 0.207441975308642*t[17]*t[182]*t[5]*t[23];
725 t[188] = -2.144421450022617*t[182]*t[5]*t[14]*t[18];
726 t[189] = .1600266666666667*t[17]*t[34]*t[5]*t[9]*t[27]*t[32];
727 t[190] = -.04000666666666666*t[17]*t[34]*t[19]*t[28]*t[9]*t[14];
728 t[191] = -.04000666666666666*t[17]*t[34]*t[19]*t[22]*t[27]*t[14];
729 t[192] = 1.929979305020355*t[70]*t[5]*t[24]*t[14];
730 t[193] = 1/pow(t[3],4.666666666666667);
731 t[194] = -.02206058249414695*t[17]*t[193]*t[5]*t[14];
732 t[195] = -.02000333333333333*t[17]*t[34]*t[41]*t[22]*t[28]*t[23];
733 t[196] = .04000666666666666*t[17]*t[34]*t[44]*t[19]*t[23];
734 t[197] = .05600933333333333*t[17]*t[10]*t[70]*t[5]*t[14];
735 t[198] = -.06223259259259259*grada*t[29]*t[70]*t[5]*t[23];
736 t[199] = .05600933333333333*t[17]*t[11]*t[70]*t[5]*t[14];
737 t[200] = -.06223259259259259*gradb*t[29]*t[70]*t[5]*t[23];
738 t[201] = -.06223259259259259*t[29]*t[70]*t[5]*t[23];
739 t[202] = -1.157987583012213*t[10]*t[70]*t[5]*t[32]*t[18];
740 t[203] = -10.8*t[10]*t[5]*t[9]*t[27]*t[24]*t[68];
741 t[204] = -.08001333333333333*grada*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];
742 t[205] = 1.8*t[10]*t[19]*t[28]*t[9]*t[24]*t[32];
743 t[206] = 1.8*t[10]*t[19]*t[22]*t[27]*t[24]*t[32];
744 t[207] = -4.8*t[6]*t[5]*t[27]*t[24]*t[32];
745 t[208] = .9925607854390399*t[10]*t[34]*t[5]*t[24]*t[32];
746 t[209] = .02000333333333333*grada*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];
747 t[210] = .02000333333333333*grada*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];
748 t[211] = 0.45*t[10]*t[41]*t[22]*t[28]*t[24]*t[14];
749 t[212] = 1.2*t[6]*t[19]*t[28]*t[24]*t[14];
750 t[213] = -0.9*t[10]*t[44]*t[19]*t[24]*t[14];
751 t[214] = 1/pow(t[3],3.666666666666667);
752 t[215] = .04963631061183065*grada*t[29]*t[214]*t[5]*t[14];
753 t[216] = .01000166666666667*grada*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];
754 t[217] = -.02000333333333333*grada*t[29]*t[7]*t[44]*t[19]*t[23];
755 t[218] = -1.157987583012213*t[11]*t[70]*t[5]*t[32]*t[18];
756 t[219] = -10.8*t[11]*t[5]*t[9]*t[27]*t[24]*t[68];
757 t[220] = -.08001333333333333*gradb*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];
758 t[221] = 1.8*t[11]*t[19]*t[28]*t[9]*t[24]*t[32];
759 t[222] = -4.8*t[26]*t[5]*t[9]*t[24]*t[32];
760 t[223] = 1.8*t[11]*t[19]*t[22]*t[27]*t[24]*t[32];
761 t[224] = .9925607854390399*t[11]*t[34]*t[5]*t[24]*t[32];
762 t[225] = .02000333333333333*gradb*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];
763 t[226] = .02000333333333333*gradb*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];
764 t[227] = 1.2*t[26]*t[19]*t[22]*t[24]*t[14];
765 t[228] = 0.45*t[11]*t[41]*t[22]*t[28]*t[24]*t[14];
766 t[229] = -0.9*t[11]*t[44]*t[19]*t[24]*t[14];
767 t[230] = .04963631061183065*gradb*t[29]*t[214]*t[5]*t[14];
768 t[231] = .01000166666666667*gradb*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];
769 t[232] = -.02000333333333333*gradb*t[29]*t[7]*t[44]*t[19]*t[23];
770 t[233] = -.08001333333333333*t[29]*t[7]*t[5]*t[9]*t[27]*t[32];
771 t[234] = .02000333333333333*t[29]*t[7]*t[19]*t[28]*t[9]*t[14];
772 t[235] = .02000333333333333*t[29]*t[7]*t[19]*t[22]*t[27]*t[14];
773 t[236] = .04963631061183065*t[29]*t[214]*t[5]*t[14];
774 t[237] = .01000166666666667*t[29]*t[7]*t[41]*t[22]*t[28]*t[23];
775 t[238] = -.02000333333333333*t[29]*t[7]*t[44]*t[19]*t[23];
776 t[239] = .04320719999999999*t[17]*t[59]*t[34]*t[5]*t[32];
777 t[240] = -0.048008*grada*t[29]*t[10]*t[34]*t[5]*t[14];
778 t[241] = -.02667111111111111*t[15]*t[60]*t[34]*t[5]*t[23];
779 t[242] = .04320719999999999*t[17]*t[10]*t[11]*t[34]*t[5]*t[32];
780 t[243] = -0.024004*gradb*t[29]*t[10]*t[34]*t[5]*t[14];
781 t[244] = -0.024004*grada*t[29]*t[11]*t[34]*t[5]*t[14];
782 t[245] = -.02667111111111111*grada*gradb*t[60]*t[34]*t[5]*t[23];
783 t[246] = -0.024004*t[29]*t[10]*t[34]*t[5]*t[14];
784 t[247] = -.02667111111111111*grada*t[60]*t[34]*t[5]*t[23];
785 t[248] = .04320719999999999*t[17]*t[63]*t[34]*t[5]*t[32];
786 t[249] = -0.048008*gradb*t[29]*t[11]*t[34]*t[5]*t[14];
787 t[250] = -.02667111111111111*t[16]*t[60]*t[34]*t[5]*t[23];
788 t[251] = -0.024004*t[29]*t[11]*t[34]*t[5]*t[14];
789 t[252] = -.02667111111111111*gradb*t[60]*t[34]*t[5]*t[23];
790 t[253] = -.02667111111111111*t[60]*t[34]*t[5]*t[23];
791 t[254] = -4.374*t[164]*t[5]*t[24]*t[68];
792 t[255] = -0.0972162*grada*t[29]*t[59]*t[7]*t[5]*t[32];
793 t[256] = -.05400899999999999*t[15]*t[60]*t[10]*t[7]*t[5]*t[14];
794 t[257] = .05400899999999999*t[29]*t[10]*t[7]*t[5]*t[14];
795 t[258] = -.06000999999999999*t[165]*t[166]*t[7]*t[5]*t[23];
796 t[259] = .06000999999999999*grada*t[60]*t[7]*t[5]*t[23];
797 t[260] = -4.374*t[59]*t[11]*t[5]*t[24]*t[68];
798 t[261] = -0.0324054*gradb*t[29]*t[59]*t[7]*t[5]*t[32];
799 t[262] = -0.0648108*grada*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];
800 t[263] = -0.036006*grada*gradb*t[60]*t[10]*t[7]*t[5]*t[14];
801 t[264] = -0.018003*t[15]*t[60]*t[11]*t[7]*t[5]*t[14];
802 t[265] = -.06000999999999999*t[15]*gradb*t[166]*t[7]*t[5]*t[23];
803 t[266] = -0.0324054*t[29]*t[59]*t[7]*t[5]*t[32];
804 t[267] = -0.036006*grada*t[60]*t[10]*t[7]*t[5]*t[14];
805 t[268] = -.06000999999999999*t[15]*t[166]*t[7]*t[5]*t[23];
806 t[269] = -4.374*t[10]*t[63]*t[5]*t[24]*t[68];
807 t[270] = -0.0324054*grada*t[29]*t[63]*t[7]*t[5]*t[32];
808 t[271] = -0.0648108*gradb*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];
809 t[272] = -0.018003*t[16]*t[60]*t[10]*t[7]*t[5]*t[14];
810 t[273] = -0.036006*grada*gradb*t[60]*t[11]*t[7]*t[5]*t[14];
811 t[274] = -.06000999999999999*grada*t[16]*t[166]*t[7]*t[5]*t[23];
812 t[275] = -0.0324054*t[29]*t[10]*t[11]*t[7]*t[5]*t[32];
813 t[276] = -0.018003*gradb*t[60]*t[10]*t[7]*t[5]*t[14];
814 t[277] = -0.018003*grada*t[60]*t[11]*t[7]*t[5]*t[14];
815 t[278] = -.06000999999999999*grada*gradb*t[166]*t[7]*t[5]*t[23];
816 t[279] = -0.018003*t[60]*t[10]*t[7]*t[5]*t[14];
817 t[280] = -.06000999999999999*grada*t[166]*t[7]*t[5]*t[23];
818 t[281] = -4.374*t[172]*t[5]*t[24]*t[68];
819 t[282] = -0.0972162*gradb*t[29]*t[63]*t[7]*t[5]*t[32];
820 t[283] = -.05400899999999999*t[16]*t[60]*t[11]*t[7]*t[5]*t[14];
821 t[284] = .05400899999999999*t[29]*t[11]*t[7]*t[5]*t[14];
822 t[285] = -.06000999999999999*t[173]*t[166]*t[7]*t[5]*t[23];
823 t[286] = .06000999999999999*gradb*t[60]*t[7]*t[5]*t[23];
824 t[287] = -0.0324054*t[29]*t[63]*t[7]*t[5]*t[32];
825 t[288] = -0.036006*gradb*t[60]*t[11]*t[7]*t[5]*t[14];
826 t[289] = -.06000999999999999*t[16]*t[166]*t[7]*t[5]*t[23];
827 t[290] = -0.018003*t[60]*t[11]*t[7]*t[5]*t[14];
828 t[291] = -.06000999999999999*gradb*t[166]*t[7]*t[5]*t[23];
829 t[292] = -.06000999999999999*t[166]*t[7]*t[5]*t[23];
830 t[293] = 1/pow(rhob,5.333333333333333);
831 t[294] = 1/pow(t[17],7.0);
832 t[295] = -0.18003*t[60]*t[12]*t[5]*t[23];
833 t[296] = 0.2916486*t[29]*t[10]*t[11]*t[12]*t[5]*t[32];
834 t[297] = 0.162027*gradb*t[60]*t[10]*t[12]*t[5]*t[14];
835 t[298] = 0.162027*grada*t[60]*t[11]*t[12]*t[5]*t[14];
836 t[299] = 0.54009*grada*gradb*t[166]*t[12]*t[5]*t[23];
837 dfdra = -1.0*t[14]*t[18]*t[5]*t[9]+t[25]+0.5*t[18]*t[19]*t[22]*t[23];
838 dfdrb = -1.0*t[14]*t[18]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[28]+t[25];
839 dfdga = 0.06001*grada*t[29]*t[12]*t[5]*t[23]-0.9*t[10]*t[5]*t[14]*t[18];
840 dfdgb = 0.06001*gradb*t[29]*t[12]*t[5]*t[23]-0.9*t[11]*t[5]*t[14]*t[18];
841 dfdab = t[30];
842 d2fdrara = -2.0*t[14]*t[24]*t[5]*t[9]-1.0*t[14]*t[18]*t[19]*t[22]*t[9]-1.0*t[14]*t[18]*t[36]*t[5]+2.0*t[18]*t[31]*t[32]*t[5]+t[43]-0.25*t[18]*t[23]*t[41]*t[42]+0.5*t[18]*t[19]*t[23]*t[40]+t[19]*t[22]*t[24]*t[23];
843 d2fdrarb = 2.0*t[18]*t[27]*t[32]*t[5]*t[9]-1.0*t[14]*t[24]*t[5]*t[9]-0.5*t[14]*t[18]*t[19]*t[28]*t[9]-1.0*t[14]*t[24]*t[27]*t[5]+0.5*t[18]*t[19]*t[23]*t[44]+t[43]-0.25*t[18]*t[22]*t[23]*t[28]*t[41]+0.5*t[19]*t[23]*t[24]*t[28]-0.5*t[14]*t[18]*t[19]*t[22]*t[27]+0.5*t[19]*t[22]*t[23]*t[24]-.2757113292886222*t[34]*t[5]*t[14]*t[18];
844 d2fdraga = t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[9]*t[14]+1.2*t[6]*t[5]*t[14]*t[18]-0.45*t[10]*t[19]*t[22]*t[14]*t[18]+1.8*t[10]*t[5]*t[9]*t[32]*t[18];
845 d2fdragb = t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[9]*t[14]-0.45*t[11]*t[19]*t[22]*t[14]*t[18]+1.8*t[11]*t[5]*t[9]*t[32]*t[18];
846 d2fdraab = t[51]+t[50]+t[49];
847 d2fdrbrb = -0.25*t[18]*t[23]*t[41]*t[56]+0.5*t[18]*t[19]*t[23]*t[55]-1.0*t[14]*t[18]*t[5]*t[54]+2.0*t[18]*t[32]*t[5]*t[52]-2.0*t[14]*t[24]*t[27]*t[5]+t[43]-1.0*t[14]*t[18]*t[19]*t[27]*t[28]+t[19]*t[28]*t[24]*t[23];
848 d2fdrbga = t[46]+0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[23]+t[45]-0.06001*grada*t[29]*t[12]*t[5]*t[27]*t[14]-0.45*t[10]*t[19]*t[28]*t[14]*t[18]+1.8*t[10]*t[5]*t[27]*t[32]*t[18];
849 d2fdrbgb = t[48]+0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[23]+t[47]-0.06001*gradb*t[29]*t[12]*t[5]*t[27]*t[14]+1.2*t[26]*t[5]*t[14]*t[18]-0.45*t[11]*t[19]*t[28]*t[14]*t[18]+1.8*t[11]*t[5]*t[27]*t[32]*t[18];
850 d2fdrbab = t[51]+t[58]+t[57];
851 d2fdgaga = t[30]-0.06001*t[15]*t[60]*t[12]*t[5]*t[23]-0.108018*grada*t[29]*t[10]*t[12]*t[5]*t[14]+1.62*t[59]*t[5]*t[32]*t[18];
852 d2fdgagb = -0.06001*grada*gradb*t[60]*t[12]*t[5]*t[23]-0.054009*grada*t[29]*t[11]*t[12]*t[5]*t[14]-0.054009*gradb*t[29]*t[10]*t[12]*t[5]*t[14]+1.62*t[10]*t[11]*t[5]*t[32]*t[18];
853 d2fdgaab = t[62]+t[61];
854 d2fdgbgb = t[30]-0.06001*t[16]*t[60]*t[12]*t[5]*t[23]-0.108018*gradb*t[29]*t[11]*t[12]*t[5]*t[14]+1.62*t[63]*t[5]*t[32]*t[18];
855 d2fdgbab = t[65]+t[64];
856 d2fdabab = t[66];
857 d3fdrarara = 6.0*t[18]*t[32]*t[36]*t[5]*t[9]+0.75*t[14]*t[18]*t[41]*t[42]*t[9]-1.5*t[14]*t[18]*t[19]*t[40]*t[9]-3.0*t[14]*t[19]*t[22]*t[24]*t[9]+t[79]+0.375*t[18]*t[23]*t[77]*t[78]+0.5*t[18]*t[19]*t[23]*t[76]-1.0*t[14]*t[18]*t[5]*t[72]-6.0*t[18]*t[5]*t[67]*t[68]-3.0*t[14]*t[24]*t[36]*t[5]+6.0*t[24]*t[31]*t[32]*t[5]-0.75*t[23]*t[24]*t[41]*t[42]-0.75*t[18]*t[22]*t[23]*t[40]*t[41]+1.5*t[19]*t[23]*t[24]*t[40]-1.5*t[14]*t[18]*t[19]*t[22]*t[36]+3.0*t[18]*t[19]*t[22]*t[31]*t[32]+.04000666666666666*t[17]*t[34]*t[19]*t[22]*t[23]-.08001333333333333*t[17]*t[34]*t[5]*t[9]*t[14];
858 d3fdrararb = -1.0*t[14]*t[18]*t[19]*t[44]*t[9]+0.5*t[14]*t[18]*t[22]*t[28]*t[41]*t[9]+2.0*t[18]*t[19]*t[22]*t[27]*t[32]*t[9]-1.0*t[14]*t[19]*t[22]*t[24]*t[9]+t[88]+t[87]+t[86]+t[85]+t[84]+t[83]+0.5*t[18]*t[19]*t[23]*t[82]+t[80]+t[79]+0.375*t[18]*t[23]*t[28]*t[42]*t[77]-6.0*t[18]*t[27]*t[31]*t[5]*t[68]+2.0*t[18]*t[27]*t[32]*t[36]*t[5]-1.0*t[14]*t[24]*t[36]*t[5]+2.0*t[24]*t[31]*t[32]*t[5]-0.5*t[18]*t[22]*t[23]*t[41]*t[44]+0.25*t[14]*t[18]*t[27]*t[41]*t[42]-0.25*t[23]*t[24]*t[41]*t[42]-0.25*t[18]*t[23]*t[28]*t[40]*t[41]-0.5*t[14]*t[18]*t[19]*t[27]*t[40]+0.5*t[19]*t[23]*t[24]*t[40]-0.5*t[14]*t[18]*t[19]*t[28]*t[36]+.01333555555555555*t[17]*t[34]*t[19]*t[28]*t[23]+.02667111111111111*t[17]*t[34]*t[19]*t[22]*t[23]+1.102845317154489*t[34]*t[5]*t[9]*t[32]*t[18]+t[19]*t[28]*t[31]*t[32]*t[18]-.2757113292886222*t[34]*t[19]*t[22]*t[14]*t[18]-.05334222222222221*t[17]*t[34]*t[5]*t[9]*t[14]-.02667111111111111*t[17]*t[34]*t[5]*t[27]*t[14];
859 d3fdraraga = -0.0150025*grada*t[29]*t[12]*t[41]*t[42]*t[23]+t[90]+0.030005*grada*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*grada*t[29]*t[7]*t[19]*t[22]*t[23]+t[89]-0.06001*grada*t[29]*t[12]*t[36]*t[5]*t[14]+2.4*t[6]*t[5]*t[24]*t[14]-0.9*t[10]*t[19]*t[22]*t[24]*t[14]+.04000666666666666*grada*t[29]*t[7]*t[5]*t[9]*t[14]-0.06001*grada*t[29]*t[12]*t[19]*t[22]*t[9]*t[14]+0.12002*grada*t[29]*t[12]*t[5]*t[31]*t[32]+3.6*t[10]*t[5]*t[9]*t[24]*t[32]+0.225*t[10]*t[41]*t[42]*t[14]*t[18]-2.8*t[33]*t[5]*t[14]*t[18]-0.45*t[10]*t[40]*t[19]*t[14]*t[18]+1.2*t[6]*t[19]*t[22]*t[14]*t[18]+1.8*t[10]*t[36]*t[5]*t[32]*t[18]-4.8*t[6]*t[5]*t[9]*t[32]*t[18]+1.8*t[10]*t[19]*t[22]*t[9]*t[32]*t[18]-5.4*t[10]*t[5]*t[31]*t[68]*t[18];
860 d3fdraragb = -0.0150025*gradb*t[29]*t[12]*t[41]*t[42]*t[23]+t[92]+0.030005*gradb*t[29]*t[12]*t[40]*t[19]*t[23]-.02000333333333333*gradb*t[29]*t[7]*t[19]*t[22]*t[23]+t[91]-0.06001*gradb*t[29]*t[12]*t[36]*t[5]*t[14]-0.9*t[11]*t[19]*t[22]*t[24]*t[14]+.04000666666666666*gradb*t[29]*t[7]*t[5]*t[9]*t[14]-0.06001*gradb*t[29]*t[12]*t[19]*t[22]*t[9]*t[14]+0.12002*gradb*t[29]*t[12]*t[5]*t[31]*t[32]+3.6*t[11]*t[5]*t[9]*t[24]*t[32]+0.225*t[11]*t[41]*t[42]*t[14]*t[18]-0.45*t[11]*t[40]*t[19]*t[14]*t[18]+1.8*t[11]*t[36]*t[5]*t[32]*t[18]+1.8*t[11]*t[19]*t[22]*t[9]*t[32]*t[18]-5.4*t[11]*t[5]*t[31]*t[68]*t[18];
861 d3fdraraab = t[100]+t[99]+t[98]+t[97]+t[96]+t[95]+t[94]+t[93];
862 d3fdrarbrb = -6.0*t[18]*t[5]*t[52]*t[68]*t[9]+0.25*t[14]*t[18]*t[41]*t[56]*t[9]-0.5*t[14]*t[18]*t[19]*t[55]*t[9]+2.0*t[18]*t[32]*t[5]*t[54]*t[9]+2.0*t[18]*t[19]*t[27]*t[28]*t[32]*t[9]+t[88]+t[87]+t[86]+t[85]+t[84]+t[83]+t[80]+t[79]+0.375*t[18]*t[22]*t[23]*t[56]*t[77]-0.25*t[23]*t[24]*t[41]*t[56]-0.25*t[18]*t[22]*t[23]*t[41]*t[55]+0.5*t[19]*t[23]*t[24]*t[55]-1.0*t[14]*t[24]*t[5]*t[54]-0.5*t[14]*t[18]*t[19]*t[22]*t[54]+2.0*t[24]*t[32]*t[5]*t[52]-0.5*t[18]*t[23]*t[28]*t[41]*t[44]-1.0*t[14]*t[18]*t[19]*t[27]*t[44]+0.5*t[14]*t[18]*t[22]*t[27]*t[28]*t[41]-1.0*t[14]*t[19]*t[24]*t[27]*t[28]+.02667111111111111*t[17]*t[34]*t[19]*t[28]*t[23]+.01333555555555555*t[17]*t[34]*t[19]*t[22]*t[23]+0.5*t[101]*t[18]*t[19]*t[23]+t[19]*t[22]*t[52]*t[32]*t[18]+1.102845317154489*t[34]*t[5]*t[27]*t[32]*t[18]-.2757113292886222*t[34]*t[19]*t[28]*t[14]*t[18]-.02667111111111111*t[17]*t[34]*t[5]*t[9]*t[14]-.05334222222222221*t[17]*t[34]*t[5]*t[27]*t[14];
863 d3fdrarbga = t[90]+0.030005*grada*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[28]*t[23]-0.0150025*grada*t[29]*t[12]*t[41]*t[22]*t[28]*t[23]-.01000166666666667*grada*t[29]*t[7]*t[19]*t[22]*t[23]+t[89]-.01654543687061022*grada*t[29]*t[102]*t[5]*t[14]+1.2*t[6]*t[5]*t[24]*t[14]-0.45*t[10]*t[19]*t[28]*t[24]*t[14]-0.45*t[10]*t[19]*t[22]*t[24]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[5]*t[27]*t[14]-0.030005*grada*t[29]*t[12]*t[19]*t[22]*t[27]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[5]*t[9]*t[14]-0.030005*grada*t[29]*t[12]*t[19]*t[28]*t[9]*t[14]+1.8*t[10]*t[5]*t[27]*t[24]*t[32]+1.8*t[10]*t[5]*t[9]*t[24]*t[32]+0.12002*grada*t[29]*t[12]*t[5]*t[9]*t[27]*t[32]-0.45*t[10]*t[44]*t[19]*t[14]*t[18]+0.6*t[6]*t[19]*t[28]*t[14]*t[18]+0.225*t[10]*t[41]*t[22]*t[28]*t[14]*t[18]+0.49628039271952*t[10]*t[34]*t[5]*t[32]*t[18]-2.4*t[6]*t[5]*t[27]*t[32]*t[18]+0.9*t[10]*t[19]*t[22]*t[27]*t[32]*t[18]+0.9*t[10]*t[19]*t[28]*t[9]*t[32]*t[18]-5.4*t[10]*t[5]*t[9]*t[27]*t[68]*t[18];
864 d3fdrarbgb = t[92]+0.030005*gradb*t[29]*t[12]*t[44]*t[19]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[28]*t[23]-0.0150025*gradb*t[29]*t[12]*t[41]*t[22]*t[28]*t[23]-.01000166666666667*gradb*t[29]*t[7]*t[19]*t[22]*t[23]+t[91]-.01654543687061022*gradb*t[29]*t[102]*t[5]*t[14]+1.2*t[26]*t[5]*t[24]*t[14]-0.45*t[11]*t[19]*t[28]*t[24]*t[14]-0.45*t[11]*t[19]*t[22]*t[24]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[5]*t[27]*t[14]-0.030005*gradb*t[29]*t[12]*t[19]*t[22]*t[27]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[5]*t[9]*t[14]-0.030005*gradb*t[29]*t[12]*t[19]*t[28]*t[9]*t[14]+1.8*t[11]*t[5]*t[27]*t[24]*t[32]+1.8*t[11]*t[5]*t[9]*t[24]*t[32]+0.12002*gradb*t[29]*t[12]*t[5]*t[9]*t[27]*t[32]-0.45*t[11]*t[44]*t[19]*t[14]*t[18]+0.225*t[11]*t[41]*t[22]*t[28]*t[14]*t[18]+0.6*t[26]*t[19]*t[22]*t[14]*t[18]+0.49628039271952*t[11]*t[34]*t[5]*t[32]*t[18]+0.9*t[11]*t[19]*t[22]*t[27]*t[32]*t[18]-2.4*t[26]*t[5]*t[9]*t[32]*t[18]+0.9*t[11]*t[19]*t[28]*t[9]*t[32]*t[18]-5.4*t[11]*t[5]*t[9]*t[27]*t[68]*t[18];
865 d3fdrarbab = t[99]+t[112]+t[111]+t[110]+t[109]+t[108]+t[107]+t[106]+t[105]+t[104]+t[103];
866 d3fdragaga = t[51]+t[116]+t[50]-0.030005*t[15]*t[60]*t[12]*t[19]*t[22]*t[23]+0.144024*grada*t[29]*t[6]*t[12]*t[5]*t[14]+t[115]+t[49]+0.06001*t[15]*t[60]*t[12]*t[5]*t[9]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[19]*t[22]*t[14]+t[114]+0.216036*grada*t[29]*t[10]*t[12]*t[5]*t[9]*t[32]-4.32*t[113]*t[5]*t[32]*t[18]+0.81*t[59]*t[19]*t[22]*t[32]*t[18]-4.86*t[59]*t[5]*t[9]*t[68]*t[18];
867 d3fdragagb = t[120]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+.07201199999999999*gradb*t[29]*t[6]*t[12]*t[5]*t[14]+t[119]+t[118]+0.06001*grada*gradb*t[60]*t[12]*t[5]*t[9]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[19]*t[22]*t[14]+t[117]+0.108018*grada*t[29]*t[11]*t[12]*t[5]*t[9]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[5]*t[9]*t[32]-2.16*t[6]*t[11]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[19]*t[22]*t[32]*t[18]-4.86*t[10]*t[11]*t[5]*t[9]*t[68]*t[18];
868 d3fdragaab = t[127]+t[126]+t[125]+t[124]+t[123]+t[122]+t[121];
869 d3fdragbgb =t[51]+t[130]+t[50]-0.030005*t[16]*t[60]*t[12]*t[19]*t[22]*t[23]+t[129]+t[49]+0.06001*t[16]*t[60]*t[12]*t[5]*t[9]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]+t[128]+0.216036*gradb*t[29]*t[11]*t[12]*t[5]*t[9]*t[32]+0.81*t[63]*t[19]*t[22]*t[32]*t[18]-4.86*t[63]*t[5]*t[9]*t[68]*t[18];
870 d3fdragbab = t[136]+t[135]+t[134]+t[133]+t[132]+t[131];
871 d3fdraabab = t[139]+t[138]+t[137];
872 d3fdrbrbrb = t[79]+0.375*t[144]*t[18]*t[23]*t[77]-6.0*t[140]*t[18]*t[5]*t[68]+0.75*t[14]*t[18]*t[27]*t[41]*t[56]-0.75*t[23]*t[24]*t[41]*t[56]-0.75*t[18]*t[23]*t[28]*t[41]*t[55]-1.5*t[14]*t[18]*t[19]*t[27]*t[55]+1.5*t[19]*t[23]*t[24]*t[55]+6.0*t[18]*t[27]*t[32]*t[5]*t[54]-3.0*t[14]*t[24]*t[5]*t[54]-1.5*t[14]*t[18]*t[19]*t[28]*t[54]+6.0*t[24]*t[32]*t[5]*t[52]+3.0*t[18]*t[19]*t[28]*t[32]*t[52]-1.0*t[14]*t[142]*t[18]*t[5]-3.0*t[14]*t[19]*t[24]*t[27]*t[28]+.04000666666666666*t[17]*t[34]*t[19]*t[28]*t[23]+0.5*t[143]*t[18]*t[19]*t[23]-.08001333333333333*t[17]*t[34]*t[5]*t[27]*t[14];
873 d3fdrbrbga = -0.0150025*grada*t[29]*t[12]*t[41]*t[56]*t[23]+t[90]+0.030005*grada*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*grada*t[29]*t[7]*t[19]*t[28]*t[23]+t[89]-0.06001*grada*t[29]*t[12]*t[54]*t[5]*t[14]-0.9*t[10]*t[19]*t[28]*t[24]*t[14]+.04000666666666666*grada*t[29]*t[7]*t[5]*t[27]*t[14]-0.06001*grada*t[29]*t[12]*t[19]*t[28]*t[27]*t[14]+0.12002*grada*t[29]*t[12]*t[5]*t[52]*t[32]+3.6*t[10]*t[5]*t[27]*t[24]*t[32]+0.225*t[10]*t[41]*t[56]*t[14]*t[18]-0.45*t[10]*t[55]*t[19]*t[14]*t[18]+1.8*t[10]*t[54]*t[5]*t[32]*t[18]+1.8*t[10]*t[19]*t[28]*t[27]*t[32]*t[18]-5.4*t[10]*t[5]*t[52]*t[68]*t[18];
874 d3fdrbrbgb = -0.0150025*gradb*t[29]*t[12]*t[41]*t[56]*t[23]+t[92]+0.030005*gradb*t[29]*t[12]*t[55]*t[19]*t[23]-.02000333333333333*gradb*t[29]*t[7]*t[19]*t[28]*t[23]+t[91]-0.06001*gradb*t[29]*t[12]*t[54]*t[5]*t[14]+2.4*t[26]*t[5]*t[24]*t[14]-0.9*t[11]*t[19]*t[28]*t[24]*t[14]+.04000666666666666*gradb*t[29]*t[7]*t[5]*t[27]*t[14]-0.06001*gradb*t[29]*t[12]*t[19]*t[28]*t[27]*t[14]+0.12002*gradb*t[29]*t[12]*t[5]*t[52]*t[32]+3.6*t[11]*t[5]*t[27]*t[24]*t[32]+0.225*t[11]*t[41]*t[56]*t[14]*t[18]-2.8*t[53]*t[5]*t[14]*t[18]-0.45*t[11]*t[55]*t[19]*t[14]*t[18]+1.2*t[26]*t[19]*t[28]*t[14]*t[18]+1.8*t[11]*t[54]*t[5]*t[32]*t[18]-4.8*t[26]*t[5]*t[27]*t[32]*t[18]+1.8*t[11]*t[19]*t[28]*t[27]*t[32]*t[18]-5.4*t[11]*t[5]*t[52]*t[68]*t[18];
875 d3fdrbrbab = t[151]+t[99]+t[150]+t[149]+t[148]+t[147]+t[146]+t[145];
876 d3fdrbgaga = t[51]+t[116]+t[58]-0.030005*t[15]*t[60]*t[12]*t[19]*t[28]*t[23]+t[115]+t[57]+0.06001*t[15]*t[60]*t[12]*t[5]*t[27]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+t[114]+0.216036*grada*t[29]*t[10]*t[12]*t[5]*t[27]*t[32]+0.81*t[59]*t[19]*t[28]*t[32]*t[18]-4.86*t[59]*t[5]*t[27]*t[68]*t[18];
877 d3fdrbgagb = t[120]-0.030005*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+.07201199999999999*grada*t[29]*t[26]*t[12]*t[5]*t[14]+t[119]+t[118]+0.06001*grada*gradb*t[60]*t[12]*t[5]*t[27]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+t[117]+0.108018*grada*t[29]*t[11]*t[12]*t[5]*t[27]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[5]*t[27]*t[32]-2.16*t[10]*t[26]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[19]*t[28]*t[32]*t[18]-4.86*t[10]*t[11]*t[5]*t[27]*t[68]*t[18];
878 d3fdrbgaab = t[127]+t[155]+t[124]+t[154]+t[153]+t[152];
879 d3fdrbgbgb = t[51]+t[130]+t[58]-0.030005*t[16]*t[60]*t[12]*t[19]*t[28]*t[23]+0.144024*gradb*t[29]*t[26]*t[12]*t[5]*t[14]+t[129]+t[57]+0.06001*t[16]*t[60]*t[12]*t[5]*t[27]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]+t[128]+0.216036*gradb*t[29]*t[11]*t[12]*t[5]*t[27]*t[32]-4.32*t[156]*t[5]*t[32]*t[18]+0.81*t[63]*t[19]*t[28]*t[32]*t[18]-4.86*t[63]*t[5]*t[27]*t[68]*t[18];
880 d3fdrbgbab = t[136]+t[161]+t[160]+t[134]+t[159]+t[158]+t[157];
881 d3fdrbabab = t[139]+t[163]+t[162];
882 d3fdgagaga = -0.18003*grada*t[60]*t[12]*t[5]*t[23]+0.18003*t[165]*t[166]*t[12]*t[5]*t[23]-0.162027*t[29]*t[10]*t[12]*t[5]*t[14]+0.162027*t[15]*t[60]*t[10]*t[12]*t[5]*t[14]+0.2916486*grada*t[29]*t[59]*t[12]*t[5]*t[32]-4.374*t[164]*t[5]*t[68]*t[18];
883 d3fdgagagb = t[65]+0.18003*t[15]*gradb*t[166]*t[12]*t[5]*t[23]+t[64]+0.054009*t[15]*t[60]*t[11]*t[12]*t[5]*t[14]+0.108018*grada*gradb*t[60]*t[10]*t[12]*t[5]*t[14]+0.1944324*grada*t[29]*t[10]*t[11]*t[12]*t[5]*t[32]+0.0972162*gradb*t[29]*t[59]*t[12]*t[5]*t[32]-4.374*t[59]*t[11]*t[5]*t[68]*t[18];
884 d3fdgagaab = t[66]+t[169]+t[168]+t[167];
885 d3fdgagbgb = t[62]+0.18003*grada*t[16]*t[166]*t[12]*t[5]*t[23]+0.108018*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[14]+t[61]+0.054009*t[16]*t[60]*t[10]*t[12]*t[5]*t[14]+0.1944324*gradb*t[29]*t[10]*t[11]*t[12]*t[5]*t[32]+0.0972162*grada*t[29]*t[63]*t[12]*t[5]*t[32]-4.374*t[10]*t[63]*t[5]*t[68]*t[18];
886 d3fdgagbab = 0.18003*grada*gradb*t[166]*t[12]*t[5]*t[23]+0.054009*grada*t[60]*t[11]*t[12]*t[5]*t[14]+0.054009*gradb*t[60]*t[10]*t[12]*t[5]*t[14]+0.0972162*t[29]*t[10]*t[11]*t[12]*t[5]*t[32];
887 d3fdgaabab = t[171]+t[170];
888 d3fdgbgbgb = -0.18003*gradb*t[60]*t[12]*t[5]*t[23]+0.18003*t[173]*t[166]*t[12]*t[5]*t[23]-0.162027*t[29]*t[11]*t[12]*t[5]*t[14]+0.162027*t[16]*t[60]*t[11]*t[12]*t[5]*t[14]+0.2916486*gradb*t[29]*t[63]*t[12]*t[5]*t[32]-4.374*t[172]*t[5]*t[68]*t[18];
889 d3fdgbgbab = t[66]+t[176]+t[175]+t[174];
890 d3fdgbabab = t[178]+t[177];
891 d3fdababab = t[179];
892 d4fdrararara = 24.0*t[18]*t[180]*t[5]*pow(t[9],4.0)-1.5*t[14]*t[18]*t[77]*t[78]*t[9]-2.0*t[14]*t[18]*t[19]*t[76]*t[9]+8.0*t[18]*t[32]*t[5]*t[72]*t[9]+24.0*t[24]*t[32]*t[36]*t[5]*t[9]+3.0*t[14]*t[24]*t[41]*t[42]*t[9]+3.0*t[14]*t[18]*t[22]*t[40]*t[41]*t[9]-6.0*t[14]*t[19]*t[24]*t[40]*t[9]+12.0*t[18]*t[19]*t[22]*t[32]*t[36]*t[9]+1.5*t[23]*t[24]*t[77]*t[78]+2.25*t[18]*t[23]*t[40]*t[42]*t[77]-1.0*t[18]*t[22]*t[23]*t[41]*t[76]+2.0*t[19]*t[23]*t[24]*t[76]+0.5*t[18]*t[19]*t[23]*(192.0*t[1]*t[73]+t[185]+t[184])-4.0*t[14]*t[24]*t[5]*t[72]-2.0*t[14]*t[18]*t[19]*t[22]*t[72]-24.0*t[24]*t[5]*t[67]*t[68]-12.0*t[18]*t[19]*t[22]*t[67]*t[68]-36.0*t[18]*t[31]*t[36]*t[5]*t[68]+6.0*t[18]*t[32]*pow(t[36],2.0)*t[5]-1.0*t[14]*t[18]*(t[183]+40.44444444444444*grada*t[181])*t[5]+1.5*t[14]*t[18]*t[36]*t[41]*t[42]-3.0*t[18]*t[31]*t[32]*t[41]*t[42]-0.75*t[18]*t[23]*pow(t[40],2.0)*t[41]-3.0*t[22]*t[23]*t[24]*t[40]*t[41]-3.0*t[14]*t[18]*t[19]*t[36]*t[40]+6.0*t[18]*t[19]*t[31]*t[32]*t[40]-6.0*t[14]*t[19]*t[22]*t[24]*t[36]+.3200533333333333*t[17]*t[34]*t[5]*t[31]*t[32]+12.0*t[19]*t[22]*t[24]*t[31]*t[32]-.04000666666666666*t[17]*t[34]*t[41]*t[42]*t[23]-0.9375*t[18]*t[186]*pow(t[22],4.0)*t[23]-.1244651851851852*t[17]*t[70]*t[19]*t[22]*t[23]+.08001333333333333*t[17]*t[34]*t[40]*t[19]*t[23]+t[187]+.2489303703703704*t[17]*t[70]*t[5]*t[9]*t[14]-.1600266666666667*t[17]*t[34]*t[19]*t[22]*t[9]*t[14]-.1600266666666667*t[17]*t[34]*t[36]*t[5]*t[14];
893 d4fdrarararb = -1.5*t[14]*t[18]*t[19]*t[82]*t[9]-1.125*t[14]*t[18]*t[28]*t[42]*t[77]*t[9]-18.0*t[18]*t[27]*t[36]*t[5]*t[68]*t[9]+6.0*t[24]*t[32]*t[36]*t[5]*t[9]+1.5*t[14]*t[18]*t[22]*t[41]*t[44]*t[9]-3.0*t[14]*t[19]*t[24]*t[44]*t[9]-1.5*t[18]*t[27]*t[32]*t[41]*t[42]*t[9]+0.75*t[14]*t[24]*t[41]*t[42]*t[9]+0.75*t[14]*t[18]*t[28]*t[40]*t[41]*t[9]+1.5*t[14]*t[22]*t[24]*t[28]*t[41]*t[9]+3.0*t[18]*t[19]*t[27]*t[32]*t[40]*t[9]-1.5*t[14]*t[19]*t[24]*t[40]*t[9]+3.0*t[18]*t[19]*t[28]*t[32]*t[36]*t[9]+6.0*t[19]*t[22]*t[24]*t[27]*t[32]*t[9]-0.75*t[18]*t[22]*t[23]*t[41]*t[82]+1.5*t[19]*t[23]*t[24]*t[82]-0.375*t[14]*t[18]*t[27]*t[77]*t[78]+0.375*t[23]*t[24]*t[77]*t[78]-0.9375*t[18]*t[186]*t[23]*t[28]*t[78]+1.125*t[18]*t[23]*t[42]*t[44]*t[77]+1.125*t[23]*t[24]*t[28]*t[42]*t[77]+1.125*t[18]*t[22]*t[23]*t[28]*t[40]*t[77]-0.25*t[18]*t[23]*t[28]*t[41]*t[76]-0.5*t[14]*t[18]*t[19]*t[27]*t[76]+0.5*t[19]*t[23]*t[24]*t[76]+0.5*t[18]*t[19]*t[23]*(96.0*t[1]*t[73]+t[184])+2.0*t[18]*t[27]*t[32]*t[5]*t[72]-1.0*t[14]*t[24]*t[5]*t[72]-0.5*t[14]*t[18]*t[19]*t[28]*t[72]-6.0*t[24]*t[5]*t[67]*t[68]-3.0*t[18]*t[19]*t[28]*t[67]*t[68]-18.0*t[24]*t[27]*t[31]*t[5]*t[68]-9.0*t[18]*t[19]*t[22]*t[27]*t[31]*t[68]+24.0*t[18]*t[180]*t[27]*t[5]*t[67]+6.0*t[24]*t[27]*t[32]*t[36]*t[5]-0.75*t[18]*t[23]*t[40]*t[41]*t[44]-1.5*t[22]*t[23]*t[24]*t[41]*t[44]-1.5*t[14]*t[18]*t[19]*t[36]*t[44]+3.0*t[18]*t[19]*t[31]*t[32]*t[44]+0.75*t[14]*t[24]*t[27]*t[41]*t[42]-0.75*t[23]*t[24]*t[28]*t[40]*t[41]+0.75*t[14]*t[18]*t[22]*t[27]*t[40]*t[41]-0.75*t[22]*t[23]*t[24]*t[40]*t[41]+0.75*t[14]*t[18]*t[22]*t[28]*t[36]*t[41]-1.5*t[18]*t[22]*t[28]*t[31]*t[32]*t[41]-1.5*t[14]*t[19]*t[24]*t[27]*t[40]+3.0*t[18]*t[19]*t[22]*t[27]*t[32]*t[36]-1.5*t[14]*t[19]*t[24]*t[28]*t[36]-1.5*t[14]*t[19]*t[22]*t[24]*t[36]+.1600266666666667*t[17]*t[34]*t[5]*t[31]*t[32]+3.0*t[19]*t[24]*t[28]*t[31]*t[32]+3.0*t[19]*t[22]*t[24]*t[31]*t[32]+3.308535951463466*t[34]*t[5]*t[9]*t[24]*t[32]-.02000333333333333*t[17]*t[34]*t[41]*t[42]*t[23]-.03111629629629629*t[17]*t[70]*t[19]*t[28]*t[23]-.09334888888888888*t[17]*t[70]*t[19]*t[22]*t[23]+.04000666666666666*t[17]*t[34]*t[40]*t[19]*t[23]+t[196]+t[195]+t[194]+t[192]+t[191]+t[190]+t[189]+t[188]+t[187]-4.9628039271952*t[34]*t[5]*t[31]*t[68]*t[18]-3.85995861004071*t[70]*t[5]*t[9]*t[32]*t[18]+1.654267975731733*t[34]*t[19]*t[22]*t[9]*t[32]*t[18]+1.654267975731733*t[34]*t[36]*t[5]*t[32]*t[18]+.2067834969664666*t[34]*t[41]*t[42]*t[14]*t[18]+.9649896525101777*t[70]*t[19]*t[22]*t[14]*t[18]-.4135669939329333*t[34]*t[40]*t[19]*t[14]*t[18]+.1866977777777778*t[17]*t[70]*t[5]*t[9]*t[14]-.08001333333333333*t[17]*t[34]*t[19]*t[22]*t[9]*t[14]-.08001333333333333*t[17]*t[34]*t[36]*t[5]*t[14]+.06223259259259259*t[17]*t[70]*t[5]*t[27]*t[14]-.8271339878658666*t[34]*t[19]*t[22]*t[24]*t[14];
894 d4fdrararaga = 0.02250375*grada*t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*grada*t[29]*t[7]*t[41]*t[42]*t[23]+t[198]-0.030005*grada*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*grada*t[29]*t[12]*t[76]*t[19]*t[23]+.04000666666666666*grada*t[29]*t[34]*t[19]*t[22]*t[23]-0.0450075*grada*t[29]*t[12]*t[40]*t[41]*t[22]*t[23]+.09601599999999999*t[17]*t[6]*t[34]*t[5]*t[14]+.06000999999999999*grada*t[29]*t[7]*t[36]*t[5]*t[14]+t[197]-0.06001*grada*t[29]*t[12]*t[72]*t[5]*t[14]+0.675*t[10]*t[41]*t[42]*t[24]*t[14]-8.4*t[33]*t[5]*t[24]*t[14]-1.35*t[10]*t[40]*t[19]*t[24]*t[14]+3.6*t[6]*t[19]*t[22]*t[24]*t[14]+0.0450075*grada*t[29]*t[12]*t[41]*t[42]*t[9]*t[14]-.08001333333333333*grada*t[29]*t[34]*t[5]*t[9]*t[14]-0.090015*grada*t[29]*t[12]*t[40]*t[19]*t[9]*t[14]+.06000999999999999*grada*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.036006*t[17]*t[10]*t[34]*t[19]*t[22]*t[14]-0.090015*grada*t[29]*t[12]*t[36]*t[19]*t[22]*t[14]-0.12002*grada*t[29]*t[7]*t[5]*t[31]*t[32]+0.18003*grada*t[29]*t[12]*t[19]*t[22]*t[31]*t[32]+5.4*t[10]*t[36]*t[5]*t[24]*t[32]-14.4*t[6]*t[5]*t[9]*t[24]*t[32]+5.4*t[10]*t[19]*t[22]*t[9]*t[24]*t[32]+0.144024*t[17]*t[10]*t[34]*t[5]*t[9]*t[32]+0.36006*grada*t[29]*t[12]*t[36]*t[5]*t[9]*t[32]-0.36006*grada*t[29]*t[12]*t[5]*t[67]*t[68]-16.2*t[10]*t[5]*t[31]*t[24]*t[68]-0.3375*t[10]*t[77]*t[78]*t[14]*t[18]-.8999999999999999*t[6]*t[41]*t[42]*t[14]*t[18]+9.333333333333334*t[69]*t[5]*t[14]*t[18]+1.8*t[6]*t[40]*t[19]*t[14]*t[18]-0.45*t[10]*t[76]*t[19]*t[14]*t[18]-4.2*t[33]*t[19]*t[22]*t[14]*t[18]+0.675*t[10]*t[40]*t[41]*t[22]*t[14]*t[18]-7.2*t[6]*t[36]*t[5]*t[32]*t[18]+1.8*t[10]*t[72]*t[5]*t[32]*t[18]-1.35*t[10]*t[41]*t[42]*t[9]*t[32]*t[18]+16.8*t[33]*t[5]*t[9]*t[32]*t[18]+2.7*t[10]*t[40]*t[19]*t[9]*t[32]*t[18]-7.2*t[6]*t[19]*t[22]*t[9]*t[32]*t[18]+2.7*t[10]*t[36]*t[19]*t[22]*t[32]*t[18]+21.6*t[6]*t[5]*t[31]*t[68]*t[18]-8.100000000000001*t[10]*t[19]*t[22]*t[31]*t[68]*t[18]-16.2*t[10]*t[36]*t[5]*t[9]*t[68]*t[18]+21.6*t[10]*t[5]*t[67]*t[180]*t[18];
895 d4fdrararagb = 0.02250375*gradb*t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*gradb*t[29]*t[7]*t[41]*t[42]*t[23]+t[200]-0.030005*gradb*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*gradb*t[29]*t[12]*t[76]*t[19]*t[23]+.04000666666666666*gradb*t[29]*t[34]*t[19]*t[22]*t[23]-0.0450075*gradb*t[29]*t[12]*t[40]*t[41]*t[22]*t[23]+.06000999999999999*gradb*t[29]*t[7]*t[36]*t[5]*t[14]+t[199]-0.06001*gradb*t[29]*t[12]*t[72]*t[5]*t[14]+0.675*t[11]*t[41]*t[42]*t[24]*t[14]-1.35*t[11]*t[40]*t[19]*t[24]*t[14]+0.0450075*gradb*t[29]*t[12]*t[41]*t[42]*t[9]*t[14]-.08001333333333333*gradb*t[29]*t[34]*t[5]*t[9]*t[14]-0.090015*gradb*t[29]*t[12]*t[40]*t[19]*t[9]*t[14]+.06000999999999999*gradb*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.036006*t[17]*t[11]*t[34]*t[19]*t[22]*t[14]-0.090015*gradb*t[29]*t[12]*t[36]*t[19]*t[22]*t[14]-0.12002*gradb*t[29]*t[7]*t[5]*t[31]*t[32]+0.18003*gradb*t[29]*t[12]*t[19]*t[22]*t[31]*t[32]+5.4*t[11]*t[36]*t[5]*t[24]*t[32]+5.4*t[11]*t[19]*t[22]*t[9]*t[24]*t[32]+0.144024*t[17]*t[11]*t[34]*t[5]*t[9]*t[32]+0.36006*gradb*t[29]*t[12]*t[36]*t[5]*t[9]*t[32]-0.36006*gradb*t[29]*t[12]*t[5]*t[67]*t[68]-16.2*t[11]*t[5]*t[31]*t[24]*t[68]-0.3375*t[11]*t[77]*t[78]*t[14]*t[18]-0.45*t[11]*t[76]*t[19]*t[14]*t[18]+0.675*t[11]*t[40]*t[41]*t[22]*t[14]*t[18]+1.8*t[11]*t[72]*t[5]*t[32]*t[18]-1.35*t[11]*t[41]*t[42]*t[9]*t[32]*t[18]+2.7*t[11]*t[40]*t[19]*t[9]*t[32]*t[18]+2.7*t[11]*t[36]*t[19]*t[22]*t[32]*t[18]-8.100000000000001*t[11]*t[19]*t[22]*t[31]*t[68]*t[18]-16.2*t[11]*t[36]*t[5]*t[9]*t[68]*t[18]+21.6*t[11]*t[5]*t[67]*t[180]*t[18];
896 d4fdrararaab = 0.02250375*t[29]*t[12]*t[77]*t[78]*t[23]+0.0150025*t[29]*t[7]*t[41]*t[42]*t[23]+t[201]-0.030005*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*t[29]*t[12]*t[76]*t[19]*t[23]+.04000666666666666*t[29]*t[34]*t[19]*t[22]*t[23]-0.0450075*t[29]*t[12]*t[40]*t[41]*t[22]*t[23]+.06000999999999999*t[29]*t[7]*t[36]*t[5]*t[14]-0.06001*t[29]*t[12]*t[72]*t[5]*t[14]+0.0450075*t[29]*t[12]*t[41]*t[42]*t[9]*t[14]-.08001333333333333*t[29]*t[34]*t[5]*t[9]*t[14]-0.090015*t[29]*t[12]*t[40]*t[19]*t[9]*t[14]+.06000999999999999*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.090015*t[29]*t[12]*t[36]*t[19]*t[22]*t[14]-0.12002*t[29]*t[7]*t[5]*t[31]*t[32]+0.18003*t[29]*t[12]*t[19]*t[22]*t[31]*t[32]+0.36006*t[29]*t[12]*t[36]*t[5]*t[9]*t[32]-0.36006*t[29]*t[12]*t[5]*t[67]*t[68];
897 d4fdrararbrb = -0.75*t[14]*t[18]*t[22]*t[56]*t[77]*t[9]-12.0*t[24]*t[5]*t[52]*t[68]*t[9]-6.0*t[18]*t[19]*t[22]*t[52]*t[68]*t[9]+0.5*t[14]*t[24]*t[41]*t[56]*t[9]+0.5*t[14]*t[18]*t[22]*t[41]*t[55]*t[9]-1.0*t[14]*t[19]*t[24]*t[55]*t[9]+4.0*t[24]*t[32]*t[5]*t[54]*t[9]+2.0*t[18]*t[19]*t[22]*t[32]*t[54]*t[9]+4.0*t[18]*t[19]*t[27]*t[32]*t[44]*t[9]-2.0*t[14]*t[19]*t[24]*t[44]*t[9]-2.0*t[18]*t[22]*t[27]*t[28]*t[32]*t[41]*t[9]+4.0*t[19]*t[24]*t[27]*t[28]*t[32]*t[9]+4.0*t[19]*t[22]*t[24]*t[27]*t[32]*t[9]-1.0*t[101]*t[14]*t[18]*t[19]*t[9]-0.5*t[18]*t[23]*t[28]*t[41]*t[82]-1.0*t[14]*t[18]*t[19]*t[27]*t[82]+0.375*t[18]*t[23]*t[40]*t[56]*t[77]+0.75*t[22]*t[23]*t[24]*t[56]*t[77]+0.375*t[18]*t[23]*t[42]*t[55]*t[77]+1.5*t[18]*t[22]*t[23]*t[28]*t[44]*t[77]-0.75*t[14]*t[18]*t[27]*t[28]*t[42]*t[77]+0.75*t[23]*t[24]*t[28]*t[42]*t[77]-6.0*t[18]*t[31]*t[5]*t[54]*t[68]-6.0*t[18]*t[36]*t[5]*t[52]*t[68]-12.0*t[24]*t[27]*t[31]*t[5]*t[68]-6.0*t[18]*t[19]*t[27]*t[28]*t[31]*t[68]-0.9375*t[18]*t[186]*t[23]*t[42]*t[56]+0.25*t[14]*t[18]*t[36]*t[41]*t[56]-0.5*t[18]*t[31]*t[32]*t[41]*t[56]-0.25*t[18]*t[23]*t[40]*t[41]*t[55]-0.5*t[22]*t[23]*t[24]*t[41]*t[55]-0.5*t[14]*t[18]*t[19]*t[36]*t[55]+2.0*t[18]*t[32]*t[36]*t[5]*t[54]+0.25*t[14]*t[18]*t[41]*t[42]*t[54]-0.5*t[14]*t[18]*t[19]*t[40]*t[54]-1.0*t[14]*t[19]*t[22]*t[24]*t[54]+24.0*t[18]*t[180]*t[31]*t[5]*t[52]-0.5*t[18]*t[32]*t[41]*t[42]*t[52]+2.0*t[19]*t[22]*t[24]*t[32]*t[52]+4.0*t[24]*t[27]*t[32]*t[36]*t[5]-0.5*t[18]*t[23]*t[41]*pow(t[44],2.0)-1.0*t[23]*t[24]*t[28]*t[41]*t[44]-1.0*t[22]*t[23]*t[24]*t[41]*t[44]-2.0*t[14]*t[19]*t[24]*t[27]*t[44]+0.5*t[14]*t[24]*t[27]*t[41]*t[42]+0.5*t[14]*t[18]*t[27]*t[28]*t[40]*t[41]-0.5*t[23]*t[24]*t[28]*t[40]*t[41]-0.5*t[101]*t[18]*t[22]*t[23]*t[41]-1.0*t[14]*t[19]*t[24]*t[27]*t[40]+0.5*t[18]*t[19]*t[23]*(24.0*t[37]+t[184])+2.0*t[18]*t[19]*t[27]*t[28]*t[32]*t[36]-1.0*t[14]*t[19]*t[24]*t[28]*t[36]+.05334222222222221*t[17]*t[34]*t[5]*t[52]*t[32]+.05334222222222221*t[17]*t[34]*t[5]*t[31]*t[32]+2.0*t[19]*t[24]*t[28]*t[31]*t[32]+.2133688888888889*t[17]*t[34]*t[5]*t[9]*t[27]*t[32]+2.205690634308977*t[34]*t[5]*t[9]*t[24]*t[32]+2.205690634308977*t[34]*t[5]*t[27]*t[24]*t[32]-.006667777777777777*t[17]*t[34]*t[41]*t[56]*t[23]-.006667777777777777*t[17]*t[34]*t[41]*t[42]*t[23]-.02667111111111111*t[17]*t[34]*t[41]*t[22]*t[28]*t[23]-.06223259259259259*t[17]*t[70]*t[19]*t[28]*t[23]+t[82]*t[19]*t[24]*t[23]+t[101]*t[19]*t[24]*t[23]-.06223259259259259*t[17]*t[70]*t[19]*t[22]*t[23]+.01333555555555555*t[17]*t[34]*t[55]*t[19]*t[23]+.05334222222222221*t[17]*t[34]*t[44]*t[19]*t[23]+.01333555555555555*t[17]*t[34]*t[40]*t[19]*t[23]+t[188]+t[187]-6.617071902926932*t[34]*t[5]*t[9]*t[27]*t[68]*t[18]-2.57330574002714*t[70]*t[5]*t[9]*t[32]*t[18]+1.102845317154489*t[34]*t[19]*t[28]*t[9]*t[32]*t[18]+t[40]*t[19]*t[52]*t[32]*t[18]+.3040669483923962*t[193]*t[5]*t[32]*t[18]+t[55]*t[19]*t[31]*t[32]*t[18]-2.57330574002714*t[70]*t[5]*t[27]*t[32]*t[18]+1.102845317154489*t[34]*t[19]*t[22]*t[27]*t[32]*t[18]+t[44]*t[41]*t[28]*t[9]*t[14]*t[18]+.2757113292886222*t[34]*t[41]*t[22]*t[28]*t[14]*t[18]+.6433264350067851*t[70]*t[19]*t[28]*t[14]*t[18]+t[44]*t[41]*t[22]*t[27]*t[14]*t[18]+.6433264350067851*t[70]*t[19]*t[22]*t[14]*t[18]-.5514226585772444*t[34]*t[44]*t[19]*t[14]*t[18]+.1244651851851852*t[17]*t[70]*t[5]*t[9]*t[14]-.05334222222222221*t[17]*t[34]*t[19]*t[28]*t[9]*t[14]-.02667111111111111*t[17]*t[34]*t[19]*t[22]*t[9]*t[14]-.02667111111111111*t[17]*t[34]*t[54]*t[5]*t[14]-.02667111111111111*t[17]*t[34]*t[36]*t[5]*t[14]-.02941410999219594*t[17]*t[193]*t[5]*t[14]+.1244651851851852*t[17]*t[70]*t[5]*t[27]*t[14]-.02667111111111111*t[17]*t[34]*t[19]*t[28]*t[27]*t[14]-.05334222222222221*t[17]*t[34]*t[19]*t[22]*t[27]*t[14]+t[41]*t[22]*t[28]*t[9]*t[24]*t[14]+2.57330574002714*t[70]*t[5]*t[24]*t[14]-.5514226585772444*t[34]*t[19]*t[28]*t[24]*t[14]+t[41]*t[22]*t[28]*t[27]*t[24]*t[14]-.5514226585772444*t[34]*t[19]*t[22]*t[24]*t[14];
898 d4fdrararbga = .005000833333333333*grada*t[29]*t[7]*t[41]*t[42]*t[23]+t[198]+t[217]-.01000166666666667*grada*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*grada*t[29]*t[12]*t[82]*t[19]*t[23]+0.02250375*grada*t[29]*t[12]*t[77]*t[42]*t[28]*t[23]+.01333555555555555*grada*t[29]*t[34]*t[19]*t[28]*t[23]-0.0150025*grada*t[29]*t[12]*t[40]*t[41]*t[28]*t[23]+t[216]+.02667111111111111*grada*t[29]*t[34]*t[19]*t[22]*t[23]-0.030005*grada*t[29]*t[12]*t[44]*t[41]*t[22]*t[23]+.06401066666666666*t[17]*t[6]*t[34]*t[5]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[36]*t[5]*t[14]+t[197]+t[215]+0.225*t[10]*t[41]*t[42]*t[24]*t[14]-2.8*t[33]*t[5]*t[24]*t[14]+t[213]-0.45*t[10]*t[40]*t[19]*t[24]*t[14]+t[212]+t[211]+1.2*t[6]*t[19]*t[22]*t[24]*t[14]+0.0150025*grada*t[29]*t[12]*t[41]*t[42]*t[27]*t[14]-.02667111111111111*grada*t[29]*t[34]*t[5]*t[27]*t[14]-0.030005*grada*t[29]*t[12]*t[40]*t[19]*t[27]*t[14]+t[210]-.05334222222222221*grada*t[29]*t[34]*t[5]*t[9]*t[14]-0.06001*grada*t[29]*t[12]*t[44]*t[19]*t[9]*t[14]+t[209]+0.030005*grada*t[29]*t[12]*t[41]*t[22]*t[28]*t[9]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.012002*t[17]*t[10]*t[34]*t[19]*t[28]*t[14]-0.030005*grada*t[29]*t[12]*t[36]*t[19]*t[28]*t[14]-0.024004*t[17]*t[10]*t[34]*t[19]*t[22]*t[14]-.01654543687061022*grada*t[29]*t[102]*t[19]*t[22]*t[14]-.04000666666666666*grada*t[29]*t[7]*t[5]*t[31]*t[32]+0.06001*grada*t[29]*t[12]*t[19]*t[28]*t[31]*t[32]+t[208]+1.8*t[10]*t[36]*t[5]*t[24]*t[32]+t[207]+t[206]-4.8*t[6]*t[5]*t[9]*t[24]*t[32]+t[205]+1.8*t[10]*t[19]*t[22]*t[9]*t[24]*t[32]+0.048008*t[17]*t[10]*t[34]*t[5]*t[27]*t[32]+0.12002*grada*t[29]*t[12]*t[36]*t[5]*t[27]*t[32]+t[204]+0.12002*grada*t[29]*t[12]*t[19]*t[22]*t[9]*t[27]*t[32]+.09601599999999999*t[17]*t[10]*t[34]*t[5]*t[9]*t[32]+.06618174748244086*grada*t[29]*t[102]*t[5]*t[9]*t[32]-5.4*t[10]*t[5]*t[31]*t[24]*t[68]+t[203]-0.36006*grada*t[29]*t[12]*t[5]*t[31]*t[27]*t[68]+1.2*t[6]*t[44]*t[19]*t[14]*t[18]-0.45*t[10]*t[82]*t[19]*t[14]*t[18]-0.3375*t[10]*t[77]*t[42]*t[28]*t[14]*t[18]-1.4*t[33]*t[19]*t[28]*t[14]*t[18]+0.225*t[10]*t[40]*t[41]*t[28]*t[14]*t[18]-0.6*t[6]*t[41]*t[22]*t[28]*t[14]*t[18]+0.45*t[10]*t[44]*t[41]*t[22]*t[14]*t[18]-1.323414380585386*t[6]*t[34]*t[5]*t[32]*t[18]+t[202]-0.45*t[10]*t[41]*t[42]*t[27]*t[32]*t[18]+5.6*t[33]*t[5]*t[27]*t[32]*t[18]+0.9*t[10]*t[40]*t[19]*t[27]*t[32]*t[18]-2.4*t[6]*t[19]*t[22]*t[27]*t[32]*t[18]+1.8*t[10]*t[44]*t[19]*t[9]*t[32]*t[18]-2.4*t[6]*t[19]*t[28]*t[9]*t[32]*t[18]-0.9*t[10]*t[41]*t[22]*t[28]*t[9]*t[32]*t[18]+0.9*t[10]*t[36]*t[19]*t[28]*t[32]*t[18]+0.49628039271952*t[10]*t[34]*t[19]*t[22]*t[32]*t[18]-2.7*t[10]*t[19]*t[28]*t[31]*t[68]*t[18]-5.4*t[10]*t[36]*t[5]*t[27]*t[68]*t[18]+14.4*t[6]*t[5]*t[9]*t[27]*t[68]*t[18]-5.4*t[10]*t[19]*t[22]*t[9]*t[27]*t[68]*t[18]-2.97768235631712*t[10]*t[34]*t[5]*t[9]*t[68]*t[18]+21.6*t[10]*t[5]*t[31]*t[27]*t[180]*t[18];
899 d4fdrararbgb = .005000833333333333*gradb*t[29]*t[7]*t[41]*t[42]*t[23]+t[200]+t[232]-.01000166666666667*gradb*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*gradb*t[29]*t[12]*t[82]*t[19]*t[23]+0.02250375*gradb*t[29]*t[12]*t[77]*t[42]*t[28]*t[23]+.01333555555555555*gradb*t[29]*t[34]*t[19]*t[28]*t[23]-0.0150025*gradb*t[29]*t[12]*t[40]*t[41]*t[28]*t[23]+t[231]+.02667111111111111*gradb*t[29]*t[34]*t[19]*t[22]*t[23]-0.030005*gradb*t[29]*t[12]*t[44]*t[41]*t[22]*t[23]+.03200533333333333*t[17]*t[26]*t[34]*t[5]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[36]*t[5]*t[14]+t[199]+t[230]+0.225*t[11]*t[41]*t[42]*t[24]*t[14]+t[229]-0.45*t[11]*t[40]*t[19]*t[24]*t[14]+t[228]+t[227]+0.0150025*gradb*t[29]*t[12]*t[41]*t[42]*t[27]*t[14]-.02667111111111111*gradb*t[29]*t[34]*t[5]*t[27]*t[14]-0.030005*gradb*t[29]*t[12]*t[40]*t[19]*t[27]*t[14]+t[226]-.05334222222222221*gradb*t[29]*t[34]*t[5]*t[9]*t[14]-0.06001*gradb*t[29]*t[12]*t[44]*t[19]*t[9]*t[14]+t[225]+0.030005*gradb*t[29]*t[12]*t[41]*t[22]*t[28]*t[9]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.012002*t[17]*t[11]*t[34]*t[19]*t[28]*t[14]-0.030005*gradb*t[29]*t[12]*t[36]*t[19]*t[28]*t[14]-0.024004*t[17]*t[11]*t[34]*t[19]*t[22]*t[14]-.01654543687061022*gradb*t[29]*t[102]*t[19]*t[22]*t[14]-.04000666666666666*gradb*t[29]*t[7]*t[5]*t[31]*t[32]+0.06001*gradb*t[29]*t[12]*t[19]*t[28]*t[31]*t[32]+t[224]+1.8*t[11]*t[36]*t[5]*t[24]*t[32]+t[223]+t[222]+t[221]+1.8*t[11]*t[19]*t[22]*t[9]*t[24]*t[32]+0.048008*t[17]*t[11]*t[34]*t[5]*t[27]*t[32]+0.12002*gradb*t[29]*t[12]*t[36]*t[5]*t[27]*t[32]+t[220]+0.12002*gradb*t[29]*t[12]*t[19]*t[22]*t[9]*t[27]*t[32]+.09601599999999999*t[17]*t[11]*t[34]*t[5]*t[9]*t[32]+.06618174748244086*gradb*t[29]*t[102]*t[5]*t[9]*t[32]-5.4*t[11]*t[5]*t[31]*t[24]*t[68]+t[219]-0.36006*gradb*t[29]*t[12]*t[5]*t[31]*t[27]*t[68]-0.3*t[26]*t[41]*t[42]*t[14]*t[18]+0.6*t[26]*t[40]*t[19]*t[14]*t[18]-0.45*t[11]*t[82]*t[19]*t[14]*t[18]-0.3375*t[11]*t[77]*t[42]*t[28]*t[14]*t[18]+0.225*t[11]*t[40]*t[41]*t[28]*t[14]*t[18]+0.45*t[11]*t[44]*t[41]*t[22]*t[14]*t[18]-2.4*t[26]*t[36]*t[5]*t[32]*t[18]+t[218]-0.45*t[11]*t[41]*t[42]*t[27]*t[32]*t[18]+0.9*t[11]*t[40]*t[19]*t[27]*t[32]*t[18]+1.8*t[11]*t[44]*t[19]*t[9]*t[32]*t[18]-0.9*t[11]*t[41]*t[22]*t[28]*t[9]*t[32]*t[18]-2.4*t[26]*t[19]*t[22]*t[9]*t[32]*t[18]+0.9*t[11]*t[36]*t[19]*t[28]*t[32]*t[18]+0.49628039271952*t[11]*t[34]*t[19]*t[22]*t[32]*t[18]+7.2*t[26]*t[5]*t[31]*t[68]*t[18]-2.7*t[11]*t[19]*t[28]*t[31]*t[68]*t[18]-5.4*t[11]*t[36]*t[5]*t[27]*t[68]*t[18]-5.4*t[11]*t[19]*t[22]*t[9]*t[27]*t[68]*t[18]-2.97768235631712*t[11]*t[34]*t[5]*t[9]*t[68]*t[18]+21.6*t[11]*t[5]*t[31]*t[27]*t[180]*t[18];
900 d4fdrararbab = .005000833333333333*t[29]*t[7]*t[41]*t[42]*t[23]+t[201]+t[238]-.01000166666666667*t[29]*t[7]*t[40]*t[19]*t[23]+0.030005*t[29]*t[12]*t[82]*t[19]*t[23]+0.02250375*t[29]*t[12]*t[77]*t[42]*t[28]*t[23]+.01333555555555555*t[29]*t[34]*t[19]*t[28]*t[23]-0.0150025*t[29]*t[12]*t[40]*t[41]*t[28]*t[23]+t[237]+.02667111111111111*t[29]*t[34]*t[19]*t[22]*t[23]-0.030005*t[29]*t[12]*t[44]*t[41]*t[22]*t[23]+.02000333333333333*t[29]*t[7]*t[36]*t[5]*t[14]+t[236]+0.0150025*t[29]*t[12]*t[41]*t[42]*t[27]*t[14]-.02667111111111111*t[29]*t[34]*t[5]*t[27]*t[14]-0.030005*t[29]*t[12]*t[40]*t[19]*t[27]*t[14]+t[235]-.05334222222222221*t[29]*t[34]*t[5]*t[9]*t[14]-0.06001*t[29]*t[12]*t[44]*t[19]*t[9]*t[14]+t[234]+0.030005*t[29]*t[12]*t[41]*t[22]*t[28]*t[9]*t[14]+.02000333333333333*t[29]*t[7]*t[19]*t[22]*t[9]*t[14]-0.030005*t[29]*t[12]*t[36]*t[19]*t[28]*t[14]-.01654543687061022*t[29]*t[102]*t[19]*t[22]*t[14]-.04000666666666666*t[29]*t[7]*t[5]*t[31]*t[32]+0.06001*t[29]*t[12]*t[19]*t[28]*t[31]*t[32]+0.12002*t[29]*t[12]*t[36]*t[5]*t[27]*t[32]+t[233]+0.12002*t[29]*t[12]*t[19]*t[22]*t[9]*t[27]*t[32]+.06618174748244086*t[29]*t[102]*t[5]*t[9]*t[32]-0.36006*t[29]*t[12]*t[5]*t[31]*t[27]*t[68];
901 d4fdraragaga = 15.84*t[18]*t[32]*t[5]/pow(rhoa,4.666666666666667)+t[99]+t[98]+t[97]+t[96]+t[95]+t[94]+t[93]-0.648108*grada*t[29]*t[10]*t[12]*t[5]*t[31]*t[68]-9.72*t[59]*t[5]*t[9]*t[24]*t[68]-0.144024*grada*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]-.5760959999999999*grada*t[29]*t[6]*t[12]*t[5]*t[9]*t[32]+0.216036*grada*t[29]*t[10]*t[12]*t[19]*t[22]*t[9]*t[32]+0.216036*grada*t[29]*t[10]*t[12]*t[36]*t[5]*t[32]-0.12002*t[15]*t[60]*t[12]*t[5]*t[31]*t[32]-8.64*t[113]*t[5]*t[24]*t[32]+1.62*t[59]*t[19]*t[22]*t[24]*t[32]+t[241]+t[240]+t[239]+0.0150025*t[15]*t[60]*t[12]*t[41]*t[42]*t[23]+.02000333333333333*t[15]*t[60]*t[7]*t[19]*t[22]*t[23]-0.030005*t[15]*t[60]*t[12]*t[40]*t[19]*t[23]+25.92*t[113]*t[5]*t[9]*t[68]*t[18]-4.86*t[59]*t[19]*t[22]*t[9]*t[68]*t[18]-4.86*t[59]*t[36]*t[5]*t[68]*t[18]-0.405*t[59]*t[41]*t[42]*t[32]*t[18]-4.32*t[113]*t[19]*t[22]*t[32]*t[18]+0.81*t[59]*t[40]*t[19]*t[32]*t[18]+19.44*t[59]*t[5]*t[31]*t[180]*t[18]-.04000666666666666*t[15]*t[60]*t[7]*t[5]*t[9]*t[14]+0.06001*t[15]*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]-.09601599999999999*grada*t[29]*t[6]*t[7]*t[5]*t[14]+0.06001*t[15]*t[60]*t[12]*t[36]*t[5]*t[14]-0.336056*grada*t[29]*t[33]*t[12]*t[5]*t[14]+0.0270045*grada*t[29]*t[10]*t[12]*t[41]*t[42]*t[14]+0.036006*grada*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]+0.144024*grada*t[29]*t[6]*t[12]*t[19]*t[22]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[40]*t[19]*t[14]+t[100];
902 d4fdraragagb = 0.0150025*grada*gradb*t[60]*t[12]*t[41]*t[42]*t[23]+t[245]-0.030005*grada*gradb*t[60]*t[12]*t[40]*t[19]*t[23]+.02000333333333333*grada*gradb*t[60]*t[7]*t[19]*t[22]*t[23]+0.01350225*grada*t[29]*t[11]*t[12]*t[41]*t[42]*t[14]+0.01350225*gradb*t[29]*t[10]*t[12]*t[41]*t[42]*t[14]-0.168028*gradb*t[29]*t[33]*t[12]*t[5]*t[14]-0.048008*gradb*t[29]*t[6]*t[7]*t[5]*t[14]+t[244]+t[243]+0.06001*grada*gradb*t[60]*t[12]*t[36]*t[5]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[40]*t[19]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[40]*t[19]*t[14]-.04000666666666666*grada*gradb*t[60]*t[7]*t[5]*t[9]*t[14]+0.06001*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]+.07201199999999999*gradb*t[29]*t[6]*t[12]*t[19]*t[22]*t[14]+0.018003*grada*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]+0.018003*gradb*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]-0.12002*grada*gradb*t[60]*t[12]*t[5]*t[31]*t[32]+t[242]+0.108018*grada*t[29]*t[11]*t[12]*t[36]*t[5]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[36]*t[5]*t[32]-4.32*t[6]*t[11]*t[5]*t[24]*t[32]+1.62*t[10]*t[11]*t[19]*t[22]*t[24]*t[32]-0.288048*gradb*t[29]*t[6]*t[12]*t[5]*t[9]*t[32]-.07201199999999999*grada*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]-.07201199999999999*gradb*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]+0.108018*grada*t[29]*t[11]*t[12]*t[19]*t[22]*t[9]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[19]*t[22]*t[9]*t[32]-0.324054*grada*t[29]*t[11]*t[12]*t[5]*t[31]*t[68]-0.324054*gradb*t[29]*t[10]*t[12]*t[5]*t[31]*t[68]-9.72*t[10]*t[11]*t[5]*t[9]*t[24]*t[68]-0.405*t[10]*t[11]*t[41]*t[42]*t[32]*t[18]+5.040000000000001*t[33]*t[11]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[40]*t[19]*t[32]*t[18]-2.16*t[6]*t[11]*t[19]*t[22]*t[32]*t[18]-4.86*t[10]*t[11]*t[36]*t[5]*t[68]*t[18]+12.96*t[6]*t[11]*t[5]*t[9]*t[68]*t[18]-4.86*t[10]*t[11]*t[19]*t[22]*t[9]*t[68]*t[18]+19.44*t[10]*t[11]*t[5]*t[31]*t[180]*t[18];
903 d4fdraragaab = 0.0150025*grada*t[60]*t[12]*t[41]*t[42]*t[23]+t[247]-0.030005*grada*t[60]*t[12]*t[40]*t[19]*t[23]+.02000333333333333*grada*t[60]*t[7]*t[19]*t[22]*t[23]+0.01350225*t[29]*t[10]*t[12]*t[41]*t[42]*t[14]-0.168028*t[29]*t[33]*t[12]*t[5]*t[14]-0.048008*t[29]*t[6]*t[7]*t[5]*t[14]+t[246]+0.06001*grada*t[60]*t[12]*t[36]*t[5]*t[14]-0.0270045*t[29]*t[10]*t[12]*t[40]*t[19]*t[14]-.04000666666666666*grada*t[60]*t[7]*t[5]*t[9]*t[14]+0.06001*grada*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]+.07201199999999999*t[29]*t[6]*t[12]*t[19]*t[22]*t[14]+0.018003*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]-0.12002*grada*t[60]*t[12]*t[5]*t[31]*t[32]+0.108018*t[29]*t[10]*t[12]*t[36]*t[5]*t[32]-0.288048*t[29]*t[6]*t[12]*t[5]*t[9]*t[32]-.07201199999999999*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]+0.108018*t[29]*t[10]*t[12]*t[19]*t[22]*t[9]*t[32]-0.324054*t[29]*t[10]*t[12]*t[5]*t[31]*t[68];
904 d4fdraragbgb = t[100]+0.0150025*t[16]*t[60]*t[12]*t[41]*t[42]*t[23]+t[99]+t[250]+t[98]-0.030005*t[16]*t[60]*t[12]*t[40]*t[19]*t[23]+t[97]+.02000333333333333*t[16]*t[60]*t[7]*t[19]*t[22]*t[23]+0.0270045*gradb*t[29]*t[11]*t[12]*t[41]*t[42]*t[14]+t[249]+t[96]+0.06001*t[16]*t[60]*t[12]*t[36]*t[5]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[40]*t[19]*t[14]+t[95]-.04000666666666666*t[16]*t[60]*t[7]*t[5]*t[9]*t[14]+t[94]+0.06001*t[16]*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]+0.036006*gradb*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]+t[93]-0.12002*t[16]*t[60]*t[12]*t[5]*t[31]*t[32]+t[248]+0.216036*gradb*t[29]*t[11]*t[12]*t[36]*t[5]*t[32]+1.62*t[63]*t[19]*t[22]*t[24]*t[32]-0.144024*gradb*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]+0.216036*gradb*t[29]*t[11]*t[12]*t[19]*t[22]*t[9]*t[32]-0.648108*gradb*t[29]*t[11]*t[12]*t[5]*t[31]*t[68]-9.72*t[63]*t[5]*t[9]*t[24]*t[68]-0.405*t[63]*t[41]*t[42]*t[32]*t[18]+0.81*t[63]*t[40]*t[19]*t[32]*t[18]-4.86*t[63]*t[36]*t[5]*t[68]*t[18]-4.86*t[63]*t[19]*t[22]*t[9]*t[68]*t[18]+19.44*t[63]*t[5]*t[31]*t[180]*t[18];
905 d4fdraragbab = 0.0150025*gradb*t[60]*t[12]*t[41]*t[42]*t[23]+t[252]-0.030005*gradb*t[60]*t[12]*t[40]*t[19]*t[23]+.02000333333333333*gradb*t[60]*t[7]*t[19]*t[22]*t[23]+0.01350225*t[29]*t[11]*t[12]*t[41]*t[42]*t[14]+t[251]+0.06001*gradb*t[60]*t[12]*t[36]*t[5]*t[14]-0.0270045*t[29]*t[11]*t[12]*t[40]*t[19]*t[14]-.04000666666666666*gradb*t[60]*t[7]*t[5]*t[9]*t[14]+0.06001*gradb*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]+0.018003*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]-0.12002*gradb*t[60]*t[12]*t[5]*t[31]*t[32]+0.108018*t[29]*t[11]*t[12]*t[36]*t[5]*t[32]-.07201199999999999*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]+0.108018*t[29]*t[11]*t[12]*t[19]*t[22]*t[9]*t[32]-0.324054*t[29]*t[11]*t[12]*t[5]*t[31]*t[68];
906 d4fdraraabab = 0.0150025*t[60]*t[12]*t[41]*t[42]*t[23]+t[253]-0.030005*t[60]*t[12]*t[40]*t[19]*t[23]+.02000333333333333*t[60]*t[7]*t[19]*t[22]*t[23]+0.06001*t[60]*t[12]*t[36]*t[5]*t[14]-.04000666666666666*t[60]*t[7]*t[5]*t[9]*t[14]+0.06001*t[60]*t[12]*t[19]*t[22]*t[9]*t[14]-0.12002*t[60]*t[12]*t[5]*t[31]*t[32];
907 d4fdrarbrbrb = -0.375*t[14]*t[144]*t[18]*t[77]*t[9]-18.0*t[18]*t[27]*t[5]*t[54]*t[68]*t[9]-18.0*t[24]*t[5]*t[52]*t[68]*t[9]-9.0*t[18]*t[19]*t[28]*t[52]*t[68]*t[9]-1.5*t[18]*t[27]*t[32]*t[41]*t[56]*t[9]+0.75*t[14]*t[24]*t[41]*t[56]*t[9]+0.75*t[14]*t[18]*t[28]*t[41]*t[55]*t[9]+3.0*t[18]*t[19]*t[27]*t[32]*t[55]*t[9]-1.5*t[14]*t[19]*t[24]*t[55]*t[9]+6.0*t[24]*t[32]*t[5]*t[54]*t[9]+3.0*t[18]*t[19]*t[28]*t[32]*t[54]*t[9]+2.0*t[142]*t[18]*t[32]*t[5]*t[9]+24.0*t[140]*t[18]*t[180]*t[5]*t[9]+6.0*t[19]*t[24]*t[27]*t[28]*t[32]*t[9]-0.5*t[14]*t[143]*t[18]*t[19]*t[9]+1.125*t[18]*t[23]*t[44]*t[56]*t[77]-1.125*t[14]*t[18]*t[22]*t[27]*t[56]*t[77]+1.125*t[22]*t[23]*t[24]*t[56]*t[77]+1.125*t[18]*t[22]*t[23]*t[28]*t[55]*t[77]+0.375*t[144]*t[23]*t[24]*t[77]+0.5*t[18]*t[19]*t[23]*(t[184]-96.0*t[1]*t[73])-6.0*t[140]*t[24]*t[5]*t[68]-3.0*t[140]*t[18]*t[19]*t[22]*t[68]+0.75*t[14]*t[24]*t[27]*t[41]*t[56]-0.75*t[18]*t[23]*t[41]*t[44]*t[55]-0.75*t[23]*t[24]*t[28]*t[41]*t[55]+0.75*t[14]*t[18]*t[22]*t[27]*t[41]*t[55]-0.75*t[22]*t[23]*t[24]*t[41]*t[55]-1.5*t[14]*t[19]*t[24]*t[27]*t[55]+6.0*t[24]*t[27]*t[32]*t[5]*t[54]-1.5*t[14]*t[18]*t[19]*t[44]*t[54]+0.75*t[14]*t[18]*t[22]*t[28]*t[41]*t[54]+3.0*t[18]*t[19]*t[22]*t[27]*t[32]*t[54]-1.5*t[14]*t[19]*t[24]*t[28]*t[54]-1.5*t[14]*t[19]*t[22]*t[24]*t[54]+3.0*t[18]*t[19]*t[32]*t[44]*t[52]-1.5*t[18]*t[22]*t[28]*t[32]*t[41]*t[52]+3.0*t[19]*t[24]*t[28]*t[32]*t[52]+3.0*t[19]*t[22]*t[24]*t[32]*t[52]-1.0*t[14]*t[142]*t[24]*t[5]+1.5*t[14]*t[18]*t[27]*t[28]*t[41]*t[44]-1.5*t[23]*t[24]*t[28]*t[41]*t[44]-3.0*t[14]*t[19]*t[24]*t[27]*t[44]+1.5*t[14]*t[22]*t[24]*t[27]*t[28]*t[41]-0.75*t[101]*t[18]*t[23]*t[28]*t[41]-0.25*t[143]*t[18]*t[22]*t[23]*t[41]+.1600266666666667*t[17]*t[34]*t[5]*t[52]*t[32]+3.308535951463466*t[34]*t[5]*t[27]*t[24]*t[32]-1.5*t[101]*t[14]*t[18]*t[19]*t[27]+0.5*t[143]*t[19]*t[23]*t[24]+1.5*t[101]*t[19]*t[23]*t[24]-.02000333333333333*t[17]*t[34]*t[41]*t[56]*t[23]-.09334888888888888*t[17]*t[70]*t[19]*t[28]*t[23]-.03111629629629629*t[17]*t[70]*t[19]*t[22]*t[23]-0.9375*t[144]*t[18]*t[186]*t[22]*t[23]+.04000666666666666*t[17]*t[34]*t[55]*t[19]*t[23]-0.5*t[14]*t[142]*t[18]*t[19]*t[22]+t[196]+t[195]+t[194]+t[192]+t[191]+t[190]+t[189]+t[188]+t[187]-4.9628039271952*t[34]*t[5]*t[52]*t[68]*t[18]+1.654267975731733*t[34]*t[54]*t[5]*t[32]*t[18]-3.85995861004071*t[70]*t[5]*t[27]*t[32]*t[18]+1.654267975731733*t[34]*t[19]*t[28]*t[27]*t[32]*t[18]+.2067834969664666*t[34]*t[41]*t[56]*t[14]*t[18]+.9649896525101777*t[70]*t[19]*t[28]*t[14]*t[18]-.4135669939329333*t[34]*t[55]*t[19]*t[14]*t[18]+.06223259259259259*t[17]*t[70]*t[5]*t[9]*t[14]-.08001333333333333*t[17]*t[34]*t[54]*t[5]*t[14]+.1866977777777778*t[17]*t[70]*t[5]*t[27]*t[14]-.08001333333333333*t[17]*t[34]*t[19]*t[28]*t[27]*t[14]-.8271339878658666*t[34]*t[19]*t[28]*t[24]*t[14];
908 d4fdrarbrbga = .005000833333333333*grada*t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*grada*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[198]+t[217]-.01000166666666667*grada*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*grada*t[29]*t[12]*t[101]*t[19]*t[23]+.02667111111111111*grada*t[29]*t[34]*t[19]*t[28]*t[23]-0.030005*grada*t[29]*t[12]*t[44]*t[41]*t[28]*t[23]+t[216]+.01333555555555555*grada*t[29]*t[34]*t[19]*t[22]*t[23]-0.0150025*grada*t[29]*t[12]*t[55]*t[41]*t[22]*t[23]+.03200533333333333*t[17]*t[6]*t[34]*t[5]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[54]*t[5]*t[14]+t[197]+t[215]+0.225*t[10]*t[41]*t[56]*t[24]*t[14]+t[213]-0.45*t[10]*t[55]*t[19]*t[24]*t[14]+t[212]+t[211]-.05334222222222221*grada*t[29]*t[34]*t[5]*t[27]*t[14]-0.06001*grada*t[29]*t[12]*t[44]*t[19]*t[27]*t[14]+.02000333333333333*grada*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]+0.030005*grada*t[29]*t[12]*t[41]*t[22]*t[28]*t[27]*t[14]+t[210]+0.0150025*grada*t[29]*t[12]*t[41]*t[56]*t[9]*t[14]-.02667111111111111*grada*t[29]*t[34]*t[5]*t[9]*t[14]-0.030005*grada*t[29]*t[12]*t[55]*t[19]*t[9]*t[14]+t[209]-0.024004*t[17]*t[10]*t[34]*t[19]*t[28]*t[14]-.01654543687061022*grada*t[29]*t[102]*t[19]*t[28]*t[14]-0.012002*t[17]*t[10]*t[34]*t[19]*t[22]*t[14]-0.030005*grada*t[29]*t[12]*t[54]*t[19]*t[22]*t[14]-.04000666666666666*grada*t[29]*t[7]*t[5]*t[52]*t[32]+0.06001*grada*t[29]*t[12]*t[19]*t[22]*t[52]*t[32]+t[208]+1.8*t[10]*t[54]*t[5]*t[24]*t[32]+t[207]+1.8*t[10]*t[19]*t[28]*t[27]*t[24]*t[32]+t[206]+t[205]+.09601599999999999*t[17]*t[10]*t[34]*t[5]*t[27]*t[32]+.06618174748244086*grada*t[29]*t[102]*t[5]*t[27]*t[32]+t[204]+0.12002*grada*t[29]*t[12]*t[19]*t[28]*t[9]*t[27]*t[32]+0.048008*t[17]*t[10]*t[34]*t[5]*t[9]*t[32]+0.12002*grada*t[29]*t[12]*t[54]*t[5]*t[9]*t[32]-0.36006*grada*t[29]*t[12]*t[5]*t[9]*t[52]*t[68]-5.4*t[10]*t[5]*t[52]*t[24]*t[68]+t[203]-0.3*t[6]*t[41]*t[56]*t[14]*t[18]-0.3375*t[10]*t[77]*t[22]*t[56]*t[14]*t[18]+0.6*t[6]*t[55]*t[19]*t[14]*t[18]-0.45*t[10]*t[101]*t[19]*t[14]*t[18]+0.45*t[10]*t[44]*t[41]*t[28]*t[14]*t[18]+0.225*t[10]*t[55]*t[41]*t[22]*t[14]*t[18]-2.4*t[6]*t[54]*t[5]*t[32]*t[18]+t[202]+1.8*t[10]*t[44]*t[19]*t[27]*t[32]*t[18]-2.4*t[6]*t[19]*t[28]*t[27]*t[32]*t[18]-0.9*t[10]*t[41]*t[22]*t[28]*t[27]*t[32]*t[18]-0.45*t[10]*t[41]*t[56]*t[9]*t[32]*t[18]+0.9*t[10]*t[55]*t[19]*t[9]*t[32]*t[18]+0.49628039271952*t[10]*t[34]*t[19]*t[28]*t[32]*t[18]+0.9*t[10]*t[54]*t[19]*t[22]*t[32]*t[18]+7.2*t[6]*t[5]*t[52]*t[68]*t[18]-2.7*t[10]*t[19]*t[22]*t[52]*t[68]*t[18]-2.97768235631712*t[10]*t[34]*t[5]*t[27]*t[68]*t[18]-5.4*t[10]*t[19]*t[28]*t[9]*t[27]*t[68]*t[18]-5.4*t[10]*t[54]*t[5]*t[9]*t[68]*t[18]+21.6*t[10]*t[5]*t[9]*t[52]*t[180]*t[18];
909 d4fdrarbrbgb = .005000833333333333*gradb*t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*gradb*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[200]+t[232]-.01000166666666667*gradb*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*gradb*t[29]*t[12]*t[101]*t[19]*t[23]+.02667111111111111*gradb*t[29]*t[34]*t[19]*t[28]*t[23]-0.030005*gradb*t[29]*t[12]*t[44]*t[41]*t[28]*t[23]+t[231]+.01333555555555555*gradb*t[29]*t[34]*t[19]*t[22]*t[23]-0.0150025*gradb*t[29]*t[12]*t[55]*t[41]*t[22]*t[23]+.06401066666666666*t[17]*t[26]*t[34]*t[5]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[54]*t[5]*t[14]+t[199]+t[230]+0.225*t[11]*t[41]*t[56]*t[24]*t[14]-2.8*t[53]*t[5]*t[24]*t[14]+t[229]-0.45*t[11]*t[55]*t[19]*t[24]*t[14]+1.2*t[26]*t[19]*t[28]*t[24]*t[14]+t[228]+t[227]-.05334222222222221*gradb*t[29]*t[34]*t[5]*t[27]*t[14]-0.06001*gradb*t[29]*t[12]*t[44]*t[19]*t[27]*t[14]+.02000333333333333*gradb*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]+0.030005*gradb*t[29]*t[12]*t[41]*t[22]*t[28]*t[27]*t[14]+t[226]+0.0150025*gradb*t[29]*t[12]*t[41]*t[56]*t[9]*t[14]-.02667111111111111*gradb*t[29]*t[34]*t[5]*t[9]*t[14]-0.030005*gradb*t[29]*t[12]*t[55]*t[19]*t[9]*t[14]+t[225]-0.024004*t[17]*t[11]*t[34]*t[19]*t[28]*t[14]-.01654543687061022*gradb*t[29]*t[102]*t[19]*t[28]*t[14]-0.012002*t[17]*t[11]*t[34]*t[19]*t[22]*t[14]-0.030005*gradb*t[29]*t[12]*t[54]*t[19]*t[22]*t[14]-.04000666666666666*gradb*t[29]*t[7]*t[5]*t[52]*t[32]+0.06001*gradb*t[29]*t[12]*t[19]*t[22]*t[52]*t[32]+t[224]+1.8*t[11]*t[54]*t[5]*t[24]*t[32]-4.8*t[26]*t[5]*t[27]*t[24]*t[32]+1.8*t[11]*t[19]*t[28]*t[27]*t[24]*t[32]+t[223]+t[222]+t[221]+.09601599999999999*t[17]*t[11]*t[34]*t[5]*t[27]*t[32]+.06618174748244086*gradb*t[29]*t[102]*t[5]*t[27]*t[32]+t[220]+0.12002*gradb*t[29]*t[12]*t[19]*t[28]*t[9]*t[27]*t[32]+0.048008*t[17]*t[11]*t[34]*t[5]*t[9]*t[32]+0.12002*gradb*t[29]*t[12]*t[54]*t[5]*t[9]*t[32]-0.36006*gradb*t[29]*t[12]*t[5]*t[9]*t[52]*t[68]-5.4*t[11]*t[5]*t[52]*t[24]*t[68]+t[219]-0.3375*t[11]*t[77]*t[22]*t[56]*t[14]*t[18]+1.2*t[26]*t[44]*t[19]*t[14]*t[18]-0.45*t[11]*t[101]*t[19]*t[14]*t[18]+0.45*t[11]*t[44]*t[41]*t[28]*t[14]*t[18]-0.6*t[26]*t[41]*t[22]*t[28]*t[14]*t[18]-1.4*t[53]*t[19]*t[22]*t[14]*t[18]+0.225*t[11]*t[55]*t[41]*t[22]*t[14]*t[18]-1.323414380585386*t[26]*t[34]*t[5]*t[32]*t[18]+t[218]+1.8*t[11]*t[44]*t[19]*t[27]*t[32]*t[18]-0.9*t[11]*t[41]*t[22]*t[28]*t[27]*t[32]*t[18]-2.4*t[26]*t[19]*t[22]*t[27]*t[32]*t[18]-0.45*t[11]*t[41]*t[56]*t[9]*t[32]*t[18]+5.6*t[53]*t[5]*t[9]*t[32]*t[18]+0.9*t[11]*t[55]*t[19]*t[9]*t[32]*t[18]-2.4*t[26]*t[19]*t[28]*t[9]*t[32]*t[18]+0.49628039271952*t[11]*t[34]*t[19]*t[28]*t[32]*t[18]+0.9*t[11]*t[54]*t[19]*t[22]*t[32]*t[18]-2.7*t[11]*t[19]*t[22]*t[52]*t[68]*t[18]-2.97768235631712*t[11]*t[34]*t[5]*t[27]*t[68]*t[18]+14.4*t[26]*t[5]*t[9]*t[27]*t[68]*t[18]-5.4*t[11]*t[19]*t[28]*t[9]*t[27]*t[68]*t[18]-5.4*t[11]*t[54]*t[5]*t[9]*t[68]*t[18]+21.6*t[11]*t[5]*t[9]*t[52]*t[180]*t[18];
910 d4fdrarbrbab = .005000833333333333*t[29]*t[7]*t[41]*t[56]*t[23]+0.02250375*t[29]*t[12]*t[77]*t[22]*t[56]*t[23]+t[201]+t[238]-.01000166666666667*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*t[29]*t[12]*t[101]*t[19]*t[23]+.02667111111111111*t[29]*t[34]*t[19]*t[28]*t[23]-0.030005*t[29]*t[12]*t[44]*t[41]*t[28]*t[23]+t[237]+.01333555555555555*t[29]*t[34]*t[19]*t[22]*t[23]-0.0150025*t[29]*t[12]*t[55]*t[41]*t[22]*t[23]+.02000333333333333*t[29]*t[7]*t[54]*t[5]*t[14]+t[236]-.05334222222222221*t[29]*t[34]*t[5]*t[27]*t[14]-0.06001*t[29]*t[12]*t[44]*t[19]*t[27]*t[14]+.02000333333333333*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]+0.030005*t[29]*t[12]*t[41]*t[22]*t[28]*t[27]*t[14]+t[235]+0.0150025*t[29]*t[12]*t[41]*t[56]*t[9]*t[14]-.02667111111111111*t[29]*t[34]*t[5]*t[9]*t[14]-0.030005*t[29]*t[12]*t[55]*t[19]*t[9]*t[14]+t[234]-.01654543687061022*t[29]*t[102]*t[19]*t[28]*t[14]-0.030005*t[29]*t[12]*t[54]*t[19]*t[22]*t[14]-.04000666666666666*t[29]*t[7]*t[5]*t[52]*t[32]+0.06001*t[29]*t[12]*t[19]*t[22]*t[52]*t[32]+.06618174748244086*t[29]*t[102]*t[5]*t[27]*t[32]+t[233]+0.12002*t[29]*t[12]*t[19]*t[28]*t[9]*t[27]*t[32]+0.12002*t[29]*t[12]*t[54]*t[5]*t[9]*t[32]-0.36006*t[29]*t[12]*t[5]*t[9]*t[52]*t[68];
911 d4fdrarbgaga = t[99]+t[241]+t[112]-0.030005*t[15]*t[60]*t[12]*t[44]*t[19]*t[23]+t[111]+.01000166666666667*t[15]*t[60]*t[7]*t[19]*t[28]*t[23]+t[110]+0.0150025*t[15]*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+t[109]+.01000166666666667*t[15]*t[60]*t[7]*t[19]*t[22]*t[23]-0.048008*grada*t[29]*t[6]*t[7]*t[5]*t[14]+t[240]+t[108]+.01654543687061022*t[15]*t[60]*t[102]*t[5]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[44]*t[19]*t[14]+t[107]-.02000333333333333*t[15]*t[60]*t[7]*t[5]*t[27]*t[14]+t[106]+0.030005*t[15]*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]+t[105]-.02000333333333333*t[15]*t[60]*t[7]*t[5]*t[9]*t[14]+t[104]+0.030005*t[15]*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]+.07201199999999999*grada*t[29]*t[6]*t[12]*t[19]*t[28]*t[14]+0.018003*grada*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]+0.0270045*grada*t[29]*t[10]*t[12]*t[41]*t[22]*t[28]*t[14]+0.018003*grada*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]+t[239]+.05956357273419678*grada*t[29]*t[10]*t[102]*t[5]*t[32]-4.32*t[113]*t[5]*t[24]*t[32]+0.81*t[59]*t[19]*t[28]*t[24]*t[32]+0.81*t[59]*t[19]*t[22]*t[24]*t[32]-0.288048*grada*t[29]*t[6]*t[12]*t[5]*t[27]*t[32]-.07201199999999999*grada*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]+t[103]-0.12002*t[15]*t[60]*t[12]*t[5]*t[9]*t[27]*t[32]+0.108018*grada*t[29]*t[10]*t[12]*t[19]*t[22]*t[27]*t[32]-.07201199999999999*grada*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]+0.108018*grada*t[29]*t[10]*t[12]*t[19]*t[28]*t[9]*t[32]-4.86*t[59]*t[5]*t[27]*t[24]*t[68]-4.86*t[59]*t[5]*t[9]*t[24]*t[68]-0.648108*grada*t[29]*t[10]*t[12]*t[5]*t[9]*t[27]*t[68]+0.81*t[59]*t[44]*t[19]*t[32]*t[18]-2.16*t[113]*t[19]*t[28]*t[32]*t[18]-0.405*t[59]*t[41]*t[22]*t[28]*t[32]*t[18]-1.339957060342704*t[59]*t[34]*t[5]*t[68]*t[18]+12.96*t[113]*t[5]*t[27]*t[68]*t[18]-2.43*t[59]*t[19]*t[22]*t[27]*t[68]*t[18]-2.43*t[59]*t[19]*t[28]*t[9]*t[68]*t[18]+19.44*t[59]*t[5]*t[9]*t[27]*t[180]*t[18];
912 d4fdrarbgagb = t[245]-0.030005*grada*gradb*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*grada*gradb*t[60]*t[7]*t[19]*t[28]*t[23]+0.0150025*grada*gradb*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+.01000166666666667*grada*gradb*t[60]*t[7]*t[19]*t[22]*t[23]-0.024004*grada*t[29]*t[26]*t[7]*t[5]*t[14]-0.024004*gradb*t[29]*t[6]*t[7]*t[5]*t[14]+t[244]+t[243]+.01654543687061022*grada*gradb*t[60]*t[102]*t[5]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[44]*t[19]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[44]*t[19]*t[14]-.02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[27]*t[14]+0.030005*grada*gradb*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]-.02000333333333333*grada*gradb*t[60]*t[7]*t[5]*t[9]*t[14]+0.030005*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]+0.036006*gradb*t[29]*t[6]*t[12]*t[19]*t[28]*t[14]+.009001499999999999*grada*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]+.009001499999999999*gradb*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]+0.01350225*grada*t[29]*t[11]*t[12]*t[41]*t[22]*t[28]*t[14]+0.01350225*gradb*t[29]*t[10]*t[12]*t[41]*t[22]*t[28]*t[14]+0.036006*grada*t[29]*t[26]*t[12]*t[19]*t[22]*t[14]+.009001499999999999*grada*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]+.009001499999999999*gradb*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]+t[242]+.02978178636709839*grada*t[29]*t[11]*t[102]*t[5]*t[32]+.02978178636709839*gradb*t[29]*t[10]*t[102]*t[5]*t[32]-2.16*t[6]*t[11]*t[5]*t[24]*t[32]-2.16*t[10]*t[26]*t[5]*t[24]*t[32]+0.81*t[10]*t[11]*t[19]*t[28]*t[24]*t[32]+0.81*t[10]*t[11]*t[19]*t[22]*t[24]*t[32]-0.144024*gradb*t[29]*t[6]*t[12]*t[5]*t[27]*t[32]-0.036006*grada*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]-0.036006*gradb*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]-0.12002*grada*gradb*t[60]*t[12]*t[5]*t[9]*t[27]*t[32]+0.054009*grada*t[29]*t[11]*t[12]*t[19]*t[22]*t[27]*t[32]+0.054009*gradb*t[29]*t[10]*t[12]*t[19]*t[22]*t[27]*t[32]-0.144024*grada*t[29]*t[26]*t[12]*t[5]*t[9]*t[32]-0.036006*grada*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]-0.036006*gradb*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]+0.054009*grada*t[29]*t[11]*t[12]*t[19]*t[28]*t[9]*t[32]+0.054009*gradb*t[29]*t[10]*t[12]*t[19]*t[28]*t[9]*t[32]-4.86*t[10]*t[11]*t[5]*t[27]*t[24]*t[68]-4.86*t[10]*t[11]*t[5]*t[9]*t[24]*t[68]-0.324054*grada*t[29]*t[11]*t[12]*t[5]*t[9]*t[27]*t[68]-0.324054*gradb*t[29]*t[10]*t[12]*t[5]*t[9]*t[27]*t[68]+2.88*t[6]*t[26]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[44]*t[19]*t[32]*t[18]-1.08*t[6]*t[11]*t[19]*t[28]*t[32]*t[18]-0.405*t[10]*t[11]*t[41]*t[22]*t[28]*t[32]*t[18]-1.08*t[10]*t[26]*t[19]*t[22]*t[32]*t[18]-1.339957060342704*t[10]*t[11]*t[34]*t[5]*t[68]*t[18]+6.48*t[6]*t[11]*t[5]*t[27]*t[68]*t[18]-2.43*t[10]*t[11]*t[19]*t[22]*t[27]*t[68]*t[18]+6.48*t[10]*t[26]*t[5]*t[9]*t[68]*t[18]-2.43*t[10]*t[11]*t[19]*t[28]*t[9]*t[68]*t[18]+19.44*t[10]*t[11]*t[5]*t[9]*t[27]*t[180]*t[18];
913 d4fdrarbgaab = t[247]-0.030005*grada*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*grada*t[60]*t[7]*t[19]*t[28]*t[23]+0.0150025*grada*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+.01000166666666667*grada*t[60]*t[7]*t[19]*t[22]*t[23]-0.024004*t[29]*t[6]*t[7]*t[5]*t[14]+t[246]+.01654543687061022*grada*t[60]*t[102]*t[5]*t[14]-0.0270045*t[29]*t[10]*t[12]*t[44]*t[19]*t[14]-.02000333333333333*grada*t[60]*t[7]*t[5]*t[27]*t[14]+0.030005*grada*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]-.02000333333333333*grada*t[60]*t[7]*t[5]*t[9]*t[14]+0.030005*grada*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]+0.036006*t[29]*t[6]*t[12]*t[19]*t[28]*t[14]+.009001499999999999*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]+0.01350225*t[29]*t[10]*t[12]*t[41]*t[22]*t[28]*t[14]+.009001499999999999*t[29]*t[10]*t[7]*t[19]*t[22]*t[14]+.02978178636709839*t[29]*t[10]*t[102]*t[5]*t[32]-0.144024*t[29]*t[6]*t[12]*t[5]*t[27]*t[32]-0.036006*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]-0.12002*grada*t[60]*t[12]*t[5]*t[9]*t[27]*t[32]+0.054009*t[29]*t[10]*t[12]*t[19]*t[22]*t[27]*t[32]-0.036006*t[29]*t[10]*t[7]*t[5]*t[9]*t[32]+0.054009*t[29]*t[10]*t[12]*t[19]*t[28]*t[9]*t[32]-0.324054*t[29]*t[10]*t[12]*t[5]*t[9]*t[27]*t[68];
914 d4fdrarbgbgb = t[99]+t[250]+t[112]-0.030005*t[16]*t[60]*t[12]*t[44]*t[19]*t[23]+t[111]+.01000166666666667*t[16]*t[60]*t[7]*t[19]*t[28]*t[23]+t[110]+0.0150025*t[16]*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+t[109]+.01000166666666667*t[16]*t[60]*t[7]*t[19]*t[22]*t[23]-0.048008*gradb*t[29]*t[26]*t[7]*t[5]*t[14]+t[249]+t[108]+.01654543687061022*t[16]*t[60]*t[102]*t[5]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[44]*t[19]*t[14]+t[107]-.02000333333333333*t[16]*t[60]*t[7]*t[5]*t[27]*t[14]+t[106]+0.030005*t[16]*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]+t[105]-.02000333333333333*t[16]*t[60]*t[7]*t[5]*t[9]*t[14]+t[104]+0.030005*t[16]*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]+0.018003*gradb*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]+0.0270045*gradb*t[29]*t[11]*t[12]*t[41]*t[22]*t[28]*t[14]+.07201199999999999*gradb*t[29]*t[26]*t[12]*t[19]*t[22]*t[14]+0.018003*gradb*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]+t[248]+.05956357273419678*gradb*t[29]*t[11]*t[102]*t[5]*t[32]-4.32*t[156]*t[5]*t[24]*t[32]+0.81*t[63]*t[19]*t[28]*t[24]*t[32]+0.81*t[63]*t[19]*t[22]*t[24]*t[32]-.07201199999999999*gradb*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]+t[103]-0.12002*t[16]*t[60]*t[12]*t[5]*t[9]*t[27]*t[32]+0.108018*gradb*t[29]*t[11]*t[12]*t[19]*t[22]*t[27]*t[32]-0.288048*gradb*t[29]*t[26]*t[12]*t[5]*t[9]*t[32]-.07201199999999999*gradb*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]+0.108018*gradb*t[29]*t[11]*t[12]*t[19]*t[28]*t[9]*t[32]-4.86*t[63]*t[5]*t[27]*t[24]*t[68]-4.86*t[63]*t[5]*t[9]*t[24]*t[68]-0.648108*gradb*t[29]*t[11]*t[12]*t[5]*t[9]*t[27]*t[68]+0.81*t[63]*t[44]*t[19]*t[32]*t[18]-0.405*t[63]*t[41]*t[22]*t[28]*t[32]*t[18]-2.16*t[156]*t[19]*t[22]*t[32]*t[18]-1.339957060342704*t[63]*t[34]*t[5]*t[68]*t[18]-2.43*t[63]*t[19]*t[22]*t[27]*t[68]*t[18]+12.96*t[156]*t[5]*t[9]*t[68]*t[18]-2.43*t[63]*t[19]*t[28]*t[9]*t[68]*t[18]+19.44*t[63]*t[5]*t[9]*t[27]*t[180]*t[18];
915 d4fdrarbgbab = t[252]-0.030005*gradb*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*gradb*t[60]*t[7]*t[19]*t[28]*t[23]+0.0150025*gradb*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+.01000166666666667*gradb*t[60]*t[7]*t[19]*t[22]*t[23]-0.024004*t[29]*t[26]*t[7]*t[5]*t[14]+t[251]+.01654543687061022*gradb*t[60]*t[102]*t[5]*t[14]-0.0270045*t[29]*t[11]*t[12]*t[44]*t[19]*t[14]-.02000333333333333*gradb*t[60]*t[7]*t[5]*t[27]*t[14]+0.030005*gradb*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]-.02000333333333333*gradb*t[60]*t[7]*t[5]*t[9]*t[14]+0.030005*gradb*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]+.009001499999999999*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]+0.01350225*t[29]*t[11]*t[12]*t[41]*t[22]*t[28]*t[14]+0.036006*t[29]*t[26]*t[12]*t[19]*t[22]*t[14]+.009001499999999999*t[29]*t[11]*t[7]*t[19]*t[22]*t[14]+.02978178636709839*t[29]*t[11]*t[102]*t[5]*t[32]-0.036006*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]-0.12002*gradb*t[60]*t[12]*t[5]*t[9]*t[27]*t[32]+0.054009*t[29]*t[11]*t[12]*t[19]*t[22]*t[27]*t[32]-0.144024*t[29]*t[26]*t[12]*t[5]*t[9]*t[32]-0.036006*t[29]*t[11]*t[7]*t[5]*t[9]*t[32]+0.054009*t[29]*t[11]*t[12]*t[19]*t[28]*t[9]*t[32]-0.324054*t[29]*t[11]*t[12]*t[5]*t[9]*t[27]*t[68];
916 d4fdrarbabab = t[253]-0.030005*t[60]*t[12]*t[44]*t[19]*t[23]+.01000166666666667*t[60]*t[7]*t[19]*t[28]*t[23]+0.0150025*t[60]*t[12]*t[41]*t[22]*t[28]*t[23]+.01000166666666667*t[60]*t[7]*t[19]*t[22]*t[23]+.01654543687061022*t[60]*t[102]*t[5]*t[14]-.02000333333333333*t[60]*t[7]*t[5]*t[27]*t[14]+0.030005*t[60]*t[12]*t[19]*t[22]*t[27]*t[14]-.02000333333333333*t[60]*t[7]*t[5]*t[9]*t[14]+0.030005*t[60]*t[12]*t[19]*t[28]*t[9]*t[14]-0.12002*t[60]*t[12]*t[5]*t[9]*t[27]*t[32];
917 d4fdragagaga = 17.496*t[18]*t[5]*t[68]/pow(rhoa,5.0)-.8749458000000001*grada*t[29]*t[59]*t[12]*t[5]*t[9]*t[68]-0.324054*t[15]*t[60]*t[10]*t[12]*t[5]*t[9]*t[32]+0.324054*t[29]*t[10]*t[12]*t[5]*t[9]*t[32]-0.7777296*grada*t[29]*t[113]*t[12]*t[5]*t[32]+0.1458243*grada*t[29]*t[59]*t[12]*t[19]*t[22]*t[32]+t[259]+t[258]+t[257]+t[256]+t[255]+t[254]-0.090015*grada*t[60]*t[12]*t[19]*t[22]*t[23]+0.090015*t[165]*t[166]*t[12]*t[19]*t[22]*t[23]-2.187*t[164]*t[19]*t[22]*t[68]*t[18]+17.496*t[164]*t[5]*t[9]*t[180]*t[18]+0.18003*grada*t[60]*t[12]*t[5]*t[9]*t[14]-0.18003*t[165]*t[166]*t[12]*t[5]*t[9]*t[14]-0.216036*t[15]*t[60]*t[6]*t[12]*t[5]*t[14]+0.216036*t[29]*t[6]*t[12]*t[5]*t[14]+0.0810135*t[15]*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.0810135*t[29]*t[10]*t[12]*t[19]*t[22]*t[14];
918 d4fdragagagb = t[136]+t[265]+t[135]+0.090015*t[15]*gradb*t[166]*t[12]*t[19]*t[22]*t[23]-0.144024*grada*gradb*t[60]*t[6]*t[12]*t[5]*t[14]+t[134]+t[264]+t[263]+t[133]-0.18003*t[15]*gradb*t[166]*t[12]*t[5]*t[9]*t[14]+t[132]+0.0270045*t[15]*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]+0.054009*grada*gradb*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.2592432*grada*t[29]*t[6]*t[11]*t[12]*t[5]*t[32]-0.2592432*gradb*t[29]*t[113]*t[12]*t[5]*t[32]+t[262]+t[261]+t[131]-0.108018*t[15]*t[60]*t[11]*t[12]*t[5]*t[9]*t[32]-0.216036*grada*gradb*t[60]*t[10]*t[12]*t[5]*t[9]*t[32]+0.0972162*grada*t[29]*t[10]*t[11]*t[12]*t[19]*t[22]*t[32]+0.0486081*gradb*t[29]*t[59]*t[12]*t[19]*t[22]*t[32]+t[260]-.5832972000000001*grada*t[29]*t[10]*t[11]*t[12]*t[5]*t[9]*t[68]-0.2916486*gradb*t[29]*t[59]*t[12]*t[5]*t[9]*t[68]+11.664*t[113]*t[11]*t[5]*t[68]*t[18]-2.187*t[59]*t[11]*t[19]*t[22]*t[68]*t[18]+17.496*t[59]*t[11]*t[5]*t[9]*t[180]*t[18];
919 d4fdragagaab = t[139]+t[268]+t[138]+0.090015*t[15]*t[166]*t[12]*t[19]*t[22]*t[23]-0.144024*grada*t[60]*t[6]*t[12]*t[5]*t[14]+t[267]+t[137]-0.18003*t[15]*t[166]*t[12]*t[5]*t[9]*t[14]+0.054009*grada*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.2592432*t[29]*t[113]*t[12]*t[5]*t[32]+t[266]-0.216036*grada*t[60]*t[10]*t[12]*t[5]*t[9]*t[32]+0.0486081*t[29]*t[59]*t[12]*t[19]*t[22]*t[32]-0.2916486*t[29]*t[59]*t[12]*t[5]*t[9]*t[68];
920 d4fdragagbgb = t[127]+t[274]+t[126]+0.090015*grada*t[16]*t[166]*t[12]*t[19]*t[22]*t[23]+t[125]-.07201199999999999*t[16]*t[60]*t[6]*t[12]*t[5]*t[14]+t[273]+t[124]+t[272]+t[123]-0.18003*grada*t[16]*t[166]*t[12]*t[5]*t[9]*t[14]+0.054009*grada*gradb*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]+t[122]+0.0270045*t[16]*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.2592432*gradb*t[29]*t[6]*t[11]*t[12]*t[5]*t[32]+t[271]+t[270]-0.216036*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[9]*t[32]+t[121]-0.108018*t[16]*t[60]*t[10]*t[12]*t[5]*t[9]*t[32]+0.0972162*gradb*t[29]*t[10]*t[11]*t[12]*t[19]*t[22]*t[32]+0.0486081*grada*t[29]*t[63]*t[12]*t[19]*t[22]*t[32]+t[269]-.5832972000000001*gradb*t[29]*t[10]*t[11]*t[12]*t[5]*t[9]*t[68]-0.2916486*grada*t[29]*t[63]*t[12]*t[5]*t[9]*t[68]+5.832*t[6]*t[63]*t[5]*t[68]*t[18]-2.187*t[10]*t[63]*t[19]*t[22]*t[68]*t[18]+17.496*t[10]*t[63]*t[5]*t[9]*t[180]*t[18];
921 d4fdragagbab = t[278]+0.090015*grada*gradb*t[166]*t[12]*t[19]*t[22]*t[23]-.07201199999999999*gradb*t[60]*t[6]*t[12]*t[5]*t[14]+t[277]+t[276]-0.18003*grada*gradb*t[166]*t[12]*t[5]*t[9]*t[14]+0.0270045*grada*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]+0.0270045*gradb*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.1296216*t[29]*t[6]*t[11]*t[12]*t[5]*t[32]+t[275]-0.108018*grada*t[60]*t[11]*t[12]*t[5]*t[9]*t[32]-0.108018*gradb*t[60]*t[10]*t[12]*t[5]*t[9]*t[32]+0.0486081*t[29]*t[10]*t[11]*t[12]*t[19]*t[22]*t[32]-0.2916486*t[29]*t[10]*t[11]*t[12]*t[5]*t[9]*t[68];
922 d4fdragaabab = t[280]+0.090015*grada*t[166]*t[12]*t[19]*t[22]*t[23]-.07201199999999999*t[60]*t[6]*t[12]*t[5]*t[14]+t[279]-0.18003*grada*t[166]*t[12]*t[5]*t[9]*t[14]+0.0270045*t[60]*t[10]*t[12]*t[19]*t[22]*t[14]-0.108018*t[60]*t[10]*t[12]*t[5]*t[9]*t[32];
923 d4fdragbgbgb = t[286]+t[285]-0.090015*gradb*t[60]*t[12]*t[19]*t[22]*t[23]+0.090015*t[173]*t[166]*t[12]*t[19]*t[22]*t[23]+t[284]+t[283]+0.18003*gradb*t[60]*t[12]*t[5]*t[9]*t[14]-0.18003*t[173]*t[166]*t[12]*t[5]*t[9]*t[14]-0.0810135*t[29]*t[11]*t[12]*t[19]*t[22]*t[14]+0.0810135*t[16]*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]+t[282]+0.324054*t[29]*t[11]*t[12]*t[5]*t[9]*t[32]-0.324054*t[16]*t[60]*t[11]*t[12]*t[5]*t[9]*t[32]+0.1458243*gradb*t[29]*t[63]*t[12]*t[19]*t[22]*t[32]+t[281]-.8749458000000001*gradb*t[29]*t[63]*t[12]*t[5]*t[9]*t[68]-2.187*t[172]*t[19]*t[22]*t[68]*t[18]+17.496*t[172]*t[5]*t[9]*t[180]*t[18];
924 d4fdragbgbab = t[139]+t[289]+t[138]+0.090015*t[16]*t[166]*t[12]*t[19]*t[22]*t[23]+t[288]+t[137]-0.18003*t[16]*t[166]*t[12]*t[5]*t[9]*t[14]+0.054009*gradb*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]+t[287]-0.216036*gradb*t[60]*t[11]*t[12]*t[5]*t[9]*t[32]+0.0486081*t[29]*t[63]*t[12]*t[19]*t[22]*t[32]-0.2916486*t[29]*t[63]*t[12]*t[5]*t[9]*t[68];
925 d4fdragbabab = t[291]+0.090015*gradb*t[166]*t[12]*t[19]*t[22]*t[23]+t[290]-0.18003*gradb*t[166]*t[12]*t[5]*t[9]*t[14]+0.0270045*t[60]*t[11]*t[12]*t[19]*t[22]*t[14]-0.108018*t[60]*t[11]*t[12]*t[5]*t[9]*t[32];
926 d4fdraababab = t[292]+0.090015*t[166]*t[12]*t[19]*t[22]*t[23]-0.18003*t[166]*t[12]*t[5]*t[9]*t[14];
927 d4fdrbrbrbrb = 2.25*t[18]*t[23]*t[55]*t[56]*t[77]-1.5*t[14]*t[144]*t[18]*t[27]*t[77]+1.5*t[144]*t[23]*t[24]*t[77]+0.5*t[18]*t[19]*t[23]*(-192.0*t[1]*t[73]+t[185]+t[184])-36.0*t[18]*t[5]*t[52]*t[54]*t[68]-24.0*t[140]*t[24]*t[5]*t[68]-12.0*t[140]*t[18]*t[19]*t[28]*t[68]+1.5*t[14]*t[18]*t[41]*t[54]*t[56]-3.0*t[18]*t[32]*t[41]*t[52]*t[56]+3.0*t[14]*t[24]*t[27]*t[41]*t[56]-0.75*t[18]*t[23]*t[41]*pow(t[55],2.0)-3.0*t[14]*t[18]*t[19]*t[54]*t[55]+6.0*t[18]*t[19]*t[32]*t[52]*t[55]+3.0*t[14]*t[18]*t[27]*t[28]*t[41]*t[55]-3.0*t[23]*t[24]*t[28]*t[41]*t[55]-6.0*t[14]*t[19]*t[24]*t[27]*t[55]+6.0*t[18]*t[32]*t[5]*pow(t[54],2.0)+24.0*t[24]*t[27]*t[32]*t[5]*t[54]+12.0*t[18]*t[19]*t[27]*t[28]*t[32]*t[54]-6.0*t[14]*t[19]*t[24]*t[28]*t[54]+12.0*t[19]*t[24]*t[28]*t[32]*t[52]+8.0*t[142]*t[18]*t[27]*t[32]*t[5]+24.0*t[18]*t[180]*pow(t[27],4.0)*t[5]-4.0*t[14]*t[142]*t[24]*t[5]-1.0*t[14]*t[18]*(t[183]+40.44444444444444*gradb*t[293])*t[5]-1.0*t[143]*t[18]*t[23]*t[28]*t[41]+.3200533333333333*t[17]*t[34]*t[5]*t[52]*t[32]-0.9375*t[18]*t[186]*t[23]*pow(t[28],4.0)-2.0*t[14]*t[142]*t[18]*t[19]*t[28]-2.0*t[14]*t[143]*t[18]*t[19]*t[27]+2.0*t[143]*t[19]*t[23]*t[24]-.04000666666666666*t[17]*t[34]*t[41]*t[56]*t[23]-.1244651851851852*t[17]*t[70]*t[19]*t[28]*t[23]+.08001333333333333*t[17]*t[34]*t[55]*t[19]*t[23]+t[187]-.1600266666666667*t[17]*t[34]*t[54]*t[5]*t[14]+.2489303703703704*t[17]*t[70]*t[5]*t[27]*t[14]-.1600266666666667*t[17]*t[34]*t[19]*t[28]*t[27]*t[14];
928 d4fdrbrbrbga = 0.02250375*grada*t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*grada*t[29]*t[7]*t[41]*t[56]*t[23]+t[198]-0.030005*grada*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*grada*t[29]*t[12]*t[143]*t[19]*t[23]+.04000666666666666*grada*t[29]*t[34]*t[19]*t[28]*t[23]-0.0450075*grada*t[29]*t[12]*t[55]*t[41]*t[28]*t[23]+.06000999999999999*grada*t[29]*t[7]*t[54]*t[5]*t[14]+t[197]-0.06001*grada*t[29]*t[12]*t[142]*t[5]*t[14]+0.675*t[10]*t[41]*t[56]*t[24]*t[14]-1.35*t[10]*t[55]*t[19]*t[24]*t[14]+0.0450075*grada*t[29]*t[12]*t[41]*t[56]*t[27]*t[14]-.08001333333333333*grada*t[29]*t[34]*t[5]*t[27]*t[14]-0.090015*grada*t[29]*t[12]*t[55]*t[19]*t[27]*t[14]+.06000999999999999*grada*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]-0.036006*t[17]*t[10]*t[34]*t[19]*t[28]*t[14]-0.090015*grada*t[29]*t[12]*t[54]*t[19]*t[28]*t[14]-0.12002*grada*t[29]*t[7]*t[5]*t[52]*t[32]+0.18003*grada*t[29]*t[12]*t[19]*t[28]*t[52]*t[32]+5.4*t[10]*t[54]*t[5]*t[24]*t[32]+5.4*t[10]*t[19]*t[28]*t[27]*t[24]*t[32]+0.144024*t[17]*t[10]*t[34]*t[5]*t[27]*t[32]+0.36006*grada*t[29]*t[12]*t[54]*t[5]*t[27]*t[32]-0.36006*grada*t[29]*t[12]*t[5]*t[140]*t[68]-16.2*t[10]*t[5]*t[52]*t[24]*t[68]-0.3375*t[10]*t[77]*t[144]*t[14]*t[18]-0.45*t[10]*t[143]*t[19]*t[14]*t[18]+0.675*t[10]*t[55]*t[41]*t[28]*t[14]*t[18]+1.8*t[10]*t[142]*t[5]*t[32]*t[18]-1.35*t[10]*t[41]*t[56]*t[27]*t[32]*t[18]+2.7*t[10]*t[55]*t[19]*t[27]*t[32]*t[18]+2.7*t[10]*t[54]*t[19]*t[28]*t[32]*t[18]-8.100000000000001*t[10]*t[19]*t[28]*t[52]*t[68]*t[18]-16.2*t[10]*t[54]*t[5]*t[27]*t[68]*t[18]+21.6*t[10]*t[5]*t[140]*t[180]*t[18];
929 d4fdrbrbrbgb = 0.02250375*gradb*t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*gradb*t[29]*t[7]*t[41]*t[56]*t[23]+t[200]-0.030005*gradb*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*gradb*t[29]*t[12]*t[143]*t[19]*t[23]+.04000666666666666*gradb*t[29]*t[34]*t[19]*t[28]*t[23]-0.0450075*gradb*t[29]*t[12]*t[55]*t[41]*t[28]*t[23]+.09601599999999999*t[17]*t[26]*t[34]*t[5]*t[14]+.06000999999999999*gradb*t[29]*t[7]*t[54]*t[5]*t[14]+t[199]-0.06001*gradb*t[29]*t[12]*t[142]*t[5]*t[14]+0.675*t[11]*t[41]*t[56]*t[24]*t[14]-8.4*t[53]*t[5]*t[24]*t[14]-1.35*t[11]*t[55]*t[19]*t[24]*t[14]+3.6*t[26]*t[19]*t[28]*t[24]*t[14]+0.0450075*gradb*t[29]*t[12]*t[41]*t[56]*t[27]*t[14]-.08001333333333333*gradb*t[29]*t[34]*t[5]*t[27]*t[14]-0.090015*gradb*t[29]*t[12]*t[55]*t[19]*t[27]*t[14]+.06000999999999999*gradb*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]-0.036006*t[17]*t[11]*t[34]*t[19]*t[28]*t[14]-0.090015*gradb*t[29]*t[12]*t[54]*t[19]*t[28]*t[14]-0.12002*gradb*t[29]*t[7]*t[5]*t[52]*t[32]+0.18003*gradb*t[29]*t[12]*t[19]*t[28]*t[52]*t[32]+5.4*t[11]*t[54]*t[5]*t[24]*t[32]-14.4*t[26]*t[5]*t[27]*t[24]*t[32]+5.4*t[11]*t[19]*t[28]*t[27]*t[24]*t[32]+0.144024*t[17]*t[11]*t[34]*t[5]*t[27]*t[32]+0.36006*gradb*t[29]*t[12]*t[54]*t[5]*t[27]*t[32]-0.36006*gradb*t[29]*t[12]*t[5]*t[140]*t[68]-16.2*t[11]*t[5]*t[52]*t[24]*t[68]-0.3375*t[11]*t[77]*t[144]*t[14]*t[18]-.8999999999999999*t[26]*t[41]*t[56]*t[14]*t[18]+9.333333333333334*t[141]*t[5]*t[14]*t[18]+1.8*t[26]*t[55]*t[19]*t[14]*t[18]-0.45*t[11]*t[143]*t[19]*t[14]*t[18]-4.2*t[53]*t[19]*t[28]*t[14]*t[18]+0.675*t[11]*t[55]*t[41]*t[28]*t[14]*t[18]-7.2*t[26]*t[54]*t[5]*t[32]*t[18]+1.8*t[11]*t[142]*t[5]*t[32]*t[18]-1.35*t[11]*t[41]*t[56]*t[27]*t[32]*t[18]+16.8*t[53]*t[5]*t[27]*t[32]*t[18]+2.7*t[11]*t[55]*t[19]*t[27]*t[32]*t[18]-7.2*t[26]*t[19]*t[28]*t[27]*t[32]*t[18]+2.7*t[11]*t[54]*t[19]*t[28]*t[32]*t[18]+21.6*t[26]*t[5]*t[52]*t[68]*t[18]-8.100000000000001*t[11]*t[19]*t[28]*t[52]*t[68]*t[18]-16.2*t[11]*t[54]*t[5]*t[27]*t[68]*t[18]+21.6*t[11]*t[5]*t[140]*t[180]*t[18];
930 d4fdrbrbrbab = 0.02250375*t[29]*t[12]*t[77]*t[144]*t[23]+0.0150025*t[29]*t[7]*t[41]*t[56]*t[23]+t[201]-0.030005*t[29]*t[7]*t[55]*t[19]*t[23]+0.030005*t[29]*t[12]*t[143]*t[19]*t[23]+.04000666666666666*t[29]*t[34]*t[19]*t[28]*t[23]-0.0450075*t[29]*t[12]*t[55]*t[41]*t[28]*t[23]+.06000999999999999*t[29]*t[7]*t[54]*t[5]*t[14]-0.06001*t[29]*t[12]*t[142]*t[5]*t[14]+0.0450075*t[29]*t[12]*t[41]*t[56]*t[27]*t[14]-.08001333333333333*t[29]*t[34]*t[5]*t[27]*t[14]-0.090015*t[29]*t[12]*t[55]*t[19]*t[27]*t[14]+.06000999999999999*t[29]*t[7]*t[19]*t[28]*t[27]*t[14]-0.090015*t[29]*t[12]*t[54]*t[19]*t[28]*t[14]-0.12002*t[29]*t[7]*t[5]*t[52]*t[32]+0.18003*t[29]*t[12]*t[19]*t[28]*t[52]*t[32]+0.36006*t[29]*t[12]*t[54]*t[5]*t[27]*t[32]-0.36006*t[29]*t[12]*t[5]*t[140]*t[68];
931 d4fdrbrbgaga = t[151]+0.0150025*t[15]*t[60]*t[12]*t[41]*t[56]*t[23]+t[99]+t[241]+t[150]-0.030005*t[15]*t[60]*t[12]*t[55]*t[19]*t[23]+t[149]+.02000333333333333*t[15]*t[60]*t[7]*t[19]*t[28]*t[23]+0.0270045*grada*t[29]*t[10]*t[12]*t[41]*t[56]*t[14]+t[240]+t[148]+0.06001*t[15]*t[60]*t[12]*t[54]*t[5]*t[14]-0.054009*grada*t[29]*t[10]*t[12]*t[55]*t[19]*t[14]+t[147]-.04000666666666666*t[15]*t[60]*t[7]*t[5]*t[27]*t[14]+t[146]+0.06001*t[15]*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]+0.036006*grada*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]+t[145]-0.12002*t[15]*t[60]*t[12]*t[5]*t[52]*t[32]+t[239]+0.216036*grada*t[29]*t[10]*t[12]*t[54]*t[5]*t[32]+1.62*t[59]*t[19]*t[28]*t[24]*t[32]-0.144024*grada*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]+0.216036*grada*t[29]*t[10]*t[12]*t[19]*t[28]*t[27]*t[32]-0.648108*grada*t[29]*t[10]*t[12]*t[5]*t[52]*t[68]-9.72*t[59]*t[5]*t[27]*t[24]*t[68]-0.405*t[59]*t[41]*t[56]*t[32]*t[18]+0.81*t[59]*t[55]*t[19]*t[32]*t[18]-4.86*t[59]*t[54]*t[5]*t[68]*t[18]-4.86*t[59]*t[19]*t[28]*t[27]*t[68]*t[18]+19.44*t[59]*t[5]*t[52]*t[180]*t[18];
932 d4fdrbrbgagb = 0.0150025*grada*gradb*t[60]*t[12]*t[41]*t[56]*t[23]+t[245]-0.030005*grada*gradb*t[60]*t[12]*t[55]*t[19]*t[23]+.02000333333333333*grada*gradb*t[60]*t[7]*t[19]*t[28]*t[23]+0.01350225*grada*t[29]*t[11]*t[12]*t[41]*t[56]*t[14]+0.01350225*gradb*t[29]*t[10]*t[12]*t[41]*t[56]*t[14]-0.168028*grada*t[29]*t[53]*t[12]*t[5]*t[14]-0.048008*grada*t[29]*t[26]*t[7]*t[5]*t[14]+t[244]+t[243]+0.06001*grada*gradb*t[60]*t[12]*t[54]*t[5]*t[14]-0.0270045*grada*t[29]*t[11]*t[12]*t[55]*t[19]*t[14]-0.0270045*gradb*t[29]*t[10]*t[12]*t[55]*t[19]*t[14]-.04000666666666666*grada*gradb*t[60]*t[7]*t[5]*t[27]*t[14]+0.06001*grada*gradb*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]+.07201199999999999*grada*t[29]*t[26]*t[12]*t[19]*t[28]*t[14]+0.018003*grada*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]+0.018003*gradb*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]-0.12002*grada*gradb*t[60]*t[12]*t[5]*t[52]*t[32]+t[242]+0.108018*grada*t[29]*t[11]*t[12]*t[54]*t[5]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[54]*t[5]*t[32]-4.32*t[10]*t[26]*t[5]*t[24]*t[32]+1.62*t[10]*t[11]*t[19]*t[28]*t[24]*t[32]-0.288048*grada*t[29]*t[26]*t[12]*t[5]*t[27]*t[32]-.07201199999999999*grada*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]-.07201199999999999*gradb*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]+0.108018*grada*t[29]*t[11]*t[12]*t[19]*t[28]*t[27]*t[32]+0.108018*gradb*t[29]*t[10]*t[12]*t[19]*t[28]*t[27]*t[32]-0.324054*grada*t[29]*t[11]*t[12]*t[5]*t[52]*t[68]-0.324054*gradb*t[29]*t[10]*t[12]*t[5]*t[52]*t[68]-9.72*t[10]*t[11]*t[5]*t[27]*t[24]*t[68]-0.405*t[10]*t[11]*t[41]*t[56]*t[32]*t[18]+5.040000000000001*t[10]*t[53]*t[5]*t[32]*t[18]+0.81*t[10]*t[11]*t[55]*t[19]*t[32]*t[18]-2.16*t[10]*t[26]*t[19]*t[28]*t[32]*t[18]-4.86*t[10]*t[11]*t[54]*t[5]*t[68]*t[18]+12.96*t[10]*t[26]*t[5]*t[27]*t[68]*t[18]-4.86*t[10]*t[11]*t[19]*t[28]*t[27]*t[68]*t[18]+19.44*t[10]*t[11]*t[5]*t[52]*t[180]*t[18];
933 d4fdrbrbgaab = 0.0150025*grada*t[60]*t[12]*t[41]*t[56]*t[23]+t[247]-0.030005*grada*t[60]*t[12]*t[55]*t[19]*t[23]+.02000333333333333*grada*t[60]*t[7]*t[19]*t[28]*t[23]+0.01350225*t[29]*t[10]*t[12]*t[41]*t[56]*t[14]+t[246]+0.06001*grada*t[60]*t[12]*t[54]*t[5]*t[14]-0.0270045*t[29]*t[10]*t[12]*t[55]*t[19]*t[14]-.04000666666666666*grada*t[60]*t[7]*t[5]*t[27]*t[14]+0.06001*grada*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]+0.018003*t[29]*t[10]*t[7]*t[19]*t[28]*t[14]-0.12002*grada*t[60]*t[12]*t[5]*t[52]*t[32]+0.108018*t[29]*t[10]*t[12]*t[54]*t[5]*t[32]-.07201199999999999*t[29]*t[10]*t[7]*t[5]*t[27]*t[32]+0.108018*t[29]*t[10]*t[12]*t[19]*t[28]*t[27]*t[32]-0.324054*t[29]*t[10]*t[12]*t[5]*t[52]*t[68];
934 d4fdrbrbgbgb = 15.84*t[18]*t[32]*t[5]/pow(rhob,4.666666666666667)+t[99]-0.648108*gradb*t[29]*t[11]*t[12]*t[5]*t[52]*t[68]-9.72*t[63]*t[5]*t[27]*t[24]*t[68]-0.12002*t[16]*t[60]*t[12]*t[5]*t[52]*t[32]+0.216036*gradb*t[29]*t[11]*t[12]*t[54]*t[5]*t[32]-0.144024*gradb*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]-.5760959999999999*gradb*t[29]*t[26]*t[12]*t[5]*t[27]*t[32]+0.216036*gradb*t[29]*t[11]*t[12]*t[19]*t[28]*t[27]*t[32]-8.64*t[156]*t[5]*t[24]*t[32]+1.62*t[63]*t[19]*t[28]*t[24]*t[32]+t[250]+t[249]+t[248]+0.0150025*t[16]*t[60]*t[12]*t[41]*t[56]*t[23]+.02000333333333333*t[16]*t[60]*t[7]*t[19]*t[28]*t[23]-0.030005*t[16]*t[60]*t[12]*t[55]*t[19]*t[23]-4.86*t[63]*t[54]*t[5]*t[68]*t[18]+25.92*t[156]*t[5]*t[27]*t[68]*t[18]-4.86*t[63]*t[19]*t[28]*t[27]*t[68]*t[18]-0.405*t[63]*t[41]*t[56]*t[32]*t[18]-4.32*t[156]*t[19]*t[28]*t[32]*t[18]+0.81*t[63]*t[55]*t[19]*t[32]*t[18]+19.44*t[63]*t[5]*t[52]*t[180]*t[18]+t[151]+t[150]+t[149]+t[148]+t[147]+t[146]+t[145]+0.0270045*gradb*t[29]*t[11]*t[12]*t[41]*t[56]*t[14]-.09601599999999999*gradb*t[29]*t[26]*t[7]*t[5]*t[14]+0.06001*t[16]*t[60]*t[12]*t[54]*t[5]*t[14]-0.336056*gradb*t[29]*t[53]*t[12]*t[5]*t[14]+0.036006*gradb*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]+0.144024*gradb*t[29]*t[26]*t[12]*t[19]*t[28]*t[14]-.04000666666666666*t[16]*t[60]*t[7]*t[5]*t[27]*t[14]+0.06001*t[16]*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]-0.054009*gradb*t[29]*t[11]*t[12]*t[55]*t[19]*t[14];
935 d4fdrbrbgbab = 0.0150025*gradb*t[60]*t[12]*t[41]*t[56]*t[23]+t[252]-0.030005*gradb*t[60]*t[12]*t[55]*t[19]*t[23]+.02000333333333333*gradb*t[60]*t[7]*t[19]*t[28]*t[23]+0.01350225*t[29]*t[11]*t[12]*t[41]*t[56]*t[14]-0.168028*t[29]*t[53]*t[12]*t[5]*t[14]-0.048008*t[29]*t[26]*t[7]*t[5]*t[14]+t[251]+0.06001*gradb*t[60]*t[12]*t[54]*t[5]*t[14]-0.0270045*t[29]*t[11]*t[12]*t[55]*t[19]*t[14]-.04000666666666666*gradb*t[60]*t[7]*t[5]*t[27]*t[14]+0.06001*gradb*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]+.07201199999999999*t[29]*t[26]*t[12]*t[19]*t[28]*t[14]+0.018003*t[29]*t[11]*t[7]*t[19]*t[28]*t[14]-0.12002*gradb*t[60]*t[12]*t[5]*t[52]*t[32]+0.108018*t[29]*t[11]*t[12]*t[54]*t[5]*t[32]-0.288048*t[29]*t[26]*t[12]*t[5]*t[27]*t[32]-.07201199999999999*t[29]*t[11]*t[7]*t[5]*t[27]*t[32]+0.108018*t[29]*t[11]*t[12]*t[19]*t[28]*t[27]*t[32]-0.324054*t[29]*t[11]*t[12]*t[5]*t[52]*t[68];
936 d4fdrbrbabab =0.0150025*t[60]*t[12]*t[41]*t[56]*t[23]+t[253]-0.030005*t[60]*t[12]*t[55]*t[19]*t[23]+.02000333333333333*t[60]*t[7]*t[19]*t[28]*t[23]+0.06001*t[60]*t[12]*t[54]*t[5]*t[14]-.04000666666666666*t[60]*t[7]*t[5]*t[27]*t[14]+0.06001*t[60]*t[12]*t[19]*t[28]*t[27]*t[14]-0.12002*t[60]*t[12]*t[5]*t[52]*t[32];
937 d4fdrbgagaga = t[259]+t[258]-0.090015*grada*t[60]*t[12]*t[19]*t[28]*t[23]+0.090015*t[165]*t[166]*t[12]*t[19]*t[28]*t[23]+t[257]+t[256]+0.18003*grada*t[60]*t[12]*t[5]*t[27]*t[14]-0.18003*t[165]*t[166]*t[12]*t[5]*t[27]*t[14]-0.0810135*t[29]*t[10]*t[12]*t[19]*t[28]*t[14]+0.0810135*t[15]*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]+t[255]+0.324054*t[29]*t[10]*t[12]*t[5]*t[27]*t[32]-0.324054*t[15]*t[60]*t[10]*t[12]*t[5]*t[27]*t[32]+0.1458243*grada*t[29]*t[59]*t[12]*t[19]*t[28]*t[32]+t[254]-.8749458000000001*grada*t[29]*t[59]*t[12]*t[5]*t[27]*t[68]-2.187*t[164]*t[19]*t[28]*t[68]*t[18]+17.496*t[164]*t[5]*t[27]*t[180]*t[18];
938 d4fdrbgagagb = t[136]+t[265]+t[161]+0.090015*t[15]*gradb*t[166]*t[12]*t[19]*t[28]*t[23]+t[160]-.07201199999999999*t[15]*t[60]*t[26]*t[12]*t[5]*t[14]+t[134]+t[264]+t[263]+t[159]-0.18003*t[15]*gradb*t[166]*t[12]*t[5]*t[27]*t[14]+t[158]+0.0270045*t[15]*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]+0.054009*grada*gradb*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]-0.2592432*grada*t[29]*t[10]*t[26]*t[12]*t[5]*t[32]+t[262]+t[261]+t[157]-0.108018*t[15]*t[60]*t[11]*t[12]*t[5]*t[27]*t[32]-0.216036*grada*gradb*t[60]*t[10]*t[12]*t[5]*t[27]*t[32]+0.0972162*grada*t[29]*t[10]*t[11]*t[12]*t[19]*t[28]*t[32]+0.0486081*gradb*t[29]*t[59]*t[12]*t[19]*t[28]*t[32]+t[260]-.5832972000000001*grada*t[29]*t[10]*t[11]*t[12]*t[5]*t[27]*t[68]-0.2916486*gradb*t[29]*t[59]*t[12]*t[5]*t[27]*t[68]+5.832*t[59]*t[26]*t[5]*t[68]*t[18]-2.187*t[59]*t[11]*t[19]*t[28]*t[68]*t[18]+17.496*t[59]*t[11]*t[5]*t[27]*t[180]*t[18];
939 d4fdrbgagaab = t[139]+t[268]+t[163]+0.090015*t[15]*t[166]*t[12]*t[19]*t[28]*t[23]+t[267]+t[162]-0.18003*t[15]*t[166]*t[12]*t[5]*t[27]*t[14]+0.054009*grada*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]+t[266]-0.216036*grada*t[60]*t[10]*t[12]*t[5]*t[27]*t[32]+0.0486081*t[29]*t[59]*t[12]*t[19]*t[28]*t[32]-0.2916486*t[29]*t[59]*t[12]*t[5]*t[27]*t[68];
940 d4fdrbgagbgb = t[127]+t[274]+t[155]+0.090015*grada*t[16]*t[166]*t[12]*t[19]*t[28]*t[23]-0.144024*grada*gradb*t[60]*t[26]*t[12]*t[5]*t[14]+t[273]+t[124]+t[272]+t[154]-0.18003*grada*t[16]*t[166]*t[12]*t[5]*t[27]*t[14]+0.054009*grada*gradb*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]+t[153]+0.0270045*t[16]*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]-0.2592432*gradb*t[29]*t[10]*t[26]*t[12]*t[5]*t[32]-0.2592432*grada*t[29]*t[156]*t[12]*t[5]*t[32]+t[271]+t[270]-0.216036*grada*gradb*t[60]*t[11]*t[12]*t[5]*t[27]*t[32]+t[152]-0.108018*t[16]*t[60]*t[10]*t[12]*t[5]*t[27]*t[32]+0.0972162*gradb*t[29]*t[10]*t[11]*t[12]*t[19]*t[28]*t[32]+0.0486081*grada*t[29]*t[63]*t[12]*t[19]*t[28]*t[32]+t[269]-.5832972000000001*gradb*t[29]*t[10]*t[11]*t[12]*t[5]*t[27]*t[68]-0.2916486*grada*t[29]*t[63]*t[12]*t[5]*t[27]*t[68]+11.664*t[10]*t[156]*t[5]*t[68]*t[18]-2.187*t[10]*t[63]*t[19]*t[28]*t[68]*t[18]+17.496*t[10]*t[63]*t[5]*t[27]*t[180]*t[18];
941 d4fdrbgagbab = t[278]+0.090015*grada*gradb*t[166]*t[12]*t[19]*t[28]*t[23]-.07201199999999999*grada*t[60]*t[26]*t[12]*t[5]*t[14]+t[277]+t[276]-0.18003*grada*gradb*t[166]*t[12]*t[5]*t[27]*t[14]+0.0270045*grada*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]+0.0270045*gradb*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]-0.1296216*t[29]*t[10]*t[26]*t[12]*t[5]*t[32]+t[275]-0.108018*grada*t[60]*t[11]*t[12]*t[5]*t[27]*t[32]-0.108018*gradb*t[60]*t[10]*t[12]*t[5]*t[27]*t[32]+0.0486081*t[29]*t[10]*t[11]*t[12]*t[19]*t[28]*t[32]-0.2916486*t[29]*t[10]*t[11]*t[12]*t[5]*t[27]*t[68];
942 d4fdrbgaabab = t[280]+0.090015*grada*t[166]*t[12]*t[19]*t[28]*t[23]+t[279]-0.18003*grada*t[166]*t[12]*t[5]*t[27]*t[14]+0.0270045*t[60]*t[10]*t[12]*t[19]*t[28]*t[14]-0.108018*t[60]*t[10]*t[12]*t[5]*t[27]*t[32];
943 d4fdrbgbgbgb = 17.496*t[18]*t[5]*t[68]/pow(rhob,5.0)-.8749458000000001*gradb*t[29]*t[63]*t[12]*t[5]*t[27]*t[68]-0.7777296*gradb*t[29]*t[156]*t[12]*t[5]*t[32]+0.1458243*gradb*t[29]*t[63]*t[12]*t[19]*t[28]*t[32]-0.324054*t[16]*t[60]*t[11]*t[12]*t[5]*t[27]*t[32]+0.324054*t[29]*t[11]*t[12]*t[5]*t[27]*t[32]+t[286]+t[285]+t[284]+t[283]+t[282]+t[281]-0.090015*gradb*t[60]*t[12]*t[19]*t[28]*t[23]+0.090015*t[173]*t[166]*t[12]*t[19]*t[28]*t[23]-2.187*t[172]*t[19]*t[28]*t[68]*t[18]+17.496*t[172]*t[5]*t[27]*t[180]*t[18]-0.216036*t[16]*t[60]*t[26]*t[12]*t[5]*t[14]+0.216036*t[29]*t[26]*t[12]*t[5]*t[14]+0.0810135*t[16]*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]-0.0810135*t[29]*t[11]*t[12]*t[19]*t[28]*t[14]+0.18003*gradb*t[60]*t[12]*t[5]*t[27]*t[14]-0.18003*t[173]*t[166]*t[12]*t[5]*t[27]*t[14];
944 d4fdrbgbgbab = t[139]+t[289]+t[163]+0.090015*t[16]*t[166]*t[12]*t[19]*t[28]*t[23]-0.144024*gradb*t[60]*t[26]*t[12]*t[5]*t[14]+t[288]+t[162]-0.18003*t[16]*t[166]*t[12]*t[5]*t[27]*t[14]+0.054009*gradb*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]-0.2592432*t[29]*t[156]*t[12]*t[5]*t[32]+t[287]-0.216036*gradb*t[60]*t[11]*t[12]*t[5]*t[27]*t[32]+0.0486081*t[29]*t[63]*t[12]*t[19]*t[28]*t[32]-0.2916486*t[29]*t[63]*t[12]*t[5]*t[27]*t[68];
945 d4fdrbgbabab = t[291]+0.090015*gradb*t[166]*t[12]*t[19]*t[28]*t[23]-.07201199999999999*t[60]*t[26]*t[12]*t[5]*t[14]+t[290]-0.18003*gradb*t[166]*t[12]*t[5]*t[27]*t[14]+0.0270045*t[60]*t[11]*t[12]*t[19]*t[28]*t[14]-0.108018*t[60]*t[11]*t[12]*t[5]*t[27]*t[32];
946 d4fdrbababab = t[292]+0.090015*t[166]*t[12]*t[19]*t[28]*t[23]-0.18003*t[166]*t[12]*t[5]*t[27]*t[14];
947 d4fdgagagaga = -0.90015*t[12]*t[23]*t[294]*t[5]*pow(grada,4.0)-1.04993496*grada*t[29]*t[164]*t[12]*t[5]*t[68]-.5832972000000001*t[15]*t[60]*t[59]*t[12]*t[5]*t[32]+.5832972000000001*t[29]*t[59]*t[12]*t[5]*t[32]+t[295]+1.08018*t[15]*t[166]*t[12]*t[5]*t[23]+15.7464*t[181]*t[5]*t[180]*t[18]+0.648108*grada*t[60]*t[10]*t[12]*t[5]*t[14]-0.648108*t[165]*t[166]*t[10]*t[12]*t[5]*t[14];
948 d4fdgagagagb = t[299]-0.90015*t[165]*gradb*t[294]*t[12]*t[5]*t[23]+t[298]-0.162027*t[165]*t[166]*t[11]*t[12]*t[5]*t[14]+t[297]-0.486081*t[15]*gradb*t[166]*t[10]*t[12]*t[5]*t[14]+t[296]-0.2916486*t[15]*t[60]*t[10]*t[11]*t[12]*t[5]*t[32]-0.2916486*grada*gradb*t[60]*t[59]*t[12]*t[5]*t[32]-.7874512200000001*grada*t[29]*t[59]*t[11]*t[12]*t[5]*t[68]-0.26248374*gradb*t[29]*t[164]*t[12]*t[5]*t[68]+15.7464*t[164]*t[11]*t[5]*t[180]*t[18];
949 d4fdgagagaab = 0.54009*grada*t[166]*t[12]*t[5]*t[23]-0.90015*t[165]*t[294]*t[12]*t[5]*t[23]+0.162027*t[60]*t[10]*t[12]*t[5]*t[14]-0.486081*t[15]*t[166]*t[10]*t[12]*t[5]*t[14]-0.2916486*grada*t[60]*t[59]*t[12]*t[5]*t[32]-0.26248374*t[29]*t[164]*t[12]*t[5]*t[68];
950 d4fdgagagbgb = t[66]+t[176]+t[169]-0.90015*t[15]*t[16]*t[294]*t[12]*t[5]*t[23]+t[175]-0.324054*t[15]*gradb*t[166]*t[11]*t[12]*t[5]*t[14]+t[168]-0.324054*grada*t[16]*t[166]*t[10]*t[12]*t[5]*t[14]-0.3888648*grada*gradb*t[60]*t[10]*t[11]*t[12]*t[5]*t[32]+t[174]-0.0972162*t[15]*t[60]*t[63]*t[12]*t[5]*t[32]+t[167]-0.0972162*t[16]*t[60]*t[59]*t[12]*t[5]*t[32]-0.52496748*gradb*t[29]*t[59]*t[11]*t[12]*t[5]*t[68]-0.52496748*grada*t[29]*t[10]*t[63]*t[12]*t[5]*t[68]+15.7464*t[59]*t[63]*t[5]*t[180]*t[18];
951 d4fdgagagbab = t[178]-0.90015*t[15]*gradb*t[294]*t[12]*t[5]*t[23]+t[177]-0.162027*t[15]*t[166]*t[11]*t[12]*t[5]*t[14]-0.324054*grada*gradb*t[166]*t[10]*t[12]*t[5]*t[14]-0.1944324*grada*t[60]*t[10]*t[11]*t[12]*t[5]*t[32]-0.0972162*gradb*t[60]*t[59]*t[12]*t[5]*t[32]-0.26248374*t[29]*t[59]*t[11]*t[12]*t[5]*t[68];
952 d4fdgagaabab = t[179]-0.90015*t[15]*t[294]*t[12]*t[5]*t[23]-0.324054*grada*t[166]*t[10]*t[12]*t[5]*t[14]-0.0972162*t[60]*t[59]*t[12]*t[5]*t[32];
953 d4fdgagbgbgb = t[299]-0.90015*grada*t[173]*t[294]*t[12]*t[5]*t[23]+t[298]-0.486081*grada*t[16]*t[166]*t[11]*t[12]*t[5]*t[14]+t[297]-0.162027*t[173]*t[166]*t[10]*t[12]*t[5]*t[14]+t[296]-0.2916486*t[16]*t[60]*t[10]*t[11]*t[12]*t[5]*t[32]-0.2916486*grada*gradb*t[60]*t[63]*t[12]*t[5]*t[32]-.7874512200000001*gradb*t[29]*t[10]*t[63]*t[12]*t[5]*t[68]-0.26248374*grada*t[29]*t[172]*t[12]*t[5]*t[68]+15.7464*t[10]*t[172]*t[5]*t[180]*t[18];
954 d4fdgagbgbab = t[171]-0.90015*grada*t[16]*t[294]*t[12]*t[5]*t[23]-0.324054*grada*gradb*t[166]*t[11]*t[12]*t[5]*t[14]+t[170]-0.162027*t[16]*t[166]*t[10]*t[12]*t[5]*t[14]-0.1944324*gradb*t[60]*t[10]*t[11]*t[12]*t[5]*t[32]-0.0972162*grada*t[60]*t[63]*t[12]*t[5]*t[32]-0.26248374*t[29]*t[10]*t[63]*t[12]*t[5]*t[68];
955 d4fdgagbabab = -0.90015*grada*gradb*t[294]*t[12]*t[5]*t[23]-0.162027*grada*t[166]*t[11]*t[12]*t[5]*t[14]-0.162027*gradb*t[166]*t[10]*t[12]*t[5]*t[14]-0.0972162*t[60]*t[10]*t[11]*t[12]*t[5]*t[32];
956 d4fdgaababab = -0.90015*grada*t[294]*t[12]*t[5]*t[23]-0.162027*t[166]*t[10]*t[12]*t[5]*t[14];
957 d4fdgbgbgbgb = -0.90015*t[12]*t[23]*t[294]*t[5]*pow(gradb,4.0)-1.04993496*gradb*t[29]*t[172]*t[12]*t[5]*t[68]-.5832972000000001*t[16]*t[60]*t[63]*t[12]*t[5]*t[32]+.5832972000000001*t[29]*t[63]*t[12]*t[5]*t[32]+t[295]+1.08018*t[16]*t[166]*t[12]*t[5]*t[23]+15.7464*t[293]*t[5]*t[180]*t[18]+0.648108*gradb*t[60]*t[11]*t[12]*t[5]*t[14]-0.648108*t[173]*t[166]*t[11]*t[12]*t[5]*t[14];
958 d4fdgbgbgbab = 0.54009*gradb*t[166]*t[12]*t[5]*t[23]-0.90015*t[173]*t[294]*t[12]*t[5]*t[23]+0.162027*t[60]*t[11]*t[12]*t[5]*t[14]-0.486081*t[16]*t[166]*t[11]*t[12]*t[5]*t[14]-0.2916486*gradb*t[60]*t[63]*t[12]*t[5]*t[32]-0.26248374*t[29]*t[172]*t[12]*t[5]*t[68];
959 d4fdgbgbabab = t[179]-0.90015*t[16]*t[294]*t[12]*t[5]*t[23]-0.324054*gradb*t[166]*t[11]*t[12]*t[5]*t[14]-0.0972162*t[60]*t[63]*t[12]*t[5]*t[32];
960 d4fdgbababab = -0.90015*gradb*t[294]*t[12]*t[5]*t[23]-0.162027*t[166]*t[11]*t[12]*t[5]*t[14];
961 d4fdabababab = -0.90015*t[294]*t[12]*t[5]*t[23];
962 ds->df1000 += factor*dfdra;
963 ds->df0100 += factor*dfdrb;
964 ds->df0010 += factor*dfdga;
965 ds->df0001 += factor*dfdgb;
966 ds->df00001 += factor*dfdab;
967 ds->df2000 += factor*d2fdrara;
968 ds->df1100 += factor*d2fdrarb;
969 ds->df1010 += factor*d2fdraga;
970 ds->df1001 += factor*d2fdragb;
971 ds->df10001 += factor*d2fdraab;
972 ds->df0200 += factor*d2fdrbrb;
973 ds->df0110 += factor*d2fdrbga;
974 ds->df0101 += factor*d2fdrbgb;
975 ds->df01001 += factor*d2fdrbab;
976 ds->df0020 += factor*d2fdgaga;
977 ds->df0011 += factor*d2fdgagb;
978 ds->df00101 += factor*d2fdgaab;
979 ds->df0002 += factor*d2fdgbgb;
980 ds->df00011 += factor*d2fdgbab;
981 ds->df00002 += factor*d2fdabab;
982 ds->df3000 += factor*d3fdrarara;
983 ds->df2100 += factor*d3fdrararb;
984 ds->df2010 += factor*d3fdraraga;
985 ds->df2001 += factor*d3fdraragb;
986 ds->df20001 += factor*d3fdraraab;
987 ds->df1200 += factor*d3fdrarbrb;
988 ds->df1110 += factor*d3fdrarbga;
989 ds->df1101 += factor*d3fdrarbgb;
990 ds->df11001 += factor*d3fdrarbab;
991 ds->df1020 += factor*d3fdragaga;
992 ds->df1011 += factor*d3fdragagb;
993 ds->df10101 += factor*d3fdragaab;
994 ds->df1002 += factor*d3fdragbgb;
995 ds->df10011 += factor*d3fdragbab;
996 ds->df10002 += factor*d3fdraabab;
997 ds->df0300 += factor*d3fdrbrbrb;
998 ds->df0210 += factor*d3fdrbrbga;
999 ds->df0201 += factor*d3fdrbrbgb;
1000 ds->df02001 += factor*d3fdrbrbab;
1001 ds->df0120 += factor*d3fdrbgaga;
1002 ds->df0111 += factor*d3fdrbgagb;
1003 ds->df01101 += factor*d3fdrbgaab;
1004 ds->df0102 += factor*d3fdrbgbgb;
1005 ds->df01011 += factor*d3fdrbgbab;
1006 ds->df01002 += factor*d3fdrbabab;
1007 ds->df0030 += factor*d3fdgagaga;
1008 ds->df0021 += factor*d3fdgagagb;
1009 ds->df00201 += factor*d3fdgagaab;
1010 ds->df0012 += factor*d3fdgagbgb;
1011 ds->df00111 += factor*d3fdgagbab;
1012 ds->df00102 += factor*d3fdgaabab;
1013 ds->df0003 += factor*d3fdgbgbgb;
1014 ds->df00021 += factor*d3fdgbgbab;
1015 ds->df00012 += factor*d3fdgbabab;
1016 ds->df00003 += factor*d3fdababab;
1017 ds->df4000 += factor*d4fdrararara;
1018 ds->df3100 += factor*d4fdrarararb;
1019 ds->df3010 += factor*d4fdrararaga;
1020 ds->df3001 += factor*d4fdrararagb;
1021 ds->df30001 += factor*d4fdrararaab;
1022 ds->df2200 += factor*d4fdrararbrb;
1023 ds->df2110 += factor*d4fdrararbga;
1024 ds->df2101 += factor*d4fdrararbgb;
1025 ds->df21001 += factor*d4fdrararbab;
1026 ds->df2020 += factor*d4fdraragaga;
1027 ds->df2011 += factor*d4fdraragagb;
1028 ds->df20101 += factor*d4fdraragaab;
1029 ds->df2002 += factor*d4fdraragbgb;
1030 ds->df20011 += factor*d4fdraragbab;
1031 ds->df20002 += factor*d4fdraraabab;
1032 ds->df1300 += factor*d4fdrarbrbrb;
1033 ds->df1210 += factor*d4fdrarbrbga;
1034 ds->df1201 += factor*d4fdrarbrbgb;
1035 ds->df12001 += factor*d4fdrarbrbab;
1036 ds->df1120 += factor*d4fdrarbgaga;
1037 ds->df1111 += factor*d4fdrarbgagb;
1038 ds->df11101 += factor*d4fdrarbgaab;
1039 ds->df1102 += factor*d4fdrarbgbgb;
1040 ds->df11011 += factor*d4fdrarbgbab;
1041 ds->df11002 += factor*d4fdrarbabab;
1042 ds->df1030 += factor*d4fdragagaga;
1043 ds->df1021 += factor*d4fdragagagb;
1044 ds->df10201 += factor*d4fdragagaab;
1045 ds->df1012 += factor*d4fdragagbgb;
1046 ds->df10111 += factor*d4fdragagbab;
1047 ds->df10102 += factor*d4fdragaabab;
1048 ds->df1003 += factor*d4fdragbgbgb;
1049 ds->df10021 += factor*d4fdragbgbab;
1050 ds->df10012 += factor*d4fdragbabab;
1051 ds->df10003 += factor*d4fdraababab;
1052 ds->df0400 += factor*d4fdrbrbrbrb;
1053 ds->df0310 += factor*d4fdrbrbrbga;
1054 ds->df0301 += factor*d4fdrbrbrbgb;
1055 ds->df03001 += factor*d4fdrbrbrbab;
1056 ds->df0220 += factor*d4fdrbrbgaga;
1057 ds->df0211 += factor*d4fdrbrbgagb;
1058 ds->df02101 += factor*d4fdrbrbgaab;
1059 ds->df0202 += factor*d4fdrbrbgbgb;
1060 ds->df02011 += factor*d4fdrbrbgbab;
1061 ds->df02002 += factor*d4fdrbrbabab;
1062 ds->df0130 += factor*d4fdrbgagaga;
1063 ds->df0121 += factor*d4fdrbgagagb;
1064 ds->df01201 += factor*d4fdrbgagaab;
1065 ds->df0112 += factor*d4fdrbgagbgb;
1066 ds->df01111 += factor*d4fdrbgagbab;
1067 ds->df01102 += factor*d4fdrbgaabab;
1068 ds->df0103 += factor*d4fdrbgbgbgb;
1069 ds->df01021 += factor*d4fdrbgbgbab;
1070 ds->df01012 += factor*d4fdrbgbabab;
1071 ds->df01003 += factor*d4fdrbababab;
1072 ds->df0040 += factor*d4fdgagagaga;
1073 ds->df0031 += factor*d4fdgagagagb;
1074 ds->df00301 += factor*d4fdgagagaab;
1075 ds->df0022 += factor*d4fdgagagbgb;
1076 ds->df00211 += factor*d4fdgagagbab;
1077 ds->df00202 += factor*d4fdgagaabab;
1078 ds->df0013 += factor*d4fdgagbgbgb;
1079 ds->df00121 += factor*d4fdgagbgbab;
1080 ds->df00112 += factor*d4fdgagbabab;
1081 ds->df00103 += factor*d4fdgaababab;
1082 ds->df0004 += factor*d4fdgbgbgbgb;
1083 ds->df00031 += factor*d4fdgbgbgbab;
1084 ds->df00022 += factor*d4fdgbgbabab;
1085 ds->df00013 += factor*d4fdgbababab;
1086 ds->df00004 += factor*d4fdabababab;
1087 };
1088 }
1089