-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathElasto_plastic_diaphragm_bulging.m
250 lines (187 loc) · 103 KB
/
Elasto_plastic_diaphragm_bulging.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
clear;
clc;
close all;
tic
global v E ho ro dr
% With the below parameters the code will take approx. 5 minutes
%% Geometry
ro = 0.1/2; % Outer radius: meters
dr = 0.0005*ro; % Mesh/grid size in radial domain: meters
r = 0.0000001:dr:ro;
%% Material properties
ho = 200e-6; % Initial thickness of the membrane: meters
Sy = 14.16e6; % Yield strength: Pa
v = 0.4; % Poisson's ratio
E = 1.1e9; % Elastic Modulus: Pa
Et = 96.8e6; % Tangent Modulus: Pa
H=E*Et/(E-Et); % Plastic Modulus: Pa
%% Elastic deformation
so=1e6:1e6:Sy;
for i=1:length(so)
disp(['Equivalent stress at the apex for next step: ' num2str(so(i)/1e6) ' MPa'] )
[w_max_e(i),er_center_e(i),er_edge_e(i),phi(:,i),s(:,i),er(:,i),eth(:,i),eh(:,i),sr(:,i),sth(:,i),w(:,i),dw(:,i),dphi(:,i),pj(i)] = Elastic_range (so(i));
end
jstart = length(pj);
p=pj(jstart);
deth(:,jstart) = gradient(eth(:,jstart),r);
deh(:,jstart) = gradient(eh(:,jstart),r);
ds(:,jstart) = gradient(s(:,jstart),r);
dphi(:,jstart) = gradient(phi(1:end-1,jstart),r(1:end-1));
del_s = 0.05e6; % Stress increment at the apex: Pa
s_end = 18e6; % Maximum stress at the apex: Pa
counter = jstart + round( (s_end-Sy)/del_s );
for j=jstart:counter
s(1,j+1) = s(1,j) + del_s;
disp(['Equivalent stress at the apex for next step: ' num2str(s(1,j+1)/1e6) ' MPa'] )
phi(1,j+1) = atan(1/sqrt(3));
eth(1,j+1) = eth(1,j)+((-1).*s(1,j)+s(1,j+1)).*(E.^(-1).*( ...
3.^(-1/2)+(-2).*3.^(-1/2).*v).*cos(phi(1,j))+E.^(-1).*sin(phi(1,j)))...
+((-1).*phi(1,j)+phi(1,j+1)).*(E.^(-1).*cos(phi(1,j)).*s(1,j)+E.^(-1).*((-1).*3.^(-1/2).*s(1,j)+2.*3.^(-1/2).*v.*s(1,j)).* ...
sin(phi(1,j)));
eth(end-1,j+1)=1;
p_inc = 10^(round(log10(p))-3); % Pressure increment: Pa
iter = 1;
while abs(eth(1,j+1)/eth(end-1,j+1)) < 100 || eth(end-1,j+1) < 0
for i=1:1:length(r)-1
sr(i,j+1) = (2.*3.^(-1/2).*cos(phi(i,j+1)).*s(i,j+1));
sth(i,j+1) = s(i,j+1).*(1/sqrt(3)*cos(phi(i,j+1)) + sin(phi(i,j+1)));
if s(i,j) > Sy
rp(i,j)=r(i);
er(i,j+1) = er(i,j) + (1/2).*((-1).*s(i,j)+s(i,j+1)).*((2.*3.^( ...
-1/2).*E.^(-1)+(1/2).*3.^(1/2).*H.^(-1)+(-1).*3.^(-1/2).*E.^(-1) ...
.*v).*cos(phi(i,j))+(2.*3.^(-1/2).*E.^(-1)+(1/2).*3.^(1/2).*H.^( ...
-1)+(-1).*3.^(-1/2).*E.^(-1).*v).*cos(phi(i,j+1))+((-1/2).*H.^( ...
-1)+(-1).*E.^(-1).*v).*sin(phi(i,j))+((-1/2).*H.^(-1)+(-1).*E.^( ...
-1).*v).*sin(phi(i,j+1)))+(1/2).*((-1).*phi(i,j)+phi(i,j+1)).*(( ...
-1).*E.^(-1).*v.*cos(phi(i,j)).*s(i,j)+(-1).*E.^(-1).*v.*cos( ...
phi(i,j+1)).*s(i,j+1)+((-2).*3.^(-1/2).*E.^(-1).*s(i,j)+3.^( ...
-1/2).*E.^(-1).*v.*s(i,j)).*sin(phi(i,j))+((-2).*3.^(-1/2).*E.^( ...
-1).*s(i,j+1)+3.^(-1/2).*E.^(-1).*v.*s(i,j+1)).*sin(phi(i,j+1) ...
));
eth(i,j+1) = eth(i,j)+(1/2).*((-1).*s(i,j)+s(i,j+1)).*((3.^( ...
-1/2).*E.^(-1)+(-2).*3.^(-1/2).*E.^(-1).*v).*cos(phi(i,j))+(3.^( ...
-1/2).*E.^(-1)+(-2).*3.^(-1/2).*E.^(-1).*v).*cos(phi(i,j+1))+( ...
E.^(-1)+H.^(-1)).*sin(phi(i,j))+(E.^(-1)+H.^(-1)).*sin(phi(i,j+1)))+(1/2).*((-1).*phi(i,j)+phi(i,j+1)).*(E.^(-1).*cos(phi(i,j)) ...
.*s(i,j)+E.^(-1).*cos(phi(i,j+1)).*s(i,j+1)+((-1).*3.^(-1/2).* ...
E.^(-1).*s(i,j)+2.*3.^(-1/2).*E.^(-1).*v.*s(i,j)).*sin(phi(i,j)) ...
+((-1).*3.^(-1/2).*E.^(-1).*s(i,j+1)+2.*3.^(-1/2).*E.^(-1).*v.* ...
s(i,j+1)).*sin(phi(i,j+1)));
eh(i,j+1) = (eh(i,j)+(1/2).*((-1).*s(i,j)+s(i,j+1)).*(((-1/2).* ...
3.^(1/2).*H.^(-1)+(-1).*3.^(1/2).*E.^(-1).*v).*cos(phi(i,j))+(( ...
-1/2).*3.^(1/2).*H.^(-1)+(-1).*3.^(1/2).*E.^(-1).*v).*cos(phi(i,j+1))+((-1/2).*H.^(-1)+(-1).*E.^(-1).*v).*sin(phi(i,j))+((-1/2) ...
.*H.^(-1)+(-1).*E.^(-1).*v).*sin(phi(i,j+1)))+(1/2).*((-1).*phi(i,j)+phi(i,j+1)).*((-1).*E.^(-1).*v.*cos(phi(i,j)).*s(i,j)+(-1) ...
.*E.^(-1).*v.*cos(phi(i,j+1)).*s(i,j+1)+3.^(1/2).*E.^(-1).*v.* ...
s(i,j).*sin(phi(i,j))+3.^(1/2).*E.^(-1).*v.*s(i,j+1).*sin(phi(i,j+1))));
e(i,j+1) = sqrt(2)/3*sqrt( (er(i,j+1)-eth(i,j+1)).^2 + (er(i,j+1)-eh(i,j+1)).^2 + (eh(i,j+1)-eth(i,j+1)).^2 );
dw(i,j+1) = (-1/2).*exp(-eh(i,j+1)+er(i,j+1)+eth(i,j+1)).*ho.^(-1).*p.*r(i).*sr(i,j+1).^(-1);
ds(i,j+1) = (exp((-1).*eth(i,j+1)).*r(i).^(-1).*(s(i,j).*(8.*E.*H.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*H.*cos(phi(i,j))+3.^(1/2).*(E+H).*cos(phi(i,j+1))+(-1).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1))))+(-1).*(E.*(3.*3.^(1/2).*E+3.^(1/2).*H+4.*3.^(1/2).*H.*v+2.*3.^(1/2).*H.*((-1)+2.*v).*cos(2.*phi(i,j))+3.^(1/2).*(3.*E+H+4.*H.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*H.*v.*cos(phi(i,j)+phi(i,j+1))+3.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+(-1).*E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))+H.*phi(i,j).*(4.*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))).*sr(i,j+1))+s(i,j+1).*((-8).*3.^(1/2).*E.^2.*H+(-12).*3.^(1/2).*E.*H.^2+(-4).*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(2.*phi(i,j+1))+4.*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+phi(i,j+1))+4.*E.*H.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*E.*H.^2.*sin(2.*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(2.*phi(i,j+1))+(-4).*E.*H.^2.*sin(phi(i,j)+phi(i,j+1))+8.*E.*H.^2.*v.*sin(phi(i,j)+phi(i,j+1))+3.*3.^(1/2).*E.^2.*sr(i,j+1)+3.*3.^(1/2).*E.*H.*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*sr(i,j+1)+3.*3.^(1/2).*E.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+3.*3.^(1/2).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(2.*phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(2.*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j).^2.*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j+1).^2.*sr(i,j+1)+3.*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+(-6).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+H.*phi(i,j).*(8.*E.*H.*((-1)+2.*v)+(-1).*((-6).*E+12.*E.*v+3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.*v.*(1+v).*phi(i,j+1)+3.^(1/2).*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*E.*sin(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1))+(-1).*H.*phi(i,j+1).*(8.*E.*H.*((-1)+2.*v)+(-1).*(3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*sin(phi(i,j)+(-1).*phi(i,j+1))+2.*E.*((-1)+2.*v).*(3+3.^(1/2).*sin(phi(i,j)+phi(i,j+1)))).*sr(i,j+1)))).^(-1).*(exp(eth(i,j+1)).*r(i).*s(i,j+1).^2.*(8.*E.*H.*(3.^(1/2).*(E+H).*cos(phi(i,j))+(-1).*3.^(1/2).*H.*cos(phi(i,j+1))+H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sin(phi(i,j+1))+((-3).*E.*H+6.*E.*H.*v+(-3).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+6.*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+12.*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+H.*phi(i,j).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+(-1).*H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+2.*3.^(1/2).*E.*H.*sin(2.*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(2.*phi(i,j+1))+2.*3.^(1/2).*E.*H.*sin(phi(i,j)+phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1)).*dphi(i,j)+(-1).*s(i,j).*((-4).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sth(i,j+1)+sr(i,j+1).*((-12).*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j))+6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+(-12).*E.*H.^2.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j))+4.*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-8).*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+8.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*deh(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*((-2).*H.*v.*cos(phi(i,j))+(E+2.*H.*v).*cos(phi(i,j+1))+3.^(1/2).*(2.*H.*v.*sin(phi(i,j))+(-1).*(E+2.*H.*v).*sin(phi(i,j+1)))).*deth(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))+s(i,j+1).*(4.*E.*H.*((-1).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*(3.*E.*cos(phi(i,j+1))+6.*H.*cos(phi(i,j+1))+H.*phi(i,j+1).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1)))+(-2).*3.^(1/2).*H.*sin(phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j+1))+H.*phi(i,j).*((-1).*3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+3.*sin(phi(i,j+1)))).*sth(i,j+1)+(-2).*sin(phi(i,j+1)).*((-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H+2.*exp(er(i,j+1)).*E.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*(3.^(1/2).*E.*cos(phi(i,j))+2.*3.^(1/2).*H.*cos(phi(i,j))+(-2).*H.*sin(phi(i,j))+4.*H.*v.*sin(phi(i,j))+(-1).*H.*phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+H.*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j+1)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*H.*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j+1)).*ds(i,j)))+sr(i,j+1).*(6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+24.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+24.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-24).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-8).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+16.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+24.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deh(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*cos(phi(i,j+1)).*deh(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*sin(phi(i,j+1)).*deh(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deh(i,j)+6.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deth(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*cos(phi(i,j+1)).*deth(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*r(i).*sin(phi(i,j+1)).*deth(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j+1).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+H.*phi(i,j+1).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j))+(-1).*H.*phi(i,j).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j))))));
dphi(i,j+1) = (-1).*r(i).^(-1).*(s(i,j).*(3.*H.*cos(phi(i,j))+(-3).*(E+H).*cos(phi(i,j+1))+3.^(1/2).*H.*((-1)+2.*v).*...
(sin(phi(i,j))+(-1).*sin(phi(i,j+1))))+s(i,j+1).*(3.*E.*cos(phi(i,j+1))+6.*H.*cos(phi(i,j+1))+H.*phi(i,j+1).*(3.^...
(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1)))+(-2).*3.^(1/2).*H.*sin(phi(i,j+1))+4.*3.^(1/2).*H.*v.*...
sin(phi(i,j+1))+H.*phi(i,j).*((-1).*3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+3.*sin(phi(i,j+1))))).^(-1).*(6.*E.*H+...
(-6).*exp(er(i,j+1)+(-1).*eth(i,j+1)).*E.*H.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+6.*E.*H.*r(i).*...
deth(i,j)+3.^(1/2).*H.*r(i).*(1+(-2).*v).*cos(phi(i,j)).*(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+3.*H.*...
r(i).*(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*sin(phi(i,j)).*dphi(i,j)+(-1).*r(i).*(s(i,j)+(-1).*s(i,j+1)).*(3.*(E+H)...
.*cos(phi(i,j))+3.^(1/2).*H.*((-1)+2.*v).*sin(phi(i,j))).*dphi(i,j)+(-1).*H.*r(i).*(s(i,j).*(3.*cos(phi(i,j))+3.^...
(1/2).*((-1)+2.*v).*sin(phi(i,j)))+s(i,j+1).*(3.*cos(phi(i,j+1))+3.^(1/2).*((-1)+2.*v).*sin(phi(i,j+1)))).*dphi(i,j)...
+(-3).*H.*r(i).*cos(phi(i,j)).*(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*H.*r(i).*(1+(-2).*v).*(phi(i,j)+(-1).*...
phi(i,j+1)).*sin(phi(i,j)).*ds(i,j)+3.*exp((-1).*eth(i,j+1)).*H.*cos(phi(i,j+1)).*((-1).*phi(i,j)+phi(i,j+1)).*(s(i,j)...
.*(8.*E.*H.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*H.*cos(phi(i,j))+3.^(1/2).*(E+H).*cos(phi(i,j+1))+(-1).*H.*((-1)+2.*v)...
.*(sin(phi(i,j))+(-1).*sin(phi(i,j+1))))+(-1).*(E.*(3.*3.^(1/2).*E+3.^(1/2).*H+4.*3.^(1/2).*H.*v+2.*3.^(1/2).*H.*...
((-1)+2.*v).*cos(2.*phi(i,j))+3.^(1/2).*(3.*E+H+4.*H.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*H.*cos(phi(i,j)...
+phi(i,j+1))+4.*3.^(1/2).*H.*v.*cos(phi(i,j)+phi(i,j+1))+3.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*H.*v.*sin(phi(i,j)+...
(-1).*phi(i,j+1)))+H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+(-1).*E.*((-1)+2.*v).*...
((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))+H.*phi(i,j).*(4.*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+E.*((-1)...
+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))).*sr(i,j+1))+s(i,j+1).*((-8).*3.^(1/2).*E.^2.*H+(-12).*3.^(1/2).*E.*H.^2 ...
+(-4).*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^...
(1/2).*E.*H.^2.*cos(2.*phi(i,j+1))+4.*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+phi(i,j+1))...
+4.*E.*H.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*E.*H.^2.*sin(2.*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(2.*phi(i,j+1))+(-4).*E.*H.^2.*sin(phi(i,j)+phi(i,j+1))+8.*E.*H.^2.*v.*sin(phi(i,j)+phi(i,j+1))+3.*3.^(1/2).*E.^2.*sr(i,j+1)+3.*3.^(1/2).*E.*H.*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*sr(i,j+1)+3.*3.^(1/2).*E.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+3.*3.^(1/2).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(2.*phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(2.*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j).^2.*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j+1).^2.*sr(i,j+1)+3.*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+(-6).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+H.*phi(i,j).*(8.*E.*H.*((-1)+2.*v)+(-1).*((-6).*E+12.*E.*v+3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.*v.*(1+v).*phi(i,j+1)+3.^(1/2).*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*E.*sin(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1))+(-1).*H.*phi(i,j+1).*(8.*E.*H.*((-1)+2.*v)+(-1).*(3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*sin(phi(i,j)+(-1).*phi(i,j+1))+2.*E.*((-1)+2.*v).*(3+3.^(1/2).*sin(phi(i,j)+phi(i,j+1)))).*sr(i,j+1)))).^(-1).*(exp(eth(i,j+1)).*r(i).*s(i,j+1).^2.*(8.*E.*H.*(3.^(1/2).*(E+H).*cos(phi(i,j))+(-1).*3.^(1/2).*H.*cos(phi(i,j+1))+H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sin(phi(i,j+1))+((-3).*E.*H+6.*E.*H.*v+(-3).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+6.*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+12.*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+H.*phi(i,j).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+(-1).*H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+2.*3.^(1/2).*E.*H.*sin(2.*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(2.*phi(i,j+1))+2.*3.^(1/2).*E.*H.*sin(phi(i,j)+phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1)).*dphi(i,j)+(-1).*s(i,j).*((-4).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sth(i,j+1)+sr(i,j+1).*((-12).*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j))+6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+(-12).*E.*H.^2.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j))+4.*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-8).*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+8.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*deh(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*((-2).*H.*v.*cos(phi(i,j))+(E+2.*H.*v).*cos(phi(i,j+1))+3.^(1/2).*(2.*H.*v.*sin(phi(i,j))+(-1).*(E+2.*H.*v).*sin(phi(i,j+1)))).*deth(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))+s(i,j+1).*(4.*E.*H.*((-1).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*(3.*E.*cos(phi(i,j+1))+6.*H.*cos(phi(i,j+1))+H.*phi(i,j+1).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1)))+(-2).*3.^(1/2).*H.*sin(phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j+1))+H.*phi(i,j).*((-1).*3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+3.*sin(phi(i,j+1)))).*sth(i,j+1)+(-2).*sin(phi(i,j+1)).*((-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H+2.*exp(er(i,j+1)).*E.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*(3.^(1/2).*E.*cos(phi(i,j))+2.*3.^(1/2).*H.*cos(phi(i,j))+(-2).*H.*sin(phi(i,j))+4.*H.*v.*sin(phi(i,j))+(-1).*H.*phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+H.*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j+1)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*H.*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j+1)).*ds(i,j)))+sr(i,j+1).*(6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+24.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+24.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-24).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-8).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+16.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+24.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deh(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*cos(phi(i,j+1)).*deh(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*sin(phi(i,j+1)).*deh(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deh(i,j)+6.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deth(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*cos(phi(i,j+1)).*deth(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*r(i).*sin(phi(i,j+1)).*deth(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j+1).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+H.*phi(i,j+1).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j))+(-1).*H.*phi(i,j).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))))+3.^(1/2).*exp((-1).*eth(i,j+1)).*H.*(1+(-2).*v).*(phi(i,j)+(-1).*phi(i,j+1)).*sin(phi(i,j+1)).*(s(i,j).*(8.*E.*H.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*H.*cos(phi(i,j))+3.^(1/2).*(E+H).*cos(phi(i,j+1))+(-1).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1))))+(-1).*(E.*(3.*3.^(1/2).*E+3.^(1/2).*H+4.*3.^(1/2).*H.*v+2.*3.^(1/2).*H.*((-1)+2.*v).*cos(2.*phi(i,j))+3.^(1/2).*(3.*E+H+4.*H.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*H.*v.*cos(phi(i,j)+phi(i,j+1))+3.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+(-1).*E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))+H.*phi(i,j).*(4.*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))).*sr(i,j+1))+s(i,j+1).*((-8).*3.^(1/2).*E.^2.*H+(-12).*3.^(1/2).*E.*H.^2+(-4).*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(2.*phi(i,j+1))+4.*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+phi(i,j+1))+4.*E.*H.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*E.*H.^2.*sin(2.*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(2.*phi(i,j+1))+(-4).*E.*H.^2.*sin(phi(i,j)+phi(i,j+1))+8.*E.*H.^2.*v.*sin(phi(i,j)+phi(i,j+1))+3.*3.^(1/2).*E.^2.*sr(i,j+1)+3.*3.^(1/2).*E.*H.*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*sr(i,j+1)+3.*3.^(1/2).*E.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+3.*3.^(1/2).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(2.*phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(2.*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j).^2.*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j+1).^2.*sr(i,j+1)+3.*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+(-6).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+H.*phi(i,j).*(8.*E.*H.*((-1)+2.*v)+(-1).*((-6).*E+12.*E.*v+3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.*v.*(1+v).*phi(i,j+1)+3.^(1/2).*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*E.*sin(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1))+(-1).*H.*phi(i,j+1).*(8.*E.*H.*((-1)+2.*v)+(-1).*(3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*sin(phi(i,j)+(-1).*phi(i,j+1))+2.*E.*((-1)+2.*v).*(3+3.^(1/2).*sin(phi(i,j)+phi(i,j+1)))).*sr(i,j+1)))).^(-1).*(exp(eth(i,j+1)).*r(i).*s(i,j+1).^2.*(8.*E.*H.*(3.^(1/2).*(E+H).*cos(phi(i,j))+(-1).*3.^(1/2).*H.*cos(phi(i,j+1))+H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sin(phi(i,j+1))+((-3).*E.*H+6.*E.*H.*v+(-3).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+6.*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+12.*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+H.*phi(i,j).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+(-1).*H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+2.*3.^(1/2).*E.*H.*sin(2.*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(2.*phi(i,j+1))+2.*3.^(1/2).*E.*H.*sin(phi(i,j)+phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1)).*dphi(i,j)+(-1).*s(i,j).*((-4).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sth(i,j+1)+sr(i,j+1).*((-12).*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j))+6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+(-12).*E.*H.^2.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j))+4.*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-8).*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+8.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*deh(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*((-2).*H.*v.*cos(phi(i,j))+(E+2.*H.*v).*cos(phi(i,j+1))+3.^(1/2).*(2.*H.*v.*sin(phi(i,j))+(-1).*(E+2.*H.*v).*sin(phi(i,j+1)))).*deth(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))+s(i,j+1).*(4.*E.*H.*((-1).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*(3.*E.*cos(phi(i,j+1))+6.*H.*cos(phi(i,j+1))+H.*phi(i,j+1).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1)))+(-2).*3.^(1/2).*H.*sin(phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j+1))+H.*phi(i,j).*((-1).*3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+3.*sin(phi(i,j+1)))).*sth(i,j+1)+(-2).*sin(phi(i,j+1)).*((-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H+2.*exp(er(i,j+1)).*E.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*(3.^(1/2).*E.*cos(phi(i,j))+2.*3.^(1/2).*H.*cos(phi(i,j))+(-2).*H.*sin(phi(i,j))+4.*H.*v.*sin(phi(i,j))+(-1).*H.*phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+H.*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j+1)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*H.*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j+1)).*ds(i,j)))+sr(i,j+1).*(6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+24.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+24.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-24).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-8).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+16.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+24.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deh(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*cos(phi(i,j+1)).*deh(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*sin(phi(i,j+1)).*deh(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deh(i,j)+6.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deth(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*cos(phi(i,j+1)).*deth(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*r(i).*sin(phi(i,j+1)).*deth(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j+1).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+H.*phi(i,j+1).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j))+(-1).*H.*phi(i,j).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))))+r(i).*(3.^(1/2).*H.*((-1)+2.*v).*cos(phi(i,j))+3.^(1/2).*H.*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*(E+H).*(sin(phi(i,j))+sin(phi(i,j+1)))).*(ds(i,j)+(-1).*exp((-1).*eth(i,j+1)).*r(i).^(-1).*(s(i,j).*(8.*E.*H.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*H.*cos(phi(i,j))+3.^(1/2).*(E+H).*cos(phi(i,j+1))+(-1).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1))))+(-1).*(E.*(3.*3.^(1/2).*E+3.^(1/2).*H+4.*3.^(1/2).*H.*v+2.*3.^(1/2).*H.*((-1)+2.*v).*cos(2.*phi(i,j))+3.^(1/2).*(3.*E+H+4.*H.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*H.*v.*cos(phi(i,j)+phi(i,j+1))+3.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+(-1).*E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))+H.*phi(i,j).*(4.*3.^(1/2).*H.*v.*(1+v).*sin(phi(i,j)+(-1).*phi(i,j+1))+E.*((-1)+2.*v).*((-3)+2.*3.^(1/2).*sin(2.*phi(i,j+1))))).*sr(i,j+1))+s(i,j+1).*((-8).*3.^(1/2).*E.^2.*H+(-12).*3.^(1/2).*E.*H.^2+(-4).*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+(-1).*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*cos(2.*phi(i,j+1))+4.*3.^(1/2).*E.^2.*H.*cos(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*H.^2.*cos(phi(i,j)+phi(i,j+1))+4.*E.*H.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*E.*H.^2.*sin(2.*phi(i,j+1))+(-8).*E.*H.^2.*v.*sin(2.*phi(i,j+1))+(-4).*E.*H.^2.*sin(phi(i,j)+phi(i,j+1))+8.*E.*H.^2.*v.*sin(phi(i,j)+phi(i,j+1))+3.*3.^(1/2).*E.^2.*sr(i,j+1)+3.*3.^(1/2).*E.*H.*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*sr(i,j+1)+3.*3.^(1/2).*E.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+3.*3.^(1/2).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+12.*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+8.*3.^(1/2).*H.^2.*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(2.*phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(2.*phi(i,j+1)).*sr(i,j+1)+2.*3.^(1/2).*E.*H.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+(-4).*3.^(1/2).*E.*H.*v.*cos(phi(i,j)+phi(i,j+1)).*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j).^2.*sr(i,j+1)+4.*3.^(1/2).*H.^2.*v.*(1+v).*phi(i,j+1).^2.*sr(i,j+1)+3.*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+(-6).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*sr(i,j+1)+H.*phi(i,j).*(8.*E.*H.*((-1)+2.*v)+(-1).*((-6).*E+12.*E.*v+3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.*v.*(1+v).*phi(i,j+1)+3.^(1/2).*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+4.*3.^(1/2).*H.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-2).*3.^(1/2).*E.*sin(phi(i,j)+phi(i,j+1))+4.*3.^(1/2).*E.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1))+(-1).*H.*phi(i,j+1).*(8.*E.*H.*((-1)+2.*v)+(-1).*(3.*E.*((-1)+2.*v).*cos(phi(i,j)+(-1).*phi(i,j+1))+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*sin(phi(i,j)+(-1).*phi(i,j+1))+2.*E.*((-1)+2.*v).*(3+3.^(1/2).*sin(phi(i,j)+phi(i,j+1)))).*sr(i,j+1)))).^(-1).*(exp(eth(i,j+1)).*r(i).*s(i,j+1).^2.*(8.*E.*H.*(3.^(1/2).*(E+H).*cos(phi(i,j))+(-1).*3.^(1/2).*H.*cos(phi(i,j+1))+H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sin(phi(i,j+1))+((-3).*E.*H+6.*E.*H.*v+(-3).*E.*H.*cos(phi(i,j)+(-1).*phi(i,j+1))+6.*E.*H.*v.*cos(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+3.*3.^(1/2).*E.*H.*sin(phi(i,j)+(-1).*phi(i,j+1))+12.*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.*sin(phi(i,j)+(-1).*phi(i,j+1))+8.*3.^(1/2).*H.^2.*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1))+H.*phi(i,j).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+(-1).*H.*phi(i,j+1).*((-4).*3.^(1/2).*H.*v+(-4).*3.^(1/2).*H.*v.^2+3.^(1/2).*(E+4.*E.*v+4.*H.*v.*(1+v)).*cos(phi(i,j)+(-1).*phi(i,j+1))+2.*3.^(1/2).*E.*((-1)+2.*v).*cos(phi(i,j)+phi(i,j+1))+3.*E.*sin(phi(i,j)+(-1).*phi(i,j+1))+(-6).*E.*v.*sin(phi(i,j)+(-1).*phi(i,j+1)))+2.*3.^(1/2).*E.*H.*sin(2.*phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(2.*phi(i,j+1))+2.*3.^(1/2).*E.*H.*sin(phi(i,j)+phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*v.*sin(phi(i,j)+phi(i,j+1))).*sr(i,j+1)).*dphi(i,j)+(-1).*s(i,j).*((-4).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*sth(i,j+1)+sr(i,j+1).*((-12).*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j))+6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+(-12).*E.*H.^2.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j))+4.*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-8).*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-12).*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+8.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*((-3).*H.*cos(phi(i,j))+3.*(E+H).*cos(phi(i,j+1))+(-1).*3.^(1/2).*H.*((-1)+2.*v).*(sin(phi(i,j))+(-1).*sin(phi(i,j+1)))).*deh(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*((-2).*H.*v.*cos(phi(i,j))+(E+2.*H.*v).*cos(phi(i,j+1))+3.^(1/2).*(2.*H.*v.*sin(phi(i,j))+(-1).*(E+2.*H.*v).*sin(phi(i,j+1)))).*deth(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-3).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*phi(i,j+1).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))+s(i,j+1).*(4.*E.*H.*((-1).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*(3.*E.*cos(phi(i,j+1))+6.*H.*cos(phi(i,j+1))+H.*phi(i,j+1).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1)))+(-2).*3.^(1/2).*H.*sin(phi(i,j+1))+4.*3.^(1/2).*H.*v.*sin(phi(i,j+1))+H.*phi(i,j).*((-1).*3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+3.*sin(phi(i,j+1)))).*sth(i,j+1)+(-2).*sin(phi(i,j+1)).*((-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H+2.*exp(er(i,j+1)).*E.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*(3.^(1/2).*E.*cos(phi(i,j))+2.*3.^(1/2).*H.*cos(phi(i,j))+(-2).*H.*sin(phi(i,j))+4.*H.*v.*sin(phi(i,j))+(-1).*H.*phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+H.*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j+1)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*H.*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*H.*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j+1)).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*sin(phi(i,j+1)).*ds(i,j)))+sr(i,j+1).*(6.*exp(eth(i,j+1)).*E.^2.*H.*cos(phi(i,j+1))+24.*exp(eth(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1))+12.*E.^2.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+24.*E.*H.^2.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-6).*exp(er(i,j+1)).*E.^2.*H.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-24).*exp(er(i,j+1)).*E.*H.^2.*v.*cos(phi(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*sin(phi(i,j+1))+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*v.*sin(phi(i,j+1))+(-8).*3.^(1/2).*E.*H.^2.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+16.*3.^(1/2).*E.*H.^2.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+6.*exp(er(i,j+1)).*E.^2.*H.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+24.*exp(er(i,j+1)).*E.*H.^2.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deh(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*cos(phi(i,j+1)).*deh(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*sin(phi(i,j+1)).*deh(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deh(i,j)+6.*exp(eth(i,j+1)).*E.^2.*H.*r(i).*cos(phi(i,j+1)).*deth(i,j)+24.*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*cos(phi(i,j+1)).*deth(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*H.*r(i).*sin(phi(i,j+1)).*deth(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.^2.*r(i).*v.*sin(phi(i,j+1)).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-6).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-24).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-16).*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*s(i,j).*sin(2.*phi(i,j+1)).*dphi(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.^2.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+3.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(2.*phi(i,j+1)).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*ds(i,j)+3.*exp(eth(i,j+1)).*E.*H.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-6).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+4.*3.^(1/2).*exp(eth(i,j+1)).*H.^2.*r(i).*v.*(1+v).*phi(i,j+1).^2.*(s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-1).*cos(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j))+H.*phi(i,j+1).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j))+(-1).*H.*phi(i,j).*((-12).*3.^(1/2).*exp(eth(i,j+1)).*E.*H.*v.*cos(phi(i,j+1))+(-4).*3.^(1/2).*E.*H.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+8.*3.^(1/2).*E.*H.*v.*cos(phi(i,j+1)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+12.*exp(er(i,j+1)).*E.*H.*v.*cos(phi(i,j+1)).*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-12).*exp(eth(i,j+1)).*E.*H.*v.*sin(phi(i,j+1))+(-12).*E.*H.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+12.*exp(er(i,j+1)).*E.*H.*v.*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j+1))+4.*exp(eth(i,j+1)).*E.*H.*r(i).*(3.^(1/2).*((-1)+2.*v).*cos(phi(i,j+1))+(-3).*sin(phi(i,j+1))).*deh(i,j)+(-12).*exp(eth(i,j+1)).*E.*H.*r(i).*v.*(3.^(1/2).*cos(phi(i,j+1))+sin(phi(i,j+1))).*deth(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+16.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*s(i,j).*dphi(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*cos(phi(i,j)+phi(i,j+1)).*s(i,j).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*phi(i,j+1).*s(i,j).*sin(phi(i,j)+(-1).*phi(i,j+1)).*dphi(i,j)+(-3).*exp(eth(i,j+1)).*E.*r(i).*ds(i,j)+6.*exp(eth(i,j+1)).*E.*r(i).*v.*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+(-8).*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*cos(phi(i,j)+(-1).*phi(i,j+1)).*phi(i,j+1).*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+12.*3.^(1/2).*exp(eth(i,j+1)).*H.*r(i).*v.^2.*sin(phi(i,j)+(-1).*phi(i,j+1)).*ds(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(2.*phi(i,j+1)).*ds(i,j)+4.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(2.*phi(i,j+1)).*ds(i,j)+(-4).*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)+8.*3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*v.*sin(phi(i,j)+phi(i,j+1)).*ds(i,j)))))));
s(i+1,j+1) = s(i,j+1) + ds(i,j+1)*dr;
phi(i+1,j+1) = phi(i,j+1) + dphi(i,j+1)*dr;
else
re(i,j)=r(i);
er(i,j+1) = er(i,j)+((-1).*s(i,j)+s(i,j+1)).*(E.^(-1).*(2.* ...
3.^(-1/2)+(-1).*3.^(-1/2).*v).*cos(phi(i,j))+(-1).*E.^(-1).*v.* ...
sin(phi(i,j)))+((-1).*phi(i,j)+phi(i,j+1)).*((-1).*E.^(-1).*v.* ...
cos(phi(i,j)).*s(i,j)+E.^(-1).*((-2).*3.^(-1/2).*s(i,j)+3.^(-1/2) ...
.*v.*s(i,j)).*sin(phi(i,j)));
eth(i,j+1) = eth(i,j)+((-1).*s(i,j)+s(i,j+1)).*(E.^(-1).*( ...
3.^(-1/2)+(-2).*3.^(-1/2).*v).*cos(phi(i,j))+E.^(-1).*sin(phi(i,j)))+((-1).*phi(i,j)+phi(i,j+1)).*(E.^(-1).*cos(phi(i,j)).*s(i,j)+E.^(-1).*((-1).*3.^(-1/2).*s(i,j)+2.*3.^(-1/2).*v.*s(i,j)).* ...
sin(phi(i,j)));
eh(i,j+1) = E.^(-1).*(E.*eh(i,j)+3.^(1/2).*v.*cos(phi(i,j)).* ...
s(i,j)+v.*cos(phi(i,j)).*phi(i,j).*s(i,j)+(-1).*v.*cos(phi(i,j)).* ...
phi(i,j+1).*s(i,j)+(-1).*3.^(1/2).*v.*cos(phi(i,j)).*s(i,j+1)+ ...
v.*s(i,j).*sin(phi(i,j))+(-1).*3.^(1/2).*v.*phi(i,j).*s(i,j).*sin( ...
phi(i,j))+3.^(1/2).*v.*phi(i,j+1).*s(i,j).*sin(phi(i,j))+(-1).* ...
v.*s(i,j+1).*sin(phi(i,j)));
e(i,j+1) = sqrt(2)/3*sqrt( (er(i,j+1)-eth(i,j+1)).^2 + (er(i,j+1)-eh(i,j+1)).^2 + (eh(i,j+1)-eth(i,j+1)).^2 );
dw(i,j+1) = (-1/2).*exp(-eh(i,j+1)+er(i,j+1)+eth(i,j+1)).*ho.^(-1).*p.*r(i).*sr(i,j+1).^(-1);
ds(i,j+1) = ((1/2).*exp((-1).*eth(i,j+1)).*r(i).^(-1).*((-1).*E.*s(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j))).*sin(phi(i,j+1))+s(i,j).*(E.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*cos(phi(i,j))+(1+(-2).*v).*sin(phi(i,j)))+3.^(1/2).*v.*(1+v).*sr(i,j+1))).^(-1).*(2.*exp(eth(i,j+1)).*E.*r(i).*s(i,j+1).^2.*(3.^(1/2).*cos(phi(i,j))+((-1)+2.*v).*sin(phi(i,j))).*sin(phi(i,j+1)).*dphi(i,j)+2.*E.*s(i,j+1).*sin(phi(i,j+1)).*(3.^(1/2).*exp(eth(i,j+1)).*E+(-1).*exp(er(i,j+1)).*E.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*((-2).*3.^(1/2).*cos(phi(i,j))+2.*sin(phi(i,j))+(-4).*v.*sin(phi(i,j))+phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+(-1).*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+(-1).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*r(i).*v.*cos(phi(i,j)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j))+s(i,j).*(E.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*cos(phi(i,j))+3.^(1/2).*(1+(-2).*v).*sin(phi(i,j))).*sth(i,j+1)+sr(i,j+1).*(3.*exp(eth(i,j+1)).*E.*v.*cos(phi(i,j))+3.*E.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-3).*exp(er(i,j+1)).*E.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*v.*sin(phi(i,j))+(-1).*3.^(1/2).*E.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+2.*3.^(1/2).*E.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+3.*exp(er(i,j+1)).*E.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+exp(eth(i,j+1)).*E.*r(i).*(3.*cos(phi(i,j))+3.^(1/2).*((-1)+2.*v).*sin(phi(i,j))).*deh(i,j)+3.*exp(eth(i,j+1)).*E.*r(i).*v.*(cos(phi(i,j))+(-1).*3.^(1/2).*sin(phi(i,j))).*deth(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*ds(i,j)))));
dphi(i,j+1) = ((-3).*E.*r(i).^(-1).*s(i,j).^(-1).*(3.*cos(phi(i,j))+(-1).*3.^(1/2).*sin(phi(i,j))+2.*3.^(1/2).*v.*sin(phi(i,j))).^(-1).*(1+(-1).*exp(er(i,j+1)+(-1).*eth(i,j+1)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+r(i).*deth(i,j)+(-1/3).*E.^(-1).*r(i).*s(i,j).*(3.*cos(phi(i,j))+(-1).*3.^(1/2).*sin(phi(i,j))+2.*3.^(1/2).*v.*sin(phi(i,j))).*dphi(i,j)+(-1/3).*E.^(-1).*r(i).*(s(i,j)+(-1).*s(i,j+1)).*(3.*cos(phi(i,j))+(-1).*3.^(1/2).*sin(phi(i,j))+2.*3.^(1/2).*v.*sin(phi(i,j))).*dphi(i,j)+(-1/3).*E.^(-1).*r(i).*(phi(i,j)+(-1).*phi(i,j+1)).*((-1).*3.^(1/2).*cos(phi(i,j)).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*v.*cos(phi(i,j)).*s(i,j).*dphi(i,j)+(-3).*s(i,j).*sin(phi(i,j)).*dphi(i,j)+3.*cos(phi(i,j)).*ds(i,j)+(-1).*3.^(1/2).*sin(phi(i,j)).*ds(i,j)+2.*3.^(1/2).*v.*sin(phi(i,j)).*ds(i,j))+(1/3).*E.^(-1).*r(i).*((-1).*3.^(1/2).*cos(phi(i,j))+2.*3.^(1/2).*v.*cos(phi(i,j))+(-3).*sin(phi(i,j))).*(ds(i,j)+(-1/2).*exp((-1).*eth(i,j+1)).*r(i).^(-1).*((-1).*E.*s(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j))).*sin(phi(i,j+1))+s(i,j).*(E.*cos(phi(i,j+1)).*((-1).*3.^(1/2).*cos(phi(i,j))+(1+(-2).*v).*sin(phi(i,j)))+3.^(1/2).*v.*(1+v).*sr(i,j+1))).^(-1).*(2.*exp(eth(i,j+1)).*E.*r(i).*s(i,j+1).^2.*(3.^(1/2).*cos(phi(i,j))+((-1)+2.*v).*sin(phi(i,j))).*sin(phi(i,j+1)).*dphi(i,j)+2.*E.*s(i,j+1).*sin(phi(i,j+1)).*(3.^(1/2).*exp(eth(i,j+1)).*E+(-1).*exp(er(i,j+1)).*E.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+3.^(1/2).*exp(eth(i,j+1)).*E.*r(i).*deth(i,j)+exp(eth(i,j+1)).*r(i).*s(i,j).*((-2).*3.^(1/2).*cos(phi(i,j))+2.*sin(phi(i,j))+(-4).*v.*sin(phi(i,j))+phi(i,j).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))+(-1).*phi(i,j+1).*((1+(-2).*v).*cos(phi(i,j))+3.^(1/2).*sin(phi(i,j)))).*dphi(i,j)+(-1).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*r(i).*v.*cos(phi(i,j)).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*phi(i,j).*ds(i,j)+3.^(1/2).*exp(eth(i,j+1)).*r(i).*cos(phi(i,j)).*phi(i,j+1).*ds(i,j)+(-1).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*sin(phi(i,j)).*ds(i,j)+exp(eth(i,j+1)).*r(i).*phi(i,j).*sin(phi(i,j)).*ds(i,j)+(-2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j).*sin(phi(i,j)).*ds(i,j)+(-1).*exp(eth(i,j+1)).*r(i).*phi(i,j+1).*sin(phi(i,j)).*ds(i,j)+2.*exp(eth(i,j+1)).*r(i).*v.*phi(i,j+1).*sin(phi(i,j)).*ds(i,j))+s(i,j).*(E.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*((-3).*cos(phi(i,j))+3.^(1/2).*(1+(-2).*v).*sin(phi(i,j))).*sth(i,j+1)+sr(i,j+1).*(3.*exp(eth(i,j+1)).*E.*v.*cos(phi(i,j))+3.*E.*cos(phi(i,j)).*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2)+(-3).*exp(er(i,j+1)).*E.*v.*cos(phi(i,j)).*(1+(-1).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2)+(-3).*3.^(1/2).*exp(eth(i,j+1)).*E.*v.*sin(phi(i,j))+(-1).*3.^(1/2).*E.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+2.*3.^(1/2).*E.*v.*(exp((2).*er(i,j+1))+(-1).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+3.*exp(er(i,j+1)).*E.*v.*(3+(-3).*exp((-2).*er(i,j+1)).*dw(i,j+1).^2).^(1/2).*sin(phi(i,j))+exp(eth(i,j+1)).*E.*r(i).*(3.*cos(phi(i,j))+3.^(1/2).*((-1)+2.*v).*sin(phi(i,j))).*deh(i,j)+3.*exp(eth(i,j+1)).*E.*r(i).*v.*(cos(phi(i,j))+(-1).*3.^(1/2).*sin(phi(i,j))).*deth(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j).*s(i,j).*dphi(i,j)+(-2).*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*phi(i,j).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*phi(i,j+1).*s(i,j).*dphi(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.*ds(i,j)+2.*3.^(1/2).*exp(eth(i,j+1)).*r(i).*v.^2.*ds(i,j)))))));
s(i+1,j+1) = s(i,j+1) + ds(i,j+1)*dr;
phi(i+1,j+1) = phi(i,j+1) + dphi(i,j+1)*dr;
end
if ~isreal(s(i+1,j+1))
display('Complex number occured')
break
end
end
sr(i+1,j+1) = (2.*3.^(-1/2).*cos(phi(i+1,j+1)).*s(i+1,j+1));
sth(i+1,j+1) = s(i+1,j+1).*(1/sqrt(3)*cos(phi(i+1,j+1)) + sin(phi(i+1,j+1)));
er(i+1,j+1) = er(i+1,j)+((-1).*s(i+1,j)+s(i+1,j+1)).*(E.^(-1).*(2.* ...
3.^(-1/2)+(-1).*3.^(-1/2).*v).*cos(phi(i+1,j))+(-1).*E.^(-1).*v.* ...
sin(phi(i+1,j)))+((-1).*phi(i+1,j)+phi(i+1,j+1)).*((-1).*E.^(-1).*v.* ...
cos(phi(i+1,j)).*s(i+1,j)+E.^(-1).*((-2).*3.^(-1/2).*s(i+1,j)+3.^(-1/2) ...
.*v.*s(i+1,j)).*sin(phi(i+1,j)));
eth(i+1,j+1) = eth(i+1,j)+((-1).*s(i+1,j)+s(i+1,j+1)).*(E.^(-1).*( ...
3.^(-1/2)+(-2).*3.^(-1/2).*v).*cos(phi(i+1,j))+E.^(-1).*sin(phi(i+1,j)))+((-1).*phi(i+1,j)+phi(i+1,j+1)).*(E.^(-1).*cos(phi(i+1,j)).*s(i+1,j)+E.^(-1).*((-1).*3.^(-1/2).*s(i+1,j)+2.*3.^(-1/2).*v.*s(i+1,j)).* ...
sin(phi(i+1,j)));
eh(i+1,j+1) = E.^(-1).*(E.*eh(i+1,j)+3.^(1/2).*v.*cos(phi(i+1,j)).* ...
s(i+1,j)+v.*cos(phi(i+1,j)).*phi(i+1,j).*s(i+1,j)+(-1).*v.*cos(phi(i+1,j)).* ...
phi(i+1,j+1).*s(i+1,j)+(-1).*3.^(1/2).*v.*cos(phi(i+1,j)).*s(i+1,j+1)+ ...
v.*s(i+1,j).*sin(phi(i+1,j))+(-1).*3.^(1/2).*v.*phi(i+1,j).*s(i+1,j).*sin( ...
phi(i+1,j))+3.^(1/2).*v.*phi(i+1,j+1).*s(i+1,j).*sin(phi(i+1,j))+(-1).* ...
v.*s(i+1,j+1).*sin(phi(i+1,j)));
e(i+1,j+1) = sqrt(2)/3*sqrt( (er(i+1,j+1)-eth(i+1,j+1)).^2 + (er(i+1,j+1)-eh(i+1,j+1)).^2 + (eh(i+1,j+1)-eth(i+1,j+1)).^2 );
dw(i+1,j+1) = (-1/2).*exp(-eh(i+1,j+1)+er(i+1,j+1)+eth(i+1,j+1)).*ho.^(-1).*p.*r(i+1).*sr(i+1,j+1).^(-1);
ds(i+1,j+1) = ds(i,j+1);
dphi(i+1,j+1) = dphi(i,j+1);
disp(['j: ' num2str(j) ' --- Iterations: ' num2str(iter) ' --- e_o/e_r = ' num2str(eth(1,j+1)/eth(end-1,j+1)) ' --- P = ' num2str(p)] )
if eth(end-1,j+1) < 0
p = p - p_inc;
elseif eth(end-1,j+1) > eth(1,j+1)*1e-2 || ~isreal(er(:,j+1))
p = p + p_inc;
end
iter = iter + 1;
if iter > 10000
disp(['Cannot converge at ' num2str(s(1,j+1)/1e6) ' MPa. Last time step not correct.'])
break;
end
end
w(:,j+1) = cumtrapz(r,dw(:,j+1)) - trapz(r,dw(:,j+1));
pj(j+1)=p;
deth(1:length(r),j+1) = gradient(eth(:,j+1),r);
deh(1:length(r),j+1) = gradient(eh(:,j+1),r);
de(1:length(r),j+1) = gradient(e(:,j+1),r);
end
figure;
plot(er(1,:),sr(1,:)/1e6)
xlabel('$\varepsilon_0$','Interpreter','Latex')
ylabel('$\bar{\sigma}_0$ (MPa)','Interpreter','Latex')
set(gca,'FontName','Palatino Linotype','FontSize',12)
figure;
plot(pj*ro/E/ho,w(1,:)/ro)
xlabel('$P r_0/E h_0$','Interpreter','Latex')
ylabel('$w/r_0$','Interpreter','Latex')
set(gca,'FontName','Palatino Linotype','FontSize',12)
figure;
plot(pj*ro/E/ho,er(1,:))
xlabel('$P r_0/E h_0$','Interpreter','Latex')
ylabel('$\varepsilon_0$','Interpreter','Latex')
set(gca,'FontName','Palatino Linotype','FontSize',12)
figure;
plot(pj*ro/E/ho,sr(1,:)/1e6)
xlabel('$P r_0/E h_0$','Interpreter','Latex')
ylabel('$\bar{\sigma}_0$ (MPa)','Interpreter','Latex')
set(gca,'FontName','Palatino Linotype','FontSize',12)
toc