-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathSourceTermStrainPath.m
279 lines (277 loc) · 22.7 KB
/
SourceTermStrainPath.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
function [velocity,epsiDot] = SourceTermStrainPath(h,r,z)
%SOURCETERMSTRAINPATH
% [VELOCITY,EPSIDOT] = SOURCETERMSTRAINPATH(H,R,Z)
% This function was generated by the Symbolic Math Toolbox version 8.1.
% 25-Jun-2020 20:26:01
t2 = h-z;
t3 = r-9.0./1.0e1;
t4 = t2.^2;
t5 = r+9.0./1.0e1;
t6 = t3.^2;
t7 = t4+t6;
t8 = t5.^2;
t9 = t4+t8;
t11 = 1.0./t9;
t12 = t7.*t11;
t10 = t12-1.0;
t13 = t10.^2;
t14 = t13.^2;
t15 = t14.^2;
t16 = t15.^2;
t17 = t16.^2;
t18 = t17.^2;
t19 = 1.0./r;
t20 = 1.0./sqrt(t9);
t21 = t7.*t11.*(1.0./4.0);
t22 = 1.0./t7;
t23 = t10.*t13.*(5.0./2.56e2);
t24 = t10.*t14.*6.7291259765625e-3;
t25 = t10.*t13.*t14.*3.37529182434082e-3;
t26 = t10.*t15.*2.023490378633142e-3;
t27 = t10.*t13.*t15.*1.347011206235038e-3;
t28 = t10.*t14.*t15.*9.607646266118763e-4;
t29 = t10.*t13.*t14.*t15.*7.196543711451839e-4;
t30 = t10.*t16.*5.591154616975373e-4;
t31 = t10.*t13.*t16.*4.468698562952858e-4;
t32 = t10.*t14.*t16.*3.653232323596661e-4;
t33 = t10.*t13.*t14.*t16.*3.042212573348881e-4;
t34 = t10.*t15.*t16.*2.572594774623877e-4;
t35 = t10.*t13.*t15.*t16.*2.203887786250346e-4;
t36 = t10.*t14.*t15.*t16.*1.90912137972017e-4;
t37 = t10.*t13.*t14.*t15.*t16.*1.669768928835425e-4;
t38 = t10.*t17.*1.472762938973477e-4;
t39 = t10.*t13.*t17.*1.308672553855569e-4;
t40 = t10.*t14.*t17.*1.170553344682559e-4;
t41 = t10.*t13.*t14.*t17.*1.053200198699657e-4;
t42 = t10.*t15.*t17.*9.526495798354636e-5;
t43 = t10.*t13.*t15.*t17.*8.658403778909685e-5;
t44 = t10.*t14.*t15.*t17.*7.903781203806091e-5;
t45 = t10.*t13.*t14.*t15.*t17.*7.243680514899048e-5;
t46 = t10.*t16.*t17.*6.662950138553806e-5;
t47 = t10.*t13.*t16.*t17.*6.14935714624319e-5;
t48 = t10.*t14.*t16.*t17.*5.692937686305637e-5;
t49 = t10.*t13.*t14.*t16.*t17.*5.285510029195744e-5;
t50 = t10.*t15.*t16.*t17.*4.920305312200104e-5;
t51 = t10.*t13.*t15.*t16.*t17.*4.591684558708133e-5;
t52 = t10.*t14.*t15.*t16.*t17.*4.294919680588934e-5;
t53 = t10.*t13.*t14.*t15.*t16.*t17.*4.026022448623208e-5;
t54 = t10.*t18.*3.781609788488733e-5;
t55 = t10.*t13.*t18.*3.558796845010862e-5;
t56 = t10.*t14.*t18.*3.355111460228601e-5;
t57 = t10.*t13.*t14.*t18.*3.16842530115369e-5;
t58 = t10.*t15.*t18.*2.996898033180698e-5;
t59 = t10.*t13.*t15.*t18.*2.838931789599933e-5;
t63 = t13.*(3.0./6.4e1);
t64 = t14.*1.068115234375e-2;
t65 = t13.*t14.*4.626274108886719e-3;
t66 = t15.*2.571023069322109e-3;
t67 = t13.*t15.*1.633968480746262e-3;
t68 = t14.*t15.*1.129525021895006e-3;
t69 = t13.*t14.*t15.*8.271889323507862e-4;
t70 = t16.*6.318059371675005e-4;
t71 = t13.*t16.*4.982857702628516e-4;
t72 = t14.*t16.*4.030207516463108e-4;
t73 = t13.*t14.*t16.*3.326781294680223e-4;
t74 = t15.*t16.*2.792656073191357e-4;
t75 = t13.*t15.*t16.*2.377557079062525e-4;
t76 = t14.*t15.*t16.*2.048575541109618e-4;
t77 = t13.*t14.*t15.*t16.*1.783437555555258e-4;
t78 = t17.*1.566631346072886e-4;
t79 = t13.*t17.*1.38708533720361e-4;
t80 = t14.*t17.*1.236725856739666e-4;
t81 = t13.*t14.*t17.*1.109553248292418e-4;
t82 = t15.*t17.*1.001033876357815e-4;
t83 = t13.*t15.*t17.*9.076903522156774e-5;
t84 = t14.*t15.*t17.*8.268194207778558e-5;
t85 = t13.*t14.*t15.*t17.*7.562939977507744e-5;
t86 = t16.*t17.*6.944218462362532e-5;
t87 = t13.*t16.*t17.*6.398431018053221e-5;
t88 = t14.*t16.*t17.*5.914549037755907e-5;
t89 = t13.*t14.*t16.*t17.*5.483552375312401e-5;
t90 = t15.*t16.*t17.*5.098005886737827e-5;
t91 = t13.*t15.*t16.*t17.*4.751736588290752e-5;
t92 = t14.*t15.*t16.*t17.*4.439585007700928e-5;
t93 = t13.*t14.*t15.*t16.*t17.*4.157211863046637e-5;
t94 = t18.*3.900946433831388e-5;
t95 = t13.*t18.*3.667666656087644e-5;
t96 = t14.*t18.*3.454703576566286e-5;
t97 = t13.*t14.*t18.*3.259764670261902e-5;
t98 = t15.*t18.*3.080871882096925e-5;
t99 = t13.*t15.*t18.*2.916311248048144e-5;
t60 = t21+t23+t24+t25+t26+t27+t28+t29+t30+t31+t32+t33+t34+t35+t36+t37+t38+t39+t40+t41+t42+t43+t44+t45+t46+t47+t48+t49+t50+t51+t52+t53+t54+t55+t56+t57+t58+t59-t63-t64-t65-t66-t67-t68-t69-t70-t71-t72-t73-t74-t75-t76-t77-t78-t79-t80-t81-t82-t83-t84-t85-t86-t87-t88-t89-t90-t91-t92-t93-t94-t95-t96-t97-t98-t99+3.0./4.0;
t61 = r.*2.0;
t62 = t61-9.0./5.0;
t100 = t11.*t62;
t101 = t61+9.0./5.0;
t102 = 1.0./t9.^2;
t104 = t7.*t101.*t102;
t103 = t100-t104;
t105 = t11.*t62.*(1.0./4.0);
t106 = 1.0./r.^2;
t107 = t13.*(9.0./6.4e1);
t108 = t14.*7.476806640625e-2;
t109 = t13.*t14.*5.088901519775391e-2;
t110 = t15.*3.856534603983164e-2;
t111 = t13.*t15.*3.104540113417897e-2;
t112 = t14.*t15.*2.597907550358514e-2;
t113 = t13.*t14.*t15.*2.233410117347123e-2;
t114 = t16.*1.958598405219252e-2;
t115 = t13.*t16.*1.74400019591998e-2;
t116 = t14.*t16.*1.571780931420612e-2;
t117 = t13.*t14.*t16.*1.430515956712496e-2;
t118 = t15.*t16.*1.312548354399938e-2;
t119 = t13.*t15.*t16.*1.212554110321888e-2;
t120 = t14.*t15.*t16.*1.12671654761029e-2;
t121 = t13.*t14.*t15.*t16.*1.052228157777602e-2;
t122 = t17.*9.869777480259179e-3;
t123 = t13.*t17.*9.293471759264186e-3;
t124 = t14.*t17.*8.780753582851632e-3;
t125 = t13.*t14.*t17.*8.321649362193139e-3;
t126 = t15.*t17.*7.908167623226735e-3;
t127 = t13.*t15.*t17.*7.533829923390122e-3;
t128 = t14.*t15.*t17.*7.193328960767345e-3;
t129 = t13.*t14.*t15.*t17.*6.882275379532047e-3;
t130 = t16.*t17.*6.597007539244406e-3;
t131 = t13.*t16.*t17.*6.334446707872688e-3;
t132 = t14.*t16.*t17.*6.091985508888584e-3;
t133 = t13.*t14.*t16.*t17.*5.867401041584269e-3;
t134 = t15.*t16.*t17.*5.658786534278988e-3;
t135 = t13.*t15.*t16.*t17.*5.464497076534365e-3;
t136 = t14.*t15.*t16.*t17.*5.283106159164104e-3;
t137 = t13.*t14.*t15.*t16.*t17.*5.113370591547364e-3;
t138 = t18.*4.954201970965862e-3;
t139 = t10.*t18.*4.878276627150465e-3;
t140 = t10.*t13.*t18.*4.733199803864447e-3;
t141 = t10.*t14.*t18.*4.596502700513184e-3;
t142 = t10.*t13.*t14.*t18.*4.467479674626702e-3;
t143 = t10.*t15.*t18.*4.345502148112012e-3;
t144 = t10.*t13.*t15.*t18.*4.230008366503901e-3;
t145 = t10.*t13.*(2.5e1./2.56e2);
t146 = t10.*t14.*6.05621337890625e-2;
t147 = t10.*t15.*3.439933643676341e-2;
t148 = t10.*t16.*1.845081023601873e-2;
t149 = t10.*t17.*9.572959103327599e-3;
t150 = t10.*t13.*t14.*4.387879371643066e-2;
t151 = t10.*t13.*t15.*2.82872353309358e-2;
t152 = t10.*t13.*t16.*1.653418468292557e-2;
t153 = t10.*t14.*t15.*2.401911566529691e-2;
t154 = t10.*t13.*t17.*9.029840621603424e-3;
t155 = t10.*t14.*t16.*1.497825252674631e-2;
t156 = t10.*t14.*t17.*8.545039416182678e-3;
t157 = t10.*t15.*t16.*1.2605714395657e-2;
t158 = t10.*t15.*t17.*7.716461596667256e-3;
t159 = t10.*t16.*t17.*6.463061634397192e-3;
t160 = t10.*t13.*t14.*t15.*2.086997676321033e-2;
t161 = t10.*t13.*t14.*t16.*1.368995658006996e-2;
t162 = t10.*t13.*t14.*t17.*8.109641529987362e-3;
t163 = t10.*t13.*t15.*t16.*1.168060526712683e-2;
t164 = t10.*t13.*t15.*t17.*7.359643212073232e-3;
t165 = t10.*t14.*t15.*t16.*1.088199186440497e-2;
t166 = t10.*t13.*t16.*t17.*6.210850717705622e-3;
t167 = t10.*t14.*t15.*t17.*7.03436527138742e-3;
t168 = t10.*t14.*t16.*t17.*5.977584570620919e-3;
t169 = t10.*t15.*t16.*t17.*5.559945002786118e-3;
t170 = t10.*t13.*t14.*t15.*t16.*1.018559046589609e-2;
t171 = t10.*t13.*t14.*t15.*t17.*6.736622878856115e-3;
t172 = t10.*t13.*t14.*t16.*t17.*5.761205931823361e-3;
t173 = t10.*t13.*t15.*t16.*t17.*5.372270933688516e-3;
t174 = t10.*t14.*t15.*t16.*t17.*5.19685281351261e-3;
t175 = t10.*t13.*t14.*t15.*t16.*t17.*5.032528060779011e-3;
t179 = t13.*t18.*4.804643319474814e-3;
t180 = t14.*t18.*4.663849828364487e-3;
t181 = t13.*t14.*t18.*4.531072891664044e-3;
t182 = t15.*t18.*4.405646791398602e-3;
t183 = t13.*t15.*t18.*4.286977534630772e-3;
t176 = t21-t107-t108-t109-t110-t111-t112-t113-t114-t115-t116-t117-t118-t119-t120-t121-t122-t123-t124-t125-t126-t127-t128-t129-t130-t131-t132-t133-t134-t135-t136-t137-t138+t139+t140+t141+t142+t143+t144+t145+t146+t147+t148+t149+t150+t151+t152+t153+t154+t155+t156+t157+t158+t159+t160+t161+t162+t163+t164+t165+t166+t167+t168+t169+t170+t171+t172+t173+t174+t175-t179-t180-t181-t182-t183-5.0./4.0;
t191 = t22.*t62;
t177 = t19-t191;
t178 = 1.0./t9.^(3.0./2.0);
t184 = h.*2.0;
t186 = z.*2.0;
t185 = t184-t186;
t187 = t11.*t185;
t189 = t7.*t102.*t185;
t188 = t187-t189;
t190 = 1.0./t7.^2;
t192 = t19.*t20.*t176.*pi.*(-1.432394487827058e-2)-t20.*t60.*t177.*pi.*1.432394487827058e-2;
velocity = [t192,t2.*t20.*t22.*t60.*pi.*(-2.864788975654117e-2)];
if nargout > 1
t193 = t11.*t185.*(1.0./4.0);
t194 = t13.*t188.*(1.5e1./2.56e2);
t195 = t14.*t188.*3.36456298828125e-2;
t196 = t13.*t14.*t188.*2.362704277038574e-2;
t197 = t15.*t188.*1.821141340769827e-2;
t198 = t13.*t15.*t188.*1.481712326858542e-2;
t199 = t14.*t15.*t188.*1.248994014595439e-2;
t200 = t13.*t14.*t15.*t188.*1.079481556717776e-2;
t201 = t16.*t188.*9.504962848858134e-3;
t202 = t13.*t16.*t188.*8.49052726961043e-3;
t203 = t14.*t16.*t188.*7.671787879552988e-3;
t204 = t13.*t14.*t16.*t188.*6.997088918702426e-3;
t205 = t15.*t16.*t188.*6.431486936559693e-3;
t206 = t13.*t15.*t16.*t188.*5.950497022875934e-3;
t207 = t14.*t15.*t16.*t188.*5.536452001188493e-3;
t208 = t13.*t14.*t15.*t16.*t188.*5.176283679389816e-3;
t209 = t17.*t188.*4.860117698612473e-3;
t210 = t13.*t17.*t188.*4.58035393849449e-3;
t211 = t14.*t17.*t188.*4.331047375325467e-3;
t212 = t13.*t14.*t17.*t188.*4.107480774928664e-3;
t213 = t15.*t17.*t188.*3.905863277325401e-3;
t214 = t13.*t15.*t17.*t188.*3.723113624931165e-3;
t215 = t14.*t15.*t17.*t188.*3.556701541712741e-3;
t216 = t13.*t14.*t15.*t17.*t188.*3.404529842002553e-3;
t217 = t16.*t17.*t188.*3.264845567891365e-3;
t218 = t13.*t16.*t17.*t188.*3.136172144584027e-3;
t219 = t14.*t16.*t17.*t188.*3.017256973741988e-3;
t220 = t13.*t14.*t16.*t17.*t188.*2.907030516057659e-3;
t221 = t15.*t16.*t17.*t188.*2.804574027954059e-3;
t222 = t13.*t15.*t16.*t17.*t188.*2.709093889637799e-3;
t223 = t14.*t15.*t16.*t17.*t188.*2.61990100515925e-3;
t224 = t13.*t14.*t15.*t16.*t17.*t188.*2.536394142632621e-3;
t225 = t18.*t188.*2.458046362517676e-3;
t226 = t13.*t18.*t188.*2.384393886157278e-3;
t227 = t14.*t18.*t188.*2.315026907557735e-3;
t228 = t13.*t14.*t18.*t188.*2.24958196381912e-3;
t229 = t15.*t18.*t188.*2.187735564221909e-3;
t230 = t13.*t15.*t18.*t188.*2.12919884219995e-3;
t231 = t193+t194+t195+t196+t197+t198+t199+t200+t201+t202+t203+t204+t205+t206+t207+t208+t209+t210+t211+t212+t213+t214+t215+t216+t217+t218+t219+t220+t221+t222+t223+t224+t225+t226+t227+t228+t229+t230-t10.*t188.*(3.0./3.2e1)-t10.*t13.*t188.*4.2724609375e-2-t10.*t14.*t188.*2.775764465332031e-2-t10.*t15.*t188.*1.633968480746262e-2-t10.*t16.*t188.*8.969143864731328e-3-t10.*t17.*t188.*4.716090146492273e-3-t10.*t18.*t188.*2.420659993017845e-3-t7.*t102.*t185.*(1.0./4.0)-t10.*t13.*t14.*t188.*2.056818455457687e-2-t10.*t13.*t15.*t188.*1.355430026274007e-2-t10.*t13.*t16.*t188.*8.060415032926216e-3-t10.*t14.*t15.*t188.*1.158064505291101e-2-t10.*t13.*t17.*t188.*4.452213084262799e-3-t10.*t14.*t16.*t188.*7.318918848296491e-3-t10.*t13.*t18.*t188.*2.349198432065075e-3-t10.*t14.*t17.*t188.*4.21630234351119e-3-t10.*t15.*t16.*t188.*6.181648405562565e-3-t10.*t14.*t18.*t188.*2.281835269183331e-3-t10.*t15.*t17.*t188.*3.812299479305845e-3-t10.*t15.*t18.*t188.*2.158070323555627e-3-t10.*t16.*t17.*t188.*3.19921550902661e-3-t10.*t13.*t14.*t15.*t188.*1.010889499468001e-2-t10.*t13.*t14.*t16.*t188.*6.702374575659256e-3-t10.*t13.*t14.*t17.*t188.*4.004135505431258e-3-t10.*t13.*t15.*t16.*t188.*5.736011515106929e-3-t10.*t13.*t14.*t18.*t188.*2.218227755109786e-3-t10.*t13.*t15.*t17.*t188.*3.638005451422565e-3-t10.*t14.*t15.*t16.*t188.*5.350312666665775e-3-t10.*t13.*t16.*t17.*t188.*3.075565499633072e-3-t10.*t14.*t15.*t17.*t188.*3.478952389653562e-3-t10.*t14.*t16.*t17.*t188.*2.961118282668697e-3-t10.*t15.*t16.*t17.*t188.*2.756007221208636e-3-t10.*t13.*t14.*t15.*t16.*t188.*5.013220307433234e-3-t10.*t13.*t14.*t15.*t17.*t188.*3.333224861934016e-3-t10.*t13.*t14.*t16.*t17.*t188.*2.854883296573183e-3-t10.*t13.*t15.*t16.*t17.*t188.*2.663751004620557e-3-t10.*t14.*t15.*t16.*t17.*t188.*2.577471355088915e-3-t10.*t13.*t14.*t15.*t16.*t17.*t188.*2.496605717652088e-3;
t232 = t13.*t103.*(1.5e1./2.56e2);
t233 = t14.*t103.*3.36456298828125e-2;
t234 = t13.*t14.*t103.*2.362704277038574e-2;
t235 = t15.*t103.*1.821141340769827e-2;
t236 = t13.*t15.*t103.*1.481712326858542e-2;
t237 = t14.*t15.*t103.*1.248994014595439e-2;
t238 = t13.*t14.*t15.*t103.*1.079481556717776e-2;
t239 = t16.*t103.*9.504962848858134e-3;
t240 = t13.*t16.*t103.*8.49052726961043e-3;
t241 = t14.*t16.*t103.*7.671787879552988e-3;
t242 = t13.*t14.*t16.*t103.*6.997088918702426e-3;
t243 = t15.*t16.*t103.*6.431486936559693e-3;
t244 = t13.*t15.*t16.*t103.*5.950497022875934e-3;
t245 = t14.*t15.*t16.*t103.*5.536452001188493e-3;
t246 = t13.*t14.*t15.*t16.*t103.*5.176283679389816e-3;
t247 = t17.*t103.*4.860117698612473e-3;
t248 = t13.*t17.*t103.*4.58035393849449e-3;
t249 = t14.*t17.*t103.*4.331047375325467e-3;
t250 = t13.*t14.*t17.*t103.*4.107480774928664e-3;
t251 = t15.*t17.*t103.*3.905863277325401e-3;
t252 = t13.*t15.*t17.*t103.*3.723113624931165e-3;
t253 = t14.*t15.*t17.*t103.*3.556701541712741e-3;
t254 = t13.*t14.*t15.*t17.*t103.*3.404529842002553e-3;
t255 = t16.*t17.*t103.*3.264845567891365e-3;
t256 = t13.*t16.*t17.*t103.*3.136172144584027e-3;
t257 = t14.*t16.*t17.*t103.*3.017256973741988e-3;
t258 = t13.*t14.*t16.*t17.*t103.*2.907030516057659e-3;
t259 = t15.*t16.*t17.*t103.*2.804574027954059e-3;
t260 = t13.*t15.*t16.*t17.*t103.*2.709093889637799e-3;
t261 = t14.*t15.*t16.*t17.*t103.*2.61990100515925e-3;
t262 = t13.*t14.*t15.*t16.*t17.*t103.*2.536394142632621e-3;
t263 = t18.*t103.*2.458046362517676e-3;
t264 = t13.*t18.*t103.*2.384393886157278e-3;
t265 = t14.*t18.*t103.*2.315026907557735e-3;
t266 = t13.*t14.*t18.*t103.*2.24958196381912e-3;
t267 = t15.*t18.*t103.*2.187735564221909e-3;
t268 = t13.*t15.*t18.*t103.*2.12919884219995e-3;
t269 = t105+t232+t233+t234+t235+t236+t237+t238+t239+t240+t241+t242+t243+t244+t245+t246+t247+t248+t249+t250+t251+t252+t253+t254+t255+t256+t257+t258+t259+t260+t261+t262+t263+t264+t265+t266+t267+t268-t10.*t103.*(3.0./3.2e1)-t10.*t13.*t103.*4.2724609375e-2-t10.*t14.*t103.*2.775764465332031e-2-t10.*t15.*t103.*1.633968480746262e-2-t10.*t16.*t103.*8.969143864731328e-3-t10.*t17.*t103.*4.716090146492273e-3-t10.*t18.*t103.*2.420659993017845e-3-t7.*t101.*t102.*(1.0./4.0)-t10.*t13.*t14.*t103.*2.056818455457687e-2-t10.*t13.*t15.*t103.*1.355430026274007e-2-t10.*t13.*t16.*t103.*8.060415032926216e-3-t10.*t14.*t15.*t103.*1.158064505291101e-2-t10.*t13.*t17.*t103.*4.452213084262799e-3-t10.*t14.*t16.*t103.*7.318918848296491e-3-t10.*t13.*t18.*t103.*2.349198432065075e-3-t10.*t14.*t17.*t103.*4.21630234351119e-3-t10.*t15.*t16.*t103.*6.181648405562565e-3-t10.*t14.*t18.*t103.*2.281835269183331e-3-t10.*t15.*t17.*t103.*3.812299479305845e-3-t10.*t15.*t18.*t103.*2.158070323555627e-3-t10.*t16.*t17.*t103.*3.19921550902661e-3-t10.*t13.*t14.*t15.*t103.*1.010889499468001e-2-t10.*t13.*t14.*t16.*t103.*6.702374575659256e-3-t10.*t13.*t14.*t17.*t103.*4.004135505431258e-3-t10.*t13.*t15.*t16.*t103.*5.736011515106929e-3-t10.*t13.*t14.*t18.*t103.*2.218227755109786e-3-t10.*t13.*t15.*t17.*t103.*3.638005451422565e-3-t10.*t14.*t15.*t16.*t103.*5.350312666665775e-3-t10.*t13.*t16.*t17.*t103.*3.075565499633072e-3-t10.*t14.*t15.*t17.*t103.*3.478952389653562e-3-t10.*t14.*t16.*t17.*t103.*2.961118282668697e-3-t10.*t15.*t16.*t17.*t103.*2.756007221208636e-3-t10.*t13.*t14.*t15.*t16.*t103.*5.013220307433234e-3-t10.*t13.*t14.*t15.*t17.*t103.*3.333224861934016e-3-t10.*t13.*t14.*t16.*t17.*t103.*2.854883296573183e-3-t10.*t13.*t15.*t16.*t17.*t103.*2.663751004620557e-3-t10.*t14.*t15.*t16.*t17.*t103.*2.577471355088915e-3-t10.*t13.*t14.*t15.*t16.*t17.*t103.*2.496605717652088e-3;
epsiDot = [t20.*t106.*t176.*pi.*(-1.432394487827058e-2)+t20.*t177.*t269.*pi.*1.432394487827058e-2+t19.*t20.*pi.*(t105-t10.*t103.*(9.0./3.2e1)+t13.*t103.*(7.5e1./2.56e2)+t14.*t103.*3.028106689453125e-1+t15.*t103.*3.095940279308707e-1+t16.*t103.*3.136637740123184e-1+t17.*t103.*3.159076504098108e-1+t18.*t103.*3.170879807647802e-1-t10.*t13.*t103.*2.99072265625e-1-t10.*t14.*t103.*3.053340911865234e-1-t10.*t15.*t103.*3.104540113417897e-1-t10.*t16.*t103.*3.139200352655965e-1-t10.*t17.*t103.*3.159780398149823e-1+t13.*t14.*t103.*3.071515560150146e-1-t10.*t18.*t103.*3.171064590853377e-1+t13.*t15.*t103.*3.111595886402938e-1+t13.*t16.*t103.*3.141495089755859e-1+t14.*t15.*t103.*3.122485036488598e-1+t13.*t17.*t103.*3.160444217561198e-1+t14.*t16.*t103.*3.145433030616725e-1+t13.*t18.*t103.*3.171243868589179e-1+t14.*t17.*t103.*3.161664583987591e-1+t15.*t16.*t103.*3.15142859891425e-1+t14.*t18.*t103.*3.171586863354097e-1+t15.*t17.*t103.*3.163749254633575e-1+t15.*t18.*t103.*3.172216568121768e-1+t16.*t17.*t103.*3.166900200854624e-1-t7.*t101.*t102.*(1.0./4.0)-t10.*t13.*t14.*t103.*3.085227683186531e-1-t10.*t13.*t15.*t103.*3.117489060430216e-1-t10.*t13.*t16.*t103.*3.143561862841224e-1-t10.*t14.*t15.*t103.*3.126774164285972e-1-t10.*t13.*t17.*t103.*3.161071289826588e-1-t10.*t14.*t16.*t103.*3.147135104767491e-1-t10.*t13.*t18.*t103.*3.171417883287851e-1-t10.*t14.*t17.*t103.*3.162226757633393e-1-t10.*t15.*t16.*t103.*3.152640686836908e-1-t10.*t14.*t18.*t103.*3.17175102416483e-1-t10.*t15.*t17.*t103.*3.164208567823851e-1+t13.*t14.*t15.*t103.*3.13049651448155e-1-t10.*t15.*t18.*t103.*3.172363375626771e-1-t10.*t16.*t17.*t103.*3.167223353936344e-1+t13.*t14.*t16.*t103.*3.148690013416092e-1+t13.*t14.*t17.*t103.*3.162760196695071e-1+t13.*t15.*t16.*t103.*3.153763422124245e-1+t13.*t14.*t18.*t103.*3.171910568984959e-1+t13.*t15.*t17.*t103.*3.16464658119149e-1+t14.*t15.*t16.*t103.*3.155777640677441e-1+t13.*t15.*t18.*t103.*3.172506274877925e-1+t13.*t16.*t17.*t103.*3.167533866029867e-1+t14.*t15.*t17.*t103.*3.165464372124339e-1+t14.*t16.*t17.*t103.*3.168119822429087e-1+t15.*t16.*t17.*t103.*3.169168651588087e-1-t10.*t13.*t14.*t15.*t103.*3.133757448350803e-1-t10.*t13.*t14.*t16.*t103.*3.150116050559851e-1-t10.*t13.*t14.*t17.*t103.*3.163267049290694e-1-t10.*t13.*t15.*t16.*t103.*3.154806333308811e-1-t10.*t13.*t14.*t18.*t103.*3.172065689806993e-1-t10.*t13.*t15.*t17.*t103.*3.165064742737632e-1-t10.*t14.*t15.*t16.*t103.*3.156684473332807e-1-t10.*t13.*t16.*t17.*t103.*3.167832464622064e-1-t10.*t14.*t15.*t17.*t103.*3.165846674584741e-1-t10.*t14.*t16.*t17.*t103.*3.168396562455505e-1-t10.*t15.*t16.*t17.*t103.*3.169408304389932e-1+t13.*t14.*t15.*t16.*t103.*3.157533044427788e-1+t13.*t14.*t15.*t17.*t103.*3.166212753062374e-1+t13.*t14.*t16.*t17.*t103.*3.168663262502848e-1+t13.*t15.*t16.*t17.*t103.*3.169639850876224e-1+t14.*t15.*t16.*t17.*t103.*3.170080216242692e-1-t10.*t13.*t14.*t15.*t16.*t103.*3.158328793682937e-1-t10.*t13.*t14.*t15.*t17.*t103.*3.166563618837315e-1-t10.*t13.*t14.*t16.*t17.*t103.*3.168920459196233e-1-t10.*t13.*t15.*t16.*t17.*t103.*3.169863695498462e-1-t10.*t14.*t15.*t16.*t17.*t103.*3.170289766759366e-1+t13.*t14.*t15.*t16.*t17.*t103.*3.170492678290777e-1-t10.*t13.*t14.*t15.*t16.*t17.*t103.*3.170689261418152e-1).*1.432394487827058e-2-t20.*t60.*pi.*(t22.*2.0+t106-t62.^2.*t190).*1.432394487827058e-2-t19.*t101.*t176.*t178.*pi.*7.161972439135291e-3-t60.*t101.*t177.*t178.*pi.*7.161972439135291e-3,t20.*t22.*t60.*pi.*(-2.864788975654117e-2)-t2.*t20.*t22.*t231.*pi.*2.864788975654117e-2+t2.*t22.*t60.*t178.*t185.*pi.*1.432394487827058e-2+t2.*t20.*t60.*t185.*t190.*pi.*2.864788975654117e-2,-t19.*t192,t20.*t177.*t231.*pi.*(-7.161972439135291e-3)-t19.*t20.*pi.*(t193-t10.*t188.*(9.0./3.2e1)+t13.*t188.*(7.5e1./2.56e2)+t14.*t188.*3.028106689453125e-1+t15.*t188.*3.095940279308707e-1+t16.*t188.*3.136637740123184e-1+t17.*t188.*3.159076504098108e-1+t18.*t188.*3.170879807647802e-1-t10.*t13.*t188.*2.99072265625e-1-t10.*t14.*t188.*3.053340911865234e-1-t10.*t15.*t188.*3.104540113417897e-1-t10.*t16.*t188.*3.139200352655965e-1-t10.*t17.*t188.*3.159780398149823e-1+t13.*t14.*t188.*3.071515560150146e-1-t10.*t18.*t188.*3.171064590853377e-1+t13.*t15.*t188.*3.111595886402938e-1+t13.*t16.*t188.*3.141495089755859e-1+t14.*t15.*t188.*3.122485036488598e-1+t13.*t17.*t188.*3.160444217561198e-1+t14.*t16.*t188.*3.145433030616725e-1+t13.*t18.*t188.*3.171243868589179e-1+t14.*t17.*t188.*3.161664583987591e-1+t15.*t16.*t188.*3.15142859891425e-1+t14.*t18.*t188.*3.171586863354097e-1+t15.*t17.*t188.*3.163749254633575e-1+t15.*t18.*t188.*3.172216568121768e-1+t16.*t17.*t188.*3.166900200854624e-1-t7.*t102.*t185.*(1.0./4.0)-t10.*t13.*t14.*t188.*3.085227683186531e-1-t10.*t13.*t15.*t188.*3.117489060430216e-1-t10.*t13.*t16.*t188.*3.143561862841224e-1-t10.*t14.*t15.*t188.*3.126774164285972e-1-t10.*t13.*t17.*t188.*3.161071289826588e-1-t10.*t14.*t16.*t188.*3.147135104767491e-1-t10.*t13.*t18.*t188.*3.171417883287851e-1-t10.*t14.*t17.*t188.*3.162226757633393e-1-t10.*t15.*t16.*t188.*3.152640686836908e-1-t10.*t14.*t18.*t188.*3.17175102416483e-1-t10.*t15.*t17.*t188.*3.164208567823851e-1+t13.*t14.*t15.*t188.*3.13049651448155e-1-t10.*t15.*t18.*t188.*3.172363375626771e-1-t10.*t16.*t17.*t188.*3.167223353936344e-1+t13.*t14.*t16.*t188.*3.148690013416092e-1+t13.*t14.*t17.*t188.*3.162760196695071e-1+t13.*t15.*t16.*t188.*3.153763422124245e-1+t13.*t14.*t18.*t188.*3.171910568984959e-1+t13.*t15.*t17.*t188.*3.16464658119149e-1+t14.*t15.*t16.*t188.*3.155777640677441e-1+t13.*t15.*t18.*t188.*3.172506274877925e-1+t13.*t16.*t17.*t188.*3.167533866029867e-1+t14.*t15.*t17.*t188.*3.165464372124339e-1+t14.*t16.*t17.*t188.*3.168119822429087e-1+t15.*t16.*t17.*t188.*3.169168651588087e-1-t10.*t13.*t14.*t15.*t188.*3.133757448350803e-1-t10.*t13.*t14.*t16.*t188.*3.150116050559851e-1-t10.*t13.*t14.*t17.*t188.*3.163267049290694e-1-t10.*t13.*t15.*t16.*t188.*3.154806333308811e-1-t10.*t13.*t14.*t18.*t188.*3.172065689806993e-1-t10.*t13.*t15.*t17.*t188.*3.165064742737632e-1-t10.*t14.*t15.*t16.*t188.*3.156684473332807e-1-t10.*t13.*t16.*t17.*t188.*3.167832464622064e-1-t10.*t14.*t15.*t17.*t188.*3.165846674584741e-1-t10.*t14.*t16.*t17.*t188.*3.168396562455505e-1-t10.*t15.*t16.*t17.*t188.*3.169408304389932e-1+t13.*t14.*t15.*t16.*t188.*3.157533044427788e-1+t13.*t14.*t15.*t17.*t188.*3.166212753062374e-1+t13.*t14.*t16.*t17.*t188.*3.168663262502848e-1+t13.*t15.*t16.*t17.*t188.*3.169639850876224e-1+t14.*t15.*t16.*t17.*t188.*3.170080216242692e-1-t10.*t13.*t14.*t15.*t16.*t188.*3.158328793682937e-1-t10.*t13.*t14.*t15.*t17.*t188.*3.166563618837315e-1-t10.*t13.*t14.*t16.*t17.*t188.*3.168920459196233e-1-t10.*t13.*t15.*t16.*t17.*t188.*3.169863695498462e-1-t10.*t14.*t15.*t16.*t17.*t188.*3.170289766759366e-1+t13.*t14.*t15.*t16.*t17.*t188.*3.170492678290777e-1-t10.*t13.*t14.*t15.*t16.*t17.*t188.*3.170689261418152e-1).*7.161972439135291e-3+t19.*t178.*t185.*pi.*(t21-t107-t108-t109-t110-t111-t112-t113-t114-t115-t116-t117-t118-t119-t120-t121-t122-t123-t124-t125-t126-t127-t128-t129-t130-t131-t132-t133-t134-t135-t136-t137-t138+t139+t140+t141+t142+t143+t144+t145+t146+t147+t148+t149+t150+t151+t152+t153+t154+t155+t156+t157+t158+t159+t160+t161+t162+t163+t164+t165+t166+t167+t168+t169+t170+t171+t172+t173+t174+t175-t179-t180-t181-t182-t183-5.0./4.0).*3.580986219567646e-3+t2.*t20.*t22.*t269.*pi.*1.432394487827058e-2+t60.*t178.*t185.*pi.*(t19-t191).*3.580986219567646e-3-t2.*t20.*t60.*t62.*t190.*pi.*1.432394487827058e-2-t2.*t22.*t60.*t101.*t178.*pi.*7.161972439135291e-3-t20.*t60.*t62.*t185.*t190.*pi.*7.161972439135291e-3];
end