-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathTwinDisks.m
922 lines (881 loc) · 114 KB
/
TwinDisks.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
function TwinDisks
SolveOrdinaryDifferentialEquations
% File TwinDisks.m created by Autolev 4.1 on Wed Sep 22 17:32:55 2010
%===========================================================================
function VAR = ReadUserInput
global i11 i22 i33 l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p u1p a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55;
global DEGtoRAD RADtoDEG z;
global TINITIAL TFINAL INTEGSTP PRINTINT ABSERR RELERR;
%-------------------------------+--------------------------+-------------------+-----------------
% Quantity | Value | Units | Description
%-------------------------------|--------------------------|-------------------|-----------------
i11 = 0.005000000000000001; % UNITS Constant
i22 = 0.005000000000000001; % UNITS Constant
i33 = 0.01; % UNITS Constant
l = 1; % UNITS Constant
m = 2; % UNITS Constant
rad = .1; % UNITS Constant
cahat1 = 0; % UNITS Initial Value
cahat2 = 0; % UNITS Initial Value
cahat3 = 0; % UNITS Initial Value
cbhat1 = -0.1732050807568877; % UNITS Initial Value
cbhat2 = 0; % UNITS Initial Value
cbhat3 = 0; % UNITS Initial Value
e1 = 0.3826834323650898; % UNITS Initial Value
e2 = 0; % UNITS Initial Value
e3 = 0; % UNITS Initial Value
e4 = 0.9238795325112867; % UNITS Initial Value
posa1 = -0.08660254037844388; % UNITS Initial Value
posa2 = 0; % UNITS Initial Value
posa3 = 0.05000000000000002; % UNITS Initial Value
posb1 = -0.1732050807568877; % UNITS Initial Value
posb2 = 0; % UNITS Initial Value
posb3 = 0.1; % UNITS Initial Value
u1 = 3; % UNITS Initial Value
TINITIAL = 0.0; % UNITS Initial Time
TFINAL = 1.0; % UNITS Final Time
INTEGSTP = 0.01; % UNITS Integration Step
PRINTINT = 1; % Positive Integer Print-Integer
ABSERR = 1.0E-08; % Absolute Error
RELERR = 1.0E-07 ; % Relative Error
%-------------------------------+--------------------------+-------------------+-----------------
% Unit conversions
Pi = 3.141592653589793;
DEGtoRAD = Pi/180.0;
RADtoDEG = 180.0/Pi;
% Reserve space and initialize matrices
z = zeros(618,1);
% Evaluate constants
z(22) = l*rad;
cahat3p = 0;
cbhat3p = 0;
z(36) = z(22)^2;
z(130) = i22 + i33;
z(134) = m*z(22);
z(213) = rad*z(22);
z(293) = rad^2;
% Set the initial values of the states
VAR(1) = cahat1;
VAR(2) = cahat2;
VAR(3) = cahat3;
VAR(4) = cbhat1;
VAR(5) = cbhat2;
VAR(6) = cbhat3;
VAR(7) = e1;
VAR(8) = e2;
VAR(9) = e3;
VAR(10) = e4;
VAR(11) = posa1;
VAR(12) = posa2;
VAR(13) = posa3;
VAR(14) = posb1;
VAR(15) = posb2;
VAR(16) = posb3;
VAR(17) = u1;
%===========================================================================
function OpenOutputFilesAndWriteHeadings
FileIdentifier = fopen('TwinDisks.1', 'wt'); if( FileIdentifier == -1 ) error('Error: unable to open file TwinDisks.1'); end
fprintf( 1, '%% t ke pe ke+pe dot(v_bhat_n>, dot(v_bhat_n>, e1 e2 e3 e4 u1 cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 dot(w_diska_n> dot(w_diska_n> dot(w_diska_n> a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55\n' );
fprintf( 1, '%% (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS)\n\n' );
fprintf(FileIdentifier, '%% FILE: TwinDisks.1\n%%\n' );
fprintf(FileIdentifier, '%% t ke pe ke+pe dot(v_bhat_n>, dot(v_bhat_n>, e1 e2 e3 e4 u1 cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 dot(w_diska_n> dot(w_diska_n> dot(w_diska_n> a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55\n' );
fprintf(FileIdentifier, '%% (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS) (UNITS)\n\n' );
%===========================================================================
% Main driver loop for numerical integration of differential equations
%===========================================================================
function SolveOrdinaryDifferentialEquations
global i11 i22 i33 l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p u1p a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55;
global DEGtoRAD RADtoDEG z;
global TINITIAL TFINAL INTEGSTP PRINTINT ABSERR RELERR;
OpenOutputFilesAndWriteHeadings
VAR = ReadUserInput;
OdeMatlabOptions = odeset( 'RelTol',RELERR, 'AbsTol',ABSERR, 'MaxStep',INTEGSTP );
T = TINITIAL;
PrintCounter = 0;
mdlDerivatives(T,VAR,0);
while 1,
if( TFINAL>=TINITIAL & T+0.01*INTEGSTP>=TFINAL ) PrintCounter = -1; end
if( TFINAL<=TINITIAL & T+0.01*INTEGSTP<=TFINAL ) PrintCounter = -1; end
if( PrintCounter <= 0.01 ),
mdlOutputs(T,VAR,0);
if( PrintCounter == -1 ) break; end
PrintCounter = PRINTINT;
end
[TimeOdeArray,VarOdeArray] = ode45( @mdlDerivatives, [T T+INTEGSTP], VAR, OdeMatlabOptions, 0 );
TimeAtEndOfArray = TimeOdeArray( length(TimeOdeArray) );
if( abs(TimeAtEndOfArray - (T+INTEGSTP) ) >= abs(0.001*INTEGSTP) ) warning('numerical integration failed'); break; end
T = TimeAtEndOfArray;
VAR = VarOdeArray( length(TimeOdeArray), : );
PrintCounter = PrintCounter - 1;
end
mdlTerminate(T,VAR,0);
data = load( 'TwinDisks.1' );
figure(1), plot(data(:,1),data(:,2),data(:,1),data(:,3),data(:,1),data(:,4))
figure(2), plot(data(:,1),data(:,18),data(:,1),data(:,19),data(:,1),data(:,20))
figure(3), plot(data(:,1),data(:,7),data(:,1),data(:,8),data(:,1),data(:,9),data(:,1),data(:,10),data(:,1),data(:,11))
%===========================================================================
% mdlDerivatives: Calculates and returns the derivatives of the continuous states
%===========================================================================
function sys = mdlDerivatives(T,VAR,u)
global i11 i22 i33 l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p u1p a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55;
global DEGtoRAD RADtoDEG z;
global TINITIAL TFINAL INTEGSTP PRINTINT ABSERR RELERR;
% Update variables after integration step
cahat1 = VAR(1);
cahat2 = VAR(2);
cahat3 = VAR(3);
cbhat1 = VAR(4);
cbhat2 = VAR(5);
cbhat3 = VAR(6);
e1 = VAR(7);
e2 = VAR(8);
e3 = VAR(9);
e4 = VAR(10);
posa1 = VAR(11);
posa2 = VAR(12);
posa3 = VAR(13);
posb1 = VAR(14);
posb2 = VAR(15);
posb3 = VAR(16);
u1 = VAR(17);
z(2) = 2*e1*e2 - 2*e3*e4;
z(4) = 2*e1*e2 + 2*e3*e4;
z(9) = 1 - 2*e1^2 - 2*e2^2;
z(10) = 1 - z(9)^2;
z(11) = z(10)^0.5;
z(12) = z(9)/z(11);
z(14) = rad*z(12);
z(8) = 2*e1*e4 + 2*e2*e3;
z(16) = 1 - z(8)^2;
z(17) = z(16)^0.5;
z(18) = z(8)/z(17);
z(20) = rad*z(18);
z(13) = 1/z(11);
z(15) = rad*z(13);
z(19) = 1/z(17);
z(21) = rad*z(19);
z(35) = z(15) - z(21);
z(7) = 2*e1*e3 - 2*e2*e4;
z(37) = z(36) + z(14)^2 + z(20)^2 + z(35)^2 + 2*z(8)*z(20)*z(35) - 2*z(22)*z(7)*z(35) - 2*z(9)*z(14)*z(35);
z(38) = z(37)^0.5;
z(41) = z(22)/z(38);
z(6) = 2*e2*e3 - 2*e1*e4;
z(39) = z(14)/z(38);
z(42) = z(20)/z(38);
z(3) = 2*e1*e3 + 2*e2*e4;
z(25) = z(6)*z(15);
z(31) = z(25) - z(6)*z(21);
z(5) = 1 - 2*e1^2 - 2*e3^2;
z(1) = 1 - 2*e2^2 - 2*e3^2;
z(28) = z(3)*z(15);
z(34) = z(3)*z(21) - z(28);
z(44) = z(22)*z(2)*(z(4)*z(41)+z(6)*z(39)) + z(4)*z(20)*(z(2)*z(42)-z(3)*z(39)) + z(31)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) - z(22)*z(5)*(z(1)*z(41)+z(3)*z(39)) - z(1)*z(20)*(z(5)*z(42)-z(6)*z(39)) - z(34)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39));
z(47) = z(8)*z(14) + z(9)*z(20);
z(23) = z(4)*z(15);
z(29) = z(23) - z(4)*z(21);
z(26) = z(1)*z(15);
z(32) = z(1)*z(21) - z(26);
z(45) = z(5)*z(14)*(z(1)*z(41)+z(3)*z(39)) + z(6)*z(20)*(z(1)*z(41)-z(2)*z(42)) + z(29)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) - z(2)*z(14)*(z(4)*z(41)+z(6)*z(39)) - z(3)*z(20)*(z(4)*z(41)-z(5)*z(42)) - z(32)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39));
z(48) = -z(22)*z(8) - z(7)*z(20);
z(24) = z(5)*z(15);
z(30) = z(24) - z(5)*z(21);
z(27) = z(2)*z(15);
z(33) = z(2)*z(21) - z(27);
z(43) = z(22)*z(6)*(z(1)*z(41)-z(2)*z(42)) + z(4)*z(14)*(z(2)*z(42)-z(3)*z(39)) + z(30)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) - z(22)*z(3)*(z(4)*z(41)-z(5)*z(42)) - z(1)*z(14)*(z(5)*z(42)-z(6)*z(39)) - z(33)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39));
z(46) = z(22)*z(9) - z(7)*z(14);
z(49) = z(43)*z(48) - z(44)*z(46);
z(50) = (z(44)*z(47)-z(45)*z(48))/z(49);
u2 = z(50)*u1;
z(51) = (z(43)*z(47)-z(45)*z(46))/z(49);
u3 = -z(51)*u1;
z(101) = i33*u3;
z(100) = i22*u2;
z(104) = u2*z(101) - u3*z(100);
z(106) = i22*u3;
z(105) = i33*u2;
z(109) = u2*z(106) - u3*z(105);
z(111) = z(14)*u1;
z(110) = z(14)*u2;
z(119) = u1*z(111) + u2*z(110);
e1p = 0.5*e2*u3 + 0.5*e4*u1 - 0.5*e3*u2;
e2p = 0.5*e3*u1 + 0.5*e4*u2 - 0.5*e1*u3;
z(54) = -4*e1*e1p - 4*e2*e2p;
z(55) = z(9)*z(54)/z(10)^0.5;
z(56) = (z(9)*z(55)+z(11)*z(54))/z(11)^2;
z(115) = rad*u1*z(56);
z(120) = z(115) - u3*z(110);
z(60) = rad*z(55)/z(11)^2;
e4p = -0.5*e1*u1 - 0.5*e2*u2 - 0.5*e3*u3;
e3p = 0.5*e1*u2 + 0.5*e4*u3 - 0.5*e2*u1;
z(67) = 2*e1*e2p + 2*e2*e1p + 2*e3*e4p + 2*e4*e3p;
z(81) = z(4)*z(60) + z(15)*z(67);
z(70) = -4*e1*e1p - 4*e3*e3p;
z(71) = z(5)*z(60) + z(15)*z(70);
z(52) = 2*e2*e3p + 2*e3*e2p - 2*e1*e4p - 2*e4*e1p;
z(76) = z(6)*z(60) + z(15)*z(52);
z(116) = u1*z(81) + u2*z(71) + u3*z(76);
z(114) = rad*u2*z(56);
z(118) = -u3*z(111) - z(114);
z(53) = -4*e2*e2p - 4*e3*e3p;
z(83) = z(1)*z(60) + z(15)*z(53);
z(65) = 2*e1*e2p + 2*e2*e1p - 2*e3*e4p - 2*e4*e3p;
z(73) = z(2)*z(60) + z(15)*z(65);
z(68) = 2*e1*e3p + 2*e2*e4p + 2*e3*e1p + 2*e4*e2p;
z(78) = z(3)*z(60) + z(15)*z(68);
z(117) = -u1*z(83) - u2*z(73) - u3*z(78);
z(121) = z(14)*u1 - z(22)*u3;
z(124) = u1*z(121) + u2*z(110);
z(122) = z(22)*u2;
z(125) = z(115) - u1*z(122) - u3*z(110);
z(123) = u2*z(122) - u3*z(121) - z(114);
z(129) = z(104) + z(109) + m*(z(3)*z(23)*z(119)+z(14)*z(120)+2*z(23)*z(116)+z(1)*z(23)*z(118)+z(2)*z(23)*z(120)+2*z(5)*z(14)*z(117)-z(6)*z(26)*z(119)-z(4)*z(26)*z(118)-z(5)*z(26)*z(120)) + m*(z(3)*z(23)*z(124)+z(14)*z(125)+z(1)*z(23)*z(123)+z(2)*z(23)*z(125)+2*z(2)*z(14)*z(116)-z(6)*z(26)*z(124)-2*z(26)*z(117)-z(4)*z(26)*z(123)-z(5)*z(26)*z(125));
z(99) = i11*u1;
z(103) = u3*z(99) - u1*z(101);
z(108) = u1*z(106) - u3*z(99);
z(133) = z(103) + m*(z(3)*z(24)*z(119)+2*z(24)*z(116)+z(1)*z(24)*z(118)+z(2)*z(24)*z(120)-z(6)*z(27)*z(119)-2*z(27)*z(117)-z(14)*z(118)-2*z(1)*z(14)*z(116)-z(4)*z(27)*z(118)-z(5)*z(27)*z(120)) + m*(z(22)*z(124)+z(3)*z(24)*z(124)+z(22)*z(3)*z(116)+z(22)*z(6)*z(117)+z(1)*z(24)*z(123)+z(2)*z(24)*z(125)-z(6)*z(27)*z(124)-z(14)*z(123)-2*z(4)*z(14)*z(117)-z(4)*z(27)*z(123)-z(5)*z(27)*z(125)) - z(108);
z(102) = u1*z(100) - u2*z(99);
z(107) = u1*z(105) - u2*z(99);
z(136) = z(102) + z(107) + m*(z(3)*z(25)*z(119)+2*z(25)*z(116)+z(1)*z(25)*z(118)+z(2)*z(25)*z(120)-z(6)*z(28)*z(119)-2*z(28)*z(117)-z(4)*z(28)*z(118)-z(5)*z(28)*z(120)) + m*(z(3)*z(25)*z(124)+z(1)*z(25)*z(123)+z(2)*z(25)*z(125)-z(6)*z(28)*z(124)-z(22)*z(125)-z(22)*z(2)*z(116)-z(22)*z(5)*z(117)-z(4)*z(28)*z(123)-z(5)*z(28)*z(125));
z(140) = z(129) + z(50)*z(133) - z(51)*z(136);
z(127) = m*(2*z(23)*z(24)+2*z(26)*z(27)+z(22)*z(3)*z(23)+2*z(2)*z(14)*z(24)+2*z(4)*z(14)*z(26)-2*z(1)*z(14)*z(23)-2*z(5)*z(14)*z(27)-z(22)*z(6)*z(26));
z(131) = z(130) + 2*m*(z(14)^2+z(24)^2+z(27)^2-2*z(1)*z(14)*z(24)) - m*(2*z(22)*z(6)*z(27)-z(36)-4*z(4)*z(14)*z(27)-2*z(22)*z(3)*z(24));
z(132) = m*(z(22)*z(2)*z(24)+z(22)*z(6)*z(28)+2*z(1)*z(14)*z(25)-2*z(24)*z(25)-2*z(27)*z(28)-2*z(4)*z(14)*z(28)-z(22)*z(3)*z(25)-z(22)*z(5)*z(27));
z(138) = z(127) + z(50)*z(131) + z(51)*z(132);
z(57) = 2*e1*e4p + 2*e2*e3p + 2*e3*e2p + 2*e4*e1p;
z(58) = z(8)*z(57)/z(16)^0.5;
z(59) = (z(8)*z(58)+z(17)*z(57))/z(17)^2;
z(88) = z(14)*z(57) + z(20)*z(54) + rad*z(8)*z(56) + rad*z(9)*z(59);
z(62) = 2*e1*e3p + 2*e3*e1p - 2*e2*e4p - 2*e4*e2p;
z(87) = z(22)*z(54) - z(14)*z(62) - rad*z(7)*z(56);
z(89) = -z(22)*z(57) - z(20)*z(62) - rad*z(7)*z(59);
z(90) = u1*z(88) + u2*z(87) + u3*z(89);
z(61) = rad*z(58)/z(17)^2;
z(82) = z(81) - z(4)*z(61) - z(21)*z(67);
z(63) = 2*rad*z(14)*z(56) + 2*rad*z(20)*z(59) + 2*z(20)*z(35)*z(57) + 2*rad*z(8)*z(35)*z(59) + 2*z(35)*(z(60)-z(61)) + 2*z(8)*z(20)*(z(60)-z(61)) - 2*z(22)*z(35)*z(62) - 2*z(14)*z(35)*z(54) - 2*rad*z(9)*z(35)*z(56) - 2*z(22)*z(7)*(z(60)-z(61)) - 2*z(9)*z(14)*(z(60)-z(61));
z(64) = z(22)*z(63)/(z(37)^0.5*z(38)^2);
z(69) = (2*rad*z(38)*z(56)-z(14)*z(63)/z(37)^0.5)/z(38)^2;
z(66) = (2*rad*z(38)*z(59)-z(20)*z(63)/z(37)^0.5)/z(38)^2;
z(84) = z(1)*z(61) + z(21)*z(53) - z(83);
z(85) = z(14)*(z(1)*z(41)+z(3)*z(39))*z(70) + rad*z(5)*(z(1)*z(41)+z(3)*z(39))*z(56) + z(20)*(z(1)*z(41)-z(2)*z(42))*z(52) + rad*z(6)*(z(1)*z(41)-z(2)*z(42))*z(59) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(82) + 0.5*z(2)*z(14)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) + 0.5*z(6)*z(20)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*z(32)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - z(14)*(z(4)*z(41)+z(6)*z(39))*z(65) - rad*z(2)*(z(4)*z(41)+z(6)*z(39))*z(56) - z(20)*(z(4)*z(41)-z(5)*z(42))*z(68) - rad*z(3)*(z(4)*z(41)-z(5)*z(42))*z(59) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(84) - 0.5*z(5)*z(14)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) - 0.5*z(3)*z(20)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(29)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66));
z(72) = z(71) - z(5)*z(61) - z(21)*z(70);
z(74) = z(2)*z(61) + z(21)*z(65) - z(73);
z(75) = z(22)*(z(1)*z(41)-z(2)*z(42))*z(52) + z(14)*(z(2)*z(42)-z(3)*z(39))*z(67) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(56) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(72) + 0.5*z(1)*z(14)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(22)*z(6)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*z(33)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - z(22)*(z(4)*z(41)-z(5)*z(42))*z(68) - z(14)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(56) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(74) - 0.5*z(4)*z(14)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(22)*z(3)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(30)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66));
z(77) = z(76) - z(6)*z(61) - z(21)*z(52);
z(79) = z(3)*z(61) + z(21)*z(68) - z(78);
z(80) = z(22)*(z(4)*z(41)+z(6)*z(39))*z(65) + z(20)*(z(2)*z(42)-z(3)*z(39))*z(67) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(59) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(77) + 0.5*z(22)*z(5)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) + 0.5*z(1)*z(20)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(34)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - z(22)*(z(1)*z(41)+z(3)*z(39))*z(70) - z(20)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(59) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(79) - 0.5*z(22)*z(2)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - 0.5*z(4)*z(20)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(31)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66));
z(86) = u1*z(85) + u2*z(75) + u3*z(80);
z(91) = (z(44)*z(90)-z(48)*z(86))/z(49);
z(128) = m*(2*z(23)*z(25)+2*z(26)*z(28)+z(22)*z(5)*z(26)+2*z(2)*z(14)*z(25)-z(22)*z(14)-2*z(5)*z(14)*z(28)-z(22)*z(2)*z(23));
z(135) = z(130) + 2*m*(z(25)^2+z(28)^2) - z(134)*(2*z(2)*z(25)-z(22)-2*z(5)*z(28));
z(139) = z(128) - z(50)*z(132) - z(51)*z(135);
z(92) = (z(43)*z(90)-z(46)*z(86))/z(49);
z(142) = z(140) + z(138)*z(91) - z(139)*z(92);
z(95) = z(14)*z(50);
z(96) = z(14) + z(22)*z(51);
z(97) = z(22)*z(50);
z(98) = m*(2*z(7)*z(95)-z(8)*z(14)-z(8)*z(96)-z(9)*z(97));
z(143) = z(142) - 9.810000000000001*z(98);
z(126) = 2*i11 + 2*m*(z(14)^2+z(23)^2+z(26)^2+2*z(2)*z(14)*z(23)-2*z(5)*z(14)*z(26));
z(137) = z(126) + z(50)*z(127) - z(51)*z(128);
z(141) = z(137) + z(50)*z(138) - z(51)*z(139);
z(144) = z(143)/z(141);
u1p = -z(144);
posa1p = z(23)*u1 + z(24)*u2 + z(25)*u3 + z(2)*z(14)*u1 - z(1)*z(14)*u2;
posa2p = z(5)*z(14)*u1 - z(26)*u1 - z(27)*u2 - z(28)*u3 - z(4)*z(14)*u2;
posa3p = -z(14)*(z(7)*u2-z(8)*u1);
posb1p = z(23)*u1 + z(24)*u2 + z(25)*u3 + z(22)*z(3)*u2 - z(1)*z(14)*u2 - z(2)*(z(22)*u3-z(14)*u1);
posb2p = z(22)*z(6)*u2 - z(26)*u1 - z(27)*u2 - z(28)*u3 - z(4)*z(14)*u2 - z(5)*(z(22)*u3-z(14)*u1);
posb3p = z(22)*z(9)*u2 - z(7)*z(14)*u2 - z(8)*(z(22)*u3-z(14)*u1);
cahat1p = z(23)*u1 + z(24)*u2 + z(25)*u3 - z(4)*z(15)*u1 - z(5)*z(15)*u2;
cahat2p = z(1)*z(15)*u1 + z(2)*z(15)*u2 - z(26)*u1 - z(27)*u2 - z(28)*u3;
cbhat1p = z(23)*u1 + z(24)*u2 + z(25)*u3 + z(22)*z(3)*u2 + z(3)*z(20)*u1 - z(1)*z(14)*u2 - z(1)*z(20)*u3 - z(4)*z(21)*u1 - z(6)*z(21)*u3 - z(2)*(z(22)*u3-z(14)*u1);
cbhat2p = z(22)*z(6)*u2 + z(1)*z(21)*u1 + z(3)*z(21)*u3 + z(6)*z(20)*u1 - z(26)*u1 - z(27)*u2 - z(28)*u3 - z(4)*z(14)*u2 - z(4)*z(20)*u3 - z(5)*(z(22)*u3-z(14)*u1);
% Update derivative array prior to integration step
VARp(1) = cahat1p;
VARp(2) = cahat2p;
VARp(3) = cahat3p;
VARp(4) = cbhat1p;
VARp(5) = cbhat2p;
VARp(6) = cbhat3p;
VARp(7) = e1p;
VARp(8) = e2p;
VARp(9) = e3p;
VARp(10) = e4p;
VARp(11) = posa1p;
VARp(12) = posa2p;
VARp(13) = posa3p;
VARp(14) = posb1p;
VARp(15) = posb2p;
VARp(16) = posb3p;
VARp(17) = u1p;
sys = VARp';
%===========================================================================
% mdlOutputs: Calculates and return the outputs
%===========================================================================
function Output = mdlOutputs(T,VAR,u)
global i11 i22 i33 l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p u1p a11 a12 a13 a14 a15 a21 a22 a23 a24 a25 a31 a32 a33 a34 a35 a41 a42 a43 a44 a45 a51 a52 a53 a54 a55;
global DEGtoRAD RADtoDEG z;
global TINITIAL TFINAL INTEGSTP PRINTINT ABSERR RELERR;
% Evaluate output quantities
ke = i11*u1^2 + 0.5*i22*u2^2 + 0.5*i22*u3^2 + 0.5*i33*u2^2 + 0.5*i33*u3^2 + 0.5*m*(z(22)^2*u2^2+(z(22)*u3-z(14)*u1)^2+2*z(22)*z(3)*u2*(z(23)*u1+z(24)*u2+z(25)*u3)+2*z(22)*z(5)*u3*(z(26)*u1+z(27)*u2+z(28)*u3)-2*z(22)*z(2)*u3*(z(23)*u1+z(24)*u2+z(25)*u3)-2*z(22)*z(6)*u2*(z(26)*u1+z(27)*u2+z(28)*u3)) + 0.5*m*(z(14)^2*u1^2+2*z(14)^2*u2^2+2*(z(23)*u1+z(24)*u2+z(25)*u3)^2+2*(z(26)*u1+z(27)*u2+z(28)*u3)^2+4*z(2)*z(14)*u1*(z(23)*u1+z(24)*u2+z(25)*u3)+4*z(4)*z(14)*u2*(z(26)*u1+z(27)*u2+z(28)*u3)-4*z(1)*z(14)*u2*(z(23)*u1+z(24)*u2+z(25)*u3)-4*z(5)*z(14)*u1*(z(26)*u1+z(27)*u2+z(28)*u3));
pe = 9.810000000000001*m*(posa3+posb3);
z(145) = z(50)*z(101) + z(51)*z(100) - i22*z(50)*u3 - i33*z(51)*u2;
z(146) = z(50)*z(106) + z(51)*z(105) - i22*z(51)*u2 - i33*z(50)*u3;
z(147) = z(14)*u1 + z(14)*z(50)*u2 + z(111) + z(50)*z(110);
z(148) = 0.5*e4 - 0.5*e2*z(51) - 0.5*e3*z(50);
z(149) = 0.5*e3 + 0.5*e1*z(51) + 0.5*e4*z(50);
z(150) = -4*e1*z(148) - 4*e2*z(149);
z(151) = z(150)*(z(11)+z(9)^2/z(10)^0.5)/z(11)^2;
z(152) = rad*(z(151)*u1+z(56));
z(153) = z(51)*z(110) + z(152) - z(14)*z(50)*u3;
z(154) = rad*z(9)*z(150)/(z(10)^0.5*z(11)^2);
z(155) = 0.5*e3*z(51) - 0.5*e1 - 0.5*e2*z(50);
z(156) = 0.5*e1*z(50) - 0.5*e2 - 0.5*e4*z(51);
z(157) = 2*e1*z(149) + 2*e2*z(148) + 2*e3*z(155) + 2*e4*z(156);
z(158) = z(4)*z(154) + z(15)*z(157);
z(159) = -4*e1*z(148) - 4*e3*z(156);
z(160) = z(5)*z(154) + z(15)*z(159);
z(161) = 2*e2*z(156) + 2*e3*z(149) - 2*e1*z(155) - 2*e4*z(148);
z(162) = z(6)*z(154) + z(15)*z(161);
z(163) = z(158)*u1 + z(160)*u2 + z(162)*u3 + z(81) + z(50)*z(71) - z(51)*z(76);
z(164) = rad*(z(151)*u2+z(50)*z(56));
z(165) = z(51)*z(111) - z(14)*u3 - z(164);
z(166) = -4*e2*z(149) - 4*e3*z(156);
z(167) = z(1)*z(154) + z(15)*z(166);
z(168) = 2*e1*z(149) + 2*e2*z(148) - 2*e3*z(155) - 2*e4*z(156);
z(169) = z(2)*z(154) + z(15)*z(168);
z(170) = 2*e1*z(156) + 2*e2*z(155) + 2*e3*z(148) + 2*e4*z(149);
z(171) = z(3)*z(154) + z(15)*z(170);
z(172) = z(51)*z(78) - z(167)*u1 - z(169)*u2 - z(171)*u3 - z(83) - z(50)*z(73);
z(173) = z(14)*z(50)*u2 + (z(14)+z(22)*z(51))*u1 + z(121) + z(50)*z(110);
z(174) = z(51)*z(110) + z(152) - z(22)*z(50)*u1 - z(14)*z(50)*u3 - z(122);
z(175) = z(22)*z(50)*u2 + z(50)*z(122) + z(51)*z(121) - (z(14)+z(22)*z(51))*u3 - z(164);
z(176) = z(145) + z(146) + m*(z(3)*z(23)*z(147)+z(14)*z(153)+2*z(23)*z(163)+z(1)*z(23)*z(165)+z(2)*z(23)*z(153)+2*z(5)*z(14)*z(172)-z(6)*z(26)*z(147)-z(4)*z(26)*z(165)-z(5)*z(26)*z(153)) + m*(z(3)*z(23)*z(173)+z(14)*z(174)+z(1)*z(23)*z(175)+z(2)*z(23)*z(174)+2*z(2)*z(14)*z(163)-z(6)*z(26)*z(173)-2*z(26)*z(172)-z(4)*z(26)*z(175)-z(5)*z(26)*z(174));
z(177) = i11*u3 + i33*z(51)*u1 - z(101) - z(51)*z(99);
z(178) = z(106) + z(51)*z(99) - i11*u3 - i22*z(51)*u1;
z(179) = z(177) + m*(z(3)*z(24)*z(147)+2*z(24)*z(163)+z(1)*z(24)*z(165)+z(2)*z(24)*z(153)-z(6)*z(27)*z(147)-2*z(27)*z(172)-z(14)*z(165)-2*z(1)*z(14)*z(163)-z(4)*z(27)*z(165)-z(5)*z(27)*z(153)) + m*(z(22)*z(173)+z(3)*z(24)*z(173)+z(22)*z(3)*z(163)+z(22)*z(6)*z(172)+z(1)*z(24)*z(175)+z(2)*z(24)*z(174)-z(6)*z(27)*z(173)-z(14)*z(175)-2*z(4)*z(14)*z(172)-z(4)*z(27)*z(175)-z(5)*z(27)*z(174)) - z(178);
z(180) = i22*z(50)*u1 + z(100) - i11*u2 - z(50)*z(99);
z(181) = i33*z(50)*u1 + z(105) - i11*u2 - z(50)*z(99);
z(182) = z(180) + z(181) + m*(z(3)*z(25)*z(147)+2*z(25)*z(163)+z(1)*z(25)*z(165)+z(2)*z(25)*z(153)-z(6)*z(28)*z(147)-2*z(28)*z(172)-z(4)*z(28)*z(165)-z(5)*z(28)*z(153)) + m*(z(3)*z(25)*z(173)+z(1)*z(25)*z(175)+z(2)*z(25)*z(174)-z(6)*z(28)*z(173)-z(22)*z(174)-z(22)*z(2)*z(163)-z(22)*z(5)*z(172)-z(4)*z(28)*z(175)-z(5)*z(28)*z(174));
z(183) = z(176) + z(50)*z(179) - z(51)*z(182);
z(184) = 2*e1*z(155) + 2*e2*z(156) + 2*e3*z(149) + 2*e4*z(148);
z(185) = z(184)*(z(17)+z(8)^2/z(16)^0.5)/z(17)^2;
z(186) = z(14)*z(184) + z(20)*z(150) + rad*z(8)*z(151) + rad*z(9)*z(185);
z(187) = 2*e1*z(156) + 2*e3*z(148) - 2*e2*z(155) - 2*e4*z(149);
z(188) = z(22)*z(150) - z(14)*z(187) - rad*z(7)*z(151);
z(189) = -z(22)*z(184) - z(20)*z(187) - rad*z(7)*z(185);
z(190) = z(186)*u1 + z(188)*u2 + z(189)*u3 + z(88) + z(50)*z(87) - z(51)*z(89);
z(191) = rad*z(8)*z(184)/(z(16)^0.5*z(17)^2);
z(192) = z(158) - z(4)*z(191) - z(21)*z(157);
z(193) = 2*rad*z(14)*z(151) + 2*rad*z(20)*z(185) + 2*z(20)*z(35)*z(184) + 2*rad*z(8)*z(35)*z(185) + 2*z(35)*(z(154)-z(191)) + 2*z(8)*z(20)*(z(154)-z(191)) - 2*z(22)*z(35)*z(187) - 2*z(14)*z(35)*z(150) - 2*rad*z(9)*z(35)*z(151) - 2*z(22)*z(7)*(z(154)-z(191)) - 2*z(9)*z(14)*(z(154)-z(191));
z(194) = z(22)*z(193)/(z(37)^0.5*z(38)^2);
z(195) = (2*rad*z(38)*z(151)-z(14)*z(193)/z(37)^0.5)/z(38)^2;
z(196) = (2*rad*z(38)*z(185)-z(20)*z(193)/z(37)^0.5)/z(38)^2;
z(197) = z(1)*z(191) + z(21)*z(166) - z(167);
z(198) = z(14)*z(159)*(z(1)*z(41)+z(3)*z(39)) + rad*z(5)*z(151)*(z(1)*z(41)+z(3)*z(39)) + z(20)*z(161)*(z(1)*z(41)-z(2)*z(42)) + rad*z(6)*z(185)*(z(1)*z(41)-z(2)*z(42)) + z(192)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(2)*z(14)*(z(4)*z(194)-2*z(39)*z(161)-2*z(41)*z(157)-z(6)*z(195)) + 0.5*z(6)*z(20)*(2*z(41)*z(166)-2*z(42)*z(168)-z(1)*z(194)-z(2)*z(196)) + 0.5*z(32)*(z(3)*z(195)+2*z(39)*z(170)+2*z(41)*z(166)-2*z(42)*z(168)-z(1)*z(194)-z(2)*z(196)) - z(14)*z(168)*(z(4)*z(41)+z(6)*z(39)) - rad*z(2)*z(151)*(z(4)*z(41)+z(6)*z(39)) - z(20)*z(170)*(z(4)*z(41)-z(5)*z(42)) - rad*z(3)*z(185)*(z(4)*z(41)-z(5)*z(42)) - z(197)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(5)*z(14)*(z(1)*z(194)-2*z(39)*z(170)-2*z(41)*z(166)-z(3)*z(195)) - 0.5*z(3)*z(20)*(2*z(41)*z(157)-2*z(42)*z(159)-z(4)*z(194)-z(5)*z(196)) - 0.5*z(29)*(z(6)*z(195)+2*z(39)*z(161)+2*z(41)*z(157)-2*z(42)*z(159)-z(4)*z(194)-z(5)*z(196));
z(199) = z(160) - z(5)*z(191) - z(21)*z(159);
z(200) = z(2)*z(191) + z(21)*z(168) - z(169);
z(201) = z(22)*z(161)*(z(1)*z(41)-z(2)*z(42)) + z(14)*z(157)*(z(2)*z(42)-z(3)*z(39)) + rad*z(4)*z(151)*(z(2)*z(42)-z(3)*z(39)) + z(199)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(1)*z(14)*(z(6)*z(195)+2*z(39)*z(161)-2*z(42)*z(159)-z(5)*z(196)) + 0.5*z(22)*z(6)*(2*z(41)*z(166)-2*z(42)*z(168)-z(1)*z(194)-z(2)*z(196)) + 0.5*z(33)*(z(3)*z(195)+2*z(39)*z(170)+2*z(41)*z(166)-2*z(42)*z(168)-z(1)*z(194)-z(2)*z(196)) - z(22)*z(170)*(z(4)*z(41)-z(5)*z(42)) - z(14)*z(166)*(z(5)*z(42)-z(6)*z(39)) - rad*z(1)*z(151)*(z(5)*z(42)-z(6)*z(39)) - z(200)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(4)*z(14)*(z(3)*z(195)+2*z(39)*z(170)-2*z(42)*z(168)-z(2)*z(196)) - 0.5*z(22)*z(3)*(2*z(41)*z(157)-2*z(42)*z(159)-z(4)*z(194)-z(5)*z(196)) - 0.5*z(30)*(z(6)*z(195)+2*z(39)*z(161)+2*z(41)*z(157)-2*z(42)*z(159)-z(4)*z(194)-z(5)*z(196));
z(202) = z(162) - z(6)*z(191) - z(21)*z(161);
z(203) = z(3)*z(191) + z(21)*z(170) - z(171);
z(204) = z(22)*z(168)*(z(4)*z(41)+z(6)*z(39)) + z(20)*z(157)*(z(2)*z(42)-z(3)*z(39)) + rad*z(4)*z(185)*(z(2)*z(42)-z(3)*z(39)) + z(202)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(22)*z(5)*(z(1)*z(194)-2*z(39)*z(170)-2*z(41)*z(166)-z(3)*z(195)) + 0.5*z(1)*z(20)*(z(6)*z(195)+2*z(39)*z(161)-2*z(42)*z(159)-z(5)*z(196)) + 0.5*z(34)*(z(3)*z(195)+2*z(39)*z(170)+2*z(41)*z(166)-2*z(42)*z(168)-z(1)*z(194)-z(2)*z(196)) - z(22)*z(159)*(z(1)*z(41)+z(3)*z(39)) - z(20)*z(166)*(z(5)*z(42)-z(6)*z(39)) - rad*z(1)*z(185)*(z(5)*z(42)-z(6)*z(39)) - z(203)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(2)*(z(4)*z(194)-2*z(39)*z(161)-2*z(41)*z(157)-z(6)*z(195)) - 0.5*z(4)*z(20)*(z(3)*z(195)+2*z(39)*z(170)-2*z(42)*z(168)-z(2)*z(196)) - 0.5*z(31)*(z(6)*z(195)+2*z(39)*z(161)+2*z(41)*z(157)-2*z(42)*z(159)-z(4)*z(194)-z(5)*z(196));
z(205) = z(198)*u1 + z(201)*u2 + z(204)*u3 + z(85) + z(50)*z(75) - z(51)*z(80);
z(206) = (z(44)*z(190)-z(48)*z(205))/z(49);
z(207) = (z(43)*z(190)-z(46)*z(205))/z(49);
z(208) = z(183) + z(138)*z(206) - z(139)*z(207);
a11 = -z(208)/z(141);
z(209) = e1*(z(11)+z(9)^2/z(10)^0.5)/z(11)^2;
z(210) = e4*(z(17)+z(8)^2/z(16)^0.5)/z(17)^2;
z(211) = e1*z(9)/(z(10)^0.5*z(11)^2);
z(212) = e4*z(8)/(z(16)^0.5*z(17)^2);
z(214) = 4*rad*z(20)*z(210) + 4*e4*z(20)*z(35) + 8*e1*z(14)*z(35) + 4*rad*z(8)*z(35)*z(210) + 8*rad*z(9)*z(35)*z(209) + 4*z(213)*z(7)*(z(212)+2*z(211)) + 4*rad*z(9)*z(14)*(z(212)+2*z(211)) - 8*rad*z(14)*z(209) - 4*z(22)*e3*z(35) - 4*rad*z(35)*(z(212)+2*z(211)) - 4*rad*z(8)*z(20)*(z(212)+2*z(211));
z(215) = z(22)*z(214)/(z(37)^0.5*z(38)^2);
z(216) = (8*rad*z(38)*z(209)+z(14)*z(214)/z(37)^0.5)/z(38)^2;
z(217) = (4*rad*z(38)*z(210)-z(20)*z(214)/z(37)^0.5)/z(38)^2;
z(218) = -2*e4*z(15) - 4*rad*z(6)*z(211);
z(219) = z(218) + 2*e4*z(21) - 2*rad*z(6)*z(212);
z(220) = 2*e3*z(15) - 4*rad*z(3)*z(211);
z(221) = 2*e3*z(21) + 2*rad*z(3)*z(212) - z(220);
z(222) = 2*z(22)*e2*(z(4)*z(41)+z(6)*z(39)) + 4*z(22)*e1*(z(1)*z(41)+z(3)*z(39)) + 2*e2*z(20)*(z(2)*z(42)-z(3)*z(39)) + 2*rad*z(4)*z(210)*(z(2)*z(42)-z(3)*z(39)) + z(219)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(22)*z(2)*(4*e2*z(41)-4*e4*z(39)-z(4)*z(215)-z(6)*z(216)) + 0.5*z(1)*z(20)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216)) + 0.5*z(34)*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216)) - 2*rad*z(1)*z(210)*(z(5)*z(42)-z(6)*z(39)) - z(221)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(5)*(4*e3*z(39)-z(1)*z(215)-z(3)*z(216)) - 0.5*z(4)*z(20)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216)) - 0.5*z(31)*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216));
z(223) = 2*e4*z(14) + 2*rad*z(9)*z(210) - 4*e1*z(20) - 4*rad*z(8)*z(209);
z(224) = 2*e2*z(15) - 4*rad*z(4)*z(211);
z(225) = z(224) - 2*e2*z(21) - 2*rad*z(4)*z(212);
z(226) = rad*z(1)*(z(212)+2*z(211));
z(227) = 4*rad*z(2)*z(209)*(z(4)*z(41)+z(6)*z(39)) + 2*rad*z(6)*z(210)*(z(1)*z(41)-z(2)*z(42)) + z(225)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(5)*z(14)*(4*e3*z(39)-z(1)*z(215)-z(3)*z(216)) + 0.5*z(32)*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216)) - 4*e1*z(14)*(z(1)*z(41)+z(3)*z(39)) - 2*e2*z(14)*(z(4)*z(41)+z(6)*z(39)) - 4*rad*z(5)*z(209)*(z(1)*z(41)+z(3)*z(39)) - 2*e3*z(20)*(z(4)*z(41)-z(5)*z(42)) - 2*e4*z(20)*(z(1)*z(41)-z(2)*z(42)) - 2*rad*z(3)*z(210)*(z(4)*z(41)-z(5)*z(42)) - 2*z(226)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(6)*z(20)*(z(1)*z(215)+z(2)*z(217)+4*e2*z(42)) - 0.5*z(2)*z(14)*(4*e2*z(41)-4*e4*z(39)-z(4)*z(215)-z(6)*z(216)) - 0.5*z(3)*z(20)*(4*e2*z(41)+8*e1*z(42)-z(4)*z(215)-z(5)*z(217)) - 0.5*z(29)*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216));
z(228) = -2*z(22)*e4 - 2*e3*z(20) - 2*rad*z(7)*z(210);
z(229) = -4*e1*z(15) - 4*rad*z(5)*z(211);
z(230) = z(229) + 4*e1*z(21) - 2*rad*z(5)*z(212);
z(231) = 2*e2*z(15) - 4*rad*z(2)*z(211);
z(232) = 2*e2*z(21) + 2*rad*z(2)*z(212) - z(231);
z(233) = 2*e2*z(14)*(z(2)*z(42)-z(3)*z(39)) + 4*rad*z(1)*z(209)*(z(5)*z(42)-z(6)*z(39)) + z(230)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(1)*z(14)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216)) + 0.5*z(33)*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216)) - 2*z(22)*e3*(z(4)*z(41)-z(5)*z(42)) - 2*z(22)*e4*(z(1)*z(41)-z(2)*z(42)) - 4*rad*z(4)*z(209)*(z(2)*z(42)-z(3)*z(39)) - z(232)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(6)*(z(1)*z(215)+z(2)*z(217)+4*e2*z(42)) - 0.5*z(22)*z(3)*(4*e2*z(41)+8*e1*z(42)-z(4)*z(215)-z(5)*z(217)) - 0.5*z(4)*z(14)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216)) - 0.5*z(30)*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216));
z(234) = 4*rad*z(7)*z(209) - 4*z(22)*e1 - 2*e3*z(14);
z(235) = z(43)*z(228) + z(48)*z(233) - z(44)*z(234) - z(46)*z(222);
z(236) = (z(235)*(z(44)*z(47)-z(45)*z(48))-z(49)*(z(44)*z(223)+z(47)*z(222)-z(45)*z(228)-z(48)*z(227)))/z(49)^2;
z(237) = (z(235)*(z(43)*z(47)-z(45)*z(46))+z(49)*(z(45)*z(234)+z(46)*z(227)-z(43)*z(223)-z(47)*z(233)))/z(49)^2;
z(238) = u1*(i22*z(236)*u3+i33*z(237)*u2-z(236)*z(101)-z(237)*z(100));
z(239) = u1*(i22*z(237)*u2+i33*z(236)*u3-z(236)*z(106)-z(237)*z(105));
z(240) = -4*rad*z(209)*u2 - z(14)*z(236)*u1;
z(241) = u2*z(240) - 4*rad*z(209)*u1^2 - z(236)*u1*z(110);
z(242) = (e2*z(237)+e3*z(236))*u1;
z(243) = -0.5*u3 - 0.5*e1*z(237)*u1 - 0.5*e4*z(236)*u1;
z(244) = -4*e2*z(243) - 2*e1*z(242) - 4*e1p;
z(245) = (z(9)*z(10)*z(244)-4*e1*z(10)*z(54)-4*e1*z(9)^2*z(54))/z(10)^1.5;
z(246) = (4*e1*z(9)*(z(11)*z(54)+2*z(9)*z(55))/z(10)^0.5+z(11)*(4*e1*z(55)-z(9)*z(245)-z(11)*z(244)))/z(11)^3;
z(247) = -u3*z(240) - z(237)*u1*z(110) - rad*u1*z(246);
z(248) = rad*(z(11)*z(245)-8*e1*z(9)*z(55)/z(10)^0.5)/z(11)^3;
z(249) = (-1+e2*z(236)-e3*z(237))*u1;
z(250) = 0.5*u2 + 0.5*e4*z(237)*u1 - 0.5*e1*z(236)*u1;
z(251) = e2*z(242) + e3*z(249) + 2*e1*z(243) + 2*e4*z(250) + 2*e2p;
z(252) = z(4)*z(248) + z(15)*z(251) + 2*e2*z(60) - 4*rad*z(211)*z(67);
z(253) = -4*e3*z(250) - 2*e1*z(242) - 4*e1p;
z(254) = z(5)*z(248) + z(15)*z(253) - 4*e1*z(60) - 4*rad*z(211)*z(70);
z(255) = 2*e2*z(250) + 2*e3*z(243) - e1*z(249) - e4*z(242) - 2*e4p;
z(256) = z(6)*z(248) + z(15)*z(255) - 2*e4*z(60) - 4*rad*z(211)*z(52);
z(257) = u1*z(252) + u2*z(254) + u3*z(256) + z(237)*u1*z(76) - z(236)*u1*z(71);
z(258) = rad*(u2*z(246)+z(236)*u1*z(56));
z(259) = 4*rad*z(209)*u1*u3 + z(258) - z(237)*u1*z(111);
z(260) = -4*e2*z(243) - 4*e3*z(250);
z(261) = z(15)*z(260) + z(1)*z(248) - 4*rad*z(211)*z(53);
z(262) = e2*z(242) + 2*e1*z(243) + 2*e2p - 2*e4*z(250) - e3*z(249);
z(263) = z(2)*z(248) + z(15)*z(262) + 2*e2*z(60) - 4*rad*z(211)*z(65);
z(264) = e2*z(249) + e3*z(242) + 2*e1*z(250) + 2*e4*z(243) + 2*e3p;
z(265) = z(3)*z(248) + z(15)*z(264) + 2*e3*z(60) - 4*rad*z(211)*z(68);
z(266) = z(236)*u1*z(73) - u1*z(261) - u2*z(263) - u3*z(265) - z(237)*u1*z(78);
z(267) = (z(22)*z(237)+4*rad*z(209))*u1;
z(268) = u2*z(240) - u1*z(267) - z(236)*u1*z(110);
z(269) = z(22)*z(236)*u1^2 - u3*z(240) - z(237)*u1*z(110) - rad*u1*z(246);
z(270) = u3*z(267) + z(258) - z(22)*z(236)*u1*u2 - z(236)*u1*z(122) - z(237)*u1*z(121);
z(271) = z(238) + z(239) - m*(z(6)*z(26)*z(241)+z(4)*z(26)*z(259)+z(5)*z(26)*z(247)+2*e2*z(26)*z(118)+4*rad*z(209)*z(120)+8*e1*z(14)*z(117)+8*rad*z(5)*z(209)*z(117)-2*e3*z(23)*z(119)-2*e4*z(26)*z(119)-z(3)*z(23)*z(241)-z(3)*z(224)*z(119)-4*rad*z(1)*z(6)*z(211)*z(119)-2*z(23)*z(257)-2*z(224)*z(116)-z(14)*z(247)-4*e1*z(26)*z(120)-2*e2*z(23)*z(120)-2*z(5)*z(14)*z(266)-z(1)*z(23)*z(259)-z(1)*z(224)*z(118)-z(2)*z(23)*z(247)-z(2)*z(224)*z(120)-4*rad*z(1)*z(4)*z(211)*z(118)-4*rad*z(1)*z(5)*z(211)*z(120)) - m*(z(6)*z(26)*z(268)+2*z(26)*z(266)+z(4)*z(26)*z(270)+z(5)*z(26)*z(269)+2*e2*z(26)*z(123)+4*rad*z(209)*z(125)+8*rad*z(2)*z(209)*z(116)-2*e3*z(23)*z(124)-2*e4*z(26)*z(124)-z(3)*z(23)*z(268)-z(3)*z(224)*z(124)-4*rad*z(1)*z(6)*z(211)*z(124)-z(14)*z(269)-4*e1*z(26)*z(125)-4*e2*z(14)*z(116)-2*e2*z(23)*z(125)-2*z(2)*z(14)*z(257)-z(1)*z(23)*z(270)-z(1)*z(224)*z(123)-z(2)*z(23)*z(269)-z(2)*z(224)*z(125)-8*rad*z(1)*z(211)*z(117)-4*rad*z(1)*z(4)*z(211)*z(123)-4*rad*z(1)*z(5)*z(211)*z(125));
z(272) = z(237)*u1*(i33*u1-z(99));
z(273) = z(237)*u1*(i22*u1-z(99));
z(274) = -z(272) - z(273) - m*(z(6)*z(27)*z(241)+z(6)*z(231)*z(119)+z(14)*z(259)+2*z(27)*z(266)+2*z(231)*z(117)+z(4)*z(27)*z(259)+z(4)*z(231)*z(118)+z(5)*z(27)*z(247)+z(5)*z(231)*z(120)+2*e2*z(27)*z(118)+2*z(1)*z(14)*z(257)-2*e3*z(24)*z(119)-2*e4*z(27)*z(119)-z(3)*z(24)*z(241)-z(3)*z(229)*z(119)-2*z(24)*z(257)-2*z(229)*z(116)-4*rad*z(209)*z(118)-4*e1*z(27)*z(120)-2*e2*z(24)*z(120)-z(1)*z(24)*z(259)-z(1)*z(229)*z(118)-z(2)*z(24)*z(247)-z(2)*z(229)*z(120)-8*rad*z(1)*z(209)*z(116)) - m*(z(6)*z(27)*z(268)+z(6)*z(231)*z(124)+z(14)*z(270)+z(4)*z(27)*z(270)+z(4)*z(231)*z(123)+z(5)*z(27)*z(269)+z(5)*z(231)*z(125)+2*z(22)*e4*z(117)+2*e2*z(27)*z(123)+2*z(4)*z(14)*z(266)+4*e2*z(14)*z(117)-z(22)*z(268)-2*e3*z(24)*z(124)-2*e4*z(27)*z(124)-z(3)*z(24)*z(268)-z(3)*z(229)*z(124)-4*rad*z(209)*z(123)-4*e1*z(27)*z(125)-2*z(22)*e3*z(116)-2*e2*z(24)*z(125)-z(22)*z(3)*z(257)-z(22)*z(6)*z(266)-z(1)*z(24)*z(270)-z(1)*z(229)*z(123)-z(2)*z(24)*z(269)-z(2)*z(229)*z(125)-8*rad*z(4)*z(209)*z(117));
z(275) = z(236)*u1*(i22*u1-z(99));
z(276) = z(236)*u1*(i33*u1-z(99));
z(277) = m*(z(3)*z(25)*z(241)+z(3)*z(218)*z(119)+2*e3*z(25)*z(119)+2*e4*z(28)*z(119)+2*z(25)*z(257)+2*z(218)*z(116)+z(1)*z(25)*z(259)+z(1)*z(218)*z(118)+z(2)*z(25)*z(247)+z(2)*z(218)*z(120)+2*e2*z(25)*z(120)+4*e1*z(28)*z(120)-z(6)*z(28)*z(241)-z(6)*z(220)*z(119)-2*z(28)*z(266)-2*z(220)*z(117)-2*e2*z(28)*z(118)-z(4)*z(28)*z(259)-z(4)*z(220)*z(118)-z(5)*z(28)*z(247)-z(5)*z(220)*z(120)) + m*(z(3)*z(25)*z(268)+z(3)*z(218)*z(124)+2*e3*z(25)*z(124)+2*e4*z(28)*z(124)+z(1)*z(25)*z(270)+z(1)*z(218)*z(123)+z(2)*z(25)*z(269)+z(2)*z(218)*z(125)+2*e2*z(25)*z(125)+4*z(22)*e1*z(117)+4*e1*z(28)*z(125)-z(6)*z(28)*z(268)-z(6)*z(220)*z(124)-z(22)*z(269)-2*z(22)*e2*z(116)-2*e2*z(28)*z(123)-z(22)*z(2)*z(257)-z(22)*z(5)*z(266)-z(4)*z(28)*z(270)-z(4)*z(220)*z(123)-z(5)*z(28)*z(269)-z(5)*z(220)*z(125)) - z(275) - z(276);
z(278) = z(271) + z(50)*z(274) + z(237)*z(136) - z(51)*z(277) - z(236)*z(133);
z(279) = m*(2*z(23)*z(229)+2*z(24)*z(224)+2*z(26)*z(231)+z(22)*z(3)*z(224)+2*z(22)*e3*z(23)+2*z(22)*e4*z(26)+2*z(2)*z(14)*z(229)+4*e2*z(14)*z(24)+4*e2*z(14)*z(26)+8*e1*z(14)*z(27)+4*z(213)*z(1)*z(6)*z(211)+8*rad*z(1)*z(23)*z(209)+8*rad*z(5)*z(27)*z(209)-2*z(1)*z(14)*z(224)-2*z(5)*z(14)*z(231)-8*rad*z(1)*z(27)*z(211)-8*rad*z(2)*z(24)*z(209)-8*rad*z(4)*z(26)*z(209)-8*rad*z(1)*z(4)*z(14)*z(211));
z(280) = m*(2*z(24)*z(229)+2*z(27)*z(231)+z(22)*z(3)*z(229)+2*z(22)*e3*z(24)+2*z(22)*e4*z(27)+2*z(4)*z(14)*z(231)+4*e2*z(14)*z(27)+8*rad*z(1)*z(24)*z(209)-8*rad*z(14)*z(209)-2*z(1)*z(14)*z(229)-z(22)*z(6)*z(231)-8*rad*z(4)*z(27)*z(209));
z(281) = m*(2*z(24)*z(218)+2*z(25)*z(229)+2*z(27)*z(220)+2*z(28)*z(231)+z(22)*z(3)*z(218)+z(22)*z(5)*z(231)+2*z(22)*e3*z(25)+2*z(22)*e4*z(28)+2*z(4)*z(14)*z(220)+4*e2*z(14)*z(28)+8*rad*z(1)*z(25)*z(209)-4*z(22)*e1*z(27)-2*z(22)*e2*z(24)-2*z(1)*z(14)*z(218)-z(22)*z(2)*z(229)-z(22)*z(6)*z(220)-8*rad*z(4)*z(28)*z(209));
z(282) = z(279) + 2*z(50)*z(280) - z(51)*z(281) - z(131)*z(236) - z(132)*z(237);
z(283) = e1*z(249) + e4*z(242) + 2*e2*z(250) + 2*e3*z(243) + 2*e4p;
z(284) = (z(8)*z(16)*z(283)+2*e4*z(16)*z(57)+2*e4*z(8)^2*z(57))/z(16)^1.5;
z(285) = (z(17)*(z(8)*z(284)+z(17)*z(283)+2*e4*z(58))+2*e4*z(8)*(z(17)*z(57)+2*z(8)*z(58))/z(16)^0.5)/z(17)^3;
z(286) = z(14)*z(283) + z(20)*z(244) + rad*z(9)*z(285) + 2*rad*e4*z(56) + 2*rad*z(210)*z(54) - 4*rad*e1*z(59) - 4*rad*z(209)*z(57) - rad*z(8)*z(246);
z(287) = e3*z(242) + 2*e1*z(250) + 2*e3p - 2*e4*z(243) - e2*z(249);
z(288) = z(22)*z(244) + rad*z(7)*z(246) + 4*rad*z(209)*z(62) - z(14)*z(287) - 2*rad*e3*z(56);
z(289) = -z(22)*z(283) - z(20)*z(287) - 2*rad*e3*z(59) - 2*rad*z(210)*z(62) - rad*z(7)*z(285);
z(290) = u1*z(286) + u2*z(288) + u3*z(289) + z(237)*u1*z(89) - z(236)*u1*z(87);
z(291) = rad*(z(17)*z(284)+4*e4*z(8)*z(58)/z(16)^0.5)/z(17)^3;
z(292) = z(252) - 2*e2*z(61) - z(4)*z(291) - z(21)*z(251) - 2*rad*z(212)*z(67);
z(294) = 2*rad*z(20)*z(285) + 2*z(20)*z(35)*z(283) + 4*z(293)*z(210)*z(59) + 2*rad*z(8)*z(35)*z(285) + 2*rad*z(9)*z(35)*z(246) + 4*rad*e4*z(35)*z(59) + 4*rad*z(35)*z(210)*z(57) + 8*rad*e1*z(35)*z(56) + 8*rad*z(35)*z(209)*z(54) + 2*z(35)*(z(248)-z(291)) + 2*z(8)*z(20)*(z(248)-z(291)) + 4*z(213)*(z(212)+2*z(211))*z(62) + 4*e4*z(20)*(z(60)-z(61)) + 8*e1*z(14)*(z(60)-z(61)) + 4*rad*z(8)*z(210)*(z(60)-z(61)) + 4*rad*z(14)*(z(212)+2*z(211))*z(54) + 4*z(293)*z(9)*(z(212)+2*z(211))*z(56) + 8*rad*z(9)*z(209)*(z(60)-z(61)) - 8*z(293)*z(209)*z(56) - 2*rad*z(14)*z(246) - 2*z(22)*z(35)*z(287) - 2*z(14)*z(35)*z(244) - 4*z(22)*e3*(z(60)-z(61)) - 2*z(22)*z(7)*(z(248)-z(291)) - 2*z(9)*z(14)*(z(248)-z(291)) - 4*rad*z(20)*(z(212)+2*z(211))*z(57) - 4*z(293)*z(8)*(z(212)+2*z(211))*z(59) - 4*rad*(z(212)+2*z(211))*(z(60)-z(61));
z(295) = z(22)*(2*z(294)/(z(37)^0.5*z(38)^2)-z(214)*(2+z(38)/z(37)^0.5)*z(63)/(z(37)*z(38)^3));
z(296) = (2*z(214)*(rad*z(38)*z(56)-z(14)*z(63)/z(37)^0.5)/z(37)^0.5+z(38)*(4*rad*z(38)*z(246)+(2*z(14)*z(37)*z(294)-z(14)*z(214)*z(63)-8*rad*z(37)*z(209)*z(63))/z(37)^1.5))/z(38)^3;
z(297) = (2*z(20)*z(214)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(285)+(z(20)*z(214)*z(63)-2*z(20)*z(37)*z(294)-4*rad*z(37)*z(210)*z(63)-2*rad*z(37)*z(214)*z(59))/z(37)^1.5))/z(38)^3;
z(298) = z(21)*z(260) + z(1)*z(291) + 2*rad*z(212)*z(53) - z(261);
z(299) = z(14)*(z(1)*z(41)+z(3)*z(39))*z(253) + rad*z(2)*(z(4)*z(41)+z(6)*z(39))*z(246) + z(20)*(z(1)*z(41)-z(2)*z(42))*z(255) + 4*rad*z(209)*(z(4)*z(41)+z(6)*z(39))*z(65) + rad*z(6)*(z(1)*z(41)-z(2)*z(42))*z(285) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(292) + 2*rad*z(210)*(z(1)*z(41)-z(2)*z(42))*z(52) + 0.5*z(14)*(4*e3*z(39)-z(1)*z(215)-z(3)*z(216))*z(70) + 0.5*rad*z(5)*(4*e3*z(39)-z(1)*z(215)-z(3)*z(216))*z(56) + e2*z(14)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) + 2*e1*z(14)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) + 2*rad*z(5)*z(209)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) + rad*z(6)*z(210)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216))*z(84) + z(226)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.25*z(5)*z(14)*(4*z(41)*z(260)+4*e3*z(69)+4*z(39)*z(264)-2*z(215)*z(53)-2*z(216)*z(68)-z(1)*z(295)-z(3)*z(296)) + 0.25*z(6)*z(20)*(4*z(41)*z(260)-4*e2*z(66)-4*z(42)*z(262)-2*z(215)*z(53)-2*z(217)*z(65)-z(1)*z(295)-z(2)*z(297)) + 0.25*z(32)*(4*z(41)*z(260)+4*e3*z(69)+4*z(39)*z(264)-4*e2*z(66)-4*z(42)*z(262)-2*z(215)*z(53)-2*z(216)*z(68)-2*z(217)*z(65)-z(1)*z(295)-z(2)*z(297)-z(3)*z(296)) - z(14)*(z(4)*z(41)+z(6)*z(39))*z(262) - 4*rad*e1*(z(1)*z(41)+z(3)*z(39))*z(56) - 4*rad*z(209)*(z(1)*z(41)+z(3)*z(39))*z(70) - 2*rad*e2*(z(4)*z(41)+z(6)*z(39))*z(56) - rad*z(5)*(z(1)*z(41)+z(3)*z(39))*z(246) - z(20)*(z(4)*z(41)-z(5)*z(42))*z(264) - 2*rad*e3*(z(4)*z(41)-z(5)*z(42))*z(59) - 2*rad*e4*(z(1)*z(41)-z(2)*z(42))*z(59) - 2*rad*z(210)*(z(4)*z(41)-z(5)*z(42))*z(68) - rad*z(3)*(z(4)*z(41)-z(5)*z(42))*z(285) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(298) - 0.5*z(20)*(z(1)*z(215)+z(2)*z(217)+4*e2*z(42))*z(52) - 0.5*rad*z(6)*(z(1)*z(215)+z(2)*z(217)+4*e2*z(42))*z(59) - 2*rad*z(2)*z(209)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e3*z(20)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - e4*z(20)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - 0.5*z(14)*(4*e2*z(41)-4*e4*z(39)-z(4)*z(215)-z(6)*z(216))*z(65) - 0.5*z(20)*(4*e2*z(41)+8*e1*z(42)-z(4)*z(215)-z(5)*z(217))*z(68) - rad*z(3)*z(210)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*rad*z(2)*(4*e2*z(41)-4*e4*z(39)-z(4)*z(215)-z(6)*z(216))*z(56) - 0.5*rad*z(3)*(4*e2*z(41)+8*e1*z(42)-z(4)*z(215)-z(5)*z(217))*z(59) - 0.5*z(225)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216))*z(82) - 0.25*z(2)*z(14)*(4*z(39)*z(255)+4*z(41)*z(251)-4*e2*z(64)-4*e4*z(69)-2*z(215)*z(67)-2*z(216)*z(52)-z(4)*z(295)-z(6)*z(296)) - 0.25*z(3)*z(20)*(4*z(41)*z(251)+8*e1*z(66)-4*e2*z(64)-4*z(42)*z(253)-2*z(215)*z(67)-2*z(217)*z(70)-z(4)*z(295)-z(5)*z(297)) - 0.25*z(29)*(4*z(39)*z(255)+4*z(41)*z(251)+8*e1*z(66)-4*e2*z(64)-4*e4*z(69)-4*z(42)*z(253)-2*z(215)*z(67)-2*z(216)*z(52)-2*z(217)*z(70)-z(4)*z(295)-z(5)*z(297)-z(6)*z(296));
z(300) = z(254) + 4*e1*z(61) - z(5)*z(291) - z(21)*z(253) - 2*rad*z(212)*z(70);
z(301) = z(2)*z(291) + z(21)*z(262) + 2*e2*z(61) + 2*rad*z(212)*z(65) - z(263);
z(302) = z(22)*(z(1)*z(41)-z(2)*z(42))*z(255) + z(14)*(z(2)*z(42)-z(3)*z(39))*z(251) + rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(246) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(300) + 2*rad*e2*(z(2)*z(42)-z(3)*z(39))*z(56) + 4*rad*z(209)*(z(5)*z(42)-z(6)*z(39))*z(53) + 0.5*z(14)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216))*z(53) + 2*rad*z(4)*z(209)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) + 0.5*rad*z(1)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216))*z(56) + 0.5*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216))*z(74) + 0.5*z(232)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.25*z(22)*z(6)*(4*z(41)*z(260)-4*e2*z(66)-4*z(42)*z(262)-2*z(215)*z(53)-2*z(217)*z(65)-z(1)*z(295)-z(2)*z(297)) + 0.25*z(1)*z(14)*(4*z(39)*z(255)+8*e1*z(66)-4*e4*z(69)-4*z(42)*z(253)-2*z(216)*z(52)-2*z(217)*z(70)-z(5)*z(297)-z(6)*z(296)) + 0.25*z(33)*(4*z(41)*z(260)+4*e3*z(69)+4*z(39)*z(264)-4*e2*z(66)-4*z(42)*z(262)-2*z(215)*z(53)-2*z(216)*z(68)-2*z(217)*z(65)-z(1)*z(295)-z(2)*z(297)-z(3)*z(296)) - z(14)*(z(5)*z(42)-z(6)*z(39))*z(260) - z(22)*(z(4)*z(41)-z(5)*z(42))*z(264) - 4*rad*z(209)*(z(2)*z(42)-z(3)*z(39))*z(67) - rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(246) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(301) - 0.5*z(22)*(z(1)*z(215)+z(2)*z(217)+4*e2*z(42))*z(52) - e2*z(14)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 2*rad*z(1)*z(209)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - z(22)*e3*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - z(22)*e4*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - 0.5*z(22)*(4*e2*z(41)+8*e1*z(42)-z(4)*z(215)-z(5)*z(217))*z(68) - 0.5*z(14)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216))*z(67) - 0.5*rad*z(4)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216))*z(56) - 0.5*z(230)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216))*z(72) - 0.25*z(22)*z(3)*(4*z(41)*z(251)+8*e1*z(66)-4*e2*z(64)-4*z(42)*z(253)-2*z(215)*z(67)-2*z(217)*z(70)-z(4)*z(295)-z(5)*z(297)) - 0.25*z(4)*z(14)*(4*e3*z(69)+4*z(39)*z(264)-4*e2*z(66)-4*z(42)*z(262)-2*z(216)*z(68)-2*z(217)*z(65)-z(2)*z(297)-z(3)*z(296)) - 0.25*z(30)*(4*z(39)*z(255)+4*z(41)*z(251)+8*e1*z(66)-4*e2*z(64)-4*e4*z(69)-4*z(42)*z(253)-2*z(215)*z(67)-2*z(216)*z(52)-2*z(217)*z(70)-z(4)*z(295)-z(5)*z(297)-z(6)*z(296));
z(303) = z(256) + 2*e4*z(61) - z(6)*z(291) - z(21)*z(255) - 2*rad*z(212)*z(52);
z(304) = z(3)*z(291) + z(21)*z(264) + 2*e3*z(61) + 2*rad*z(212)*z(68) - z(265);
z(305) = z(22)*(z(4)*z(41)+z(6)*z(39))*z(262) + z(20)*(z(2)*z(42)-z(3)*z(39))*z(251) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(285) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(303) + 2*rad*e2*(z(2)*z(42)-z(3)*z(39))*z(59) + 2*rad*z(210)*(z(2)*z(42)-z(3)*z(39))*z(67) + rad*z(1)*z(210)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(22)*(4*e2*z(41)-4*e4*z(39)-z(4)*z(215)-z(6)*z(216))*z(65) + 0.5*z(20)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216))*z(53) + 0.5*rad*z(1)*(8*e1*z(42)-4*e4*z(39)-z(5)*z(217)-z(6)*z(216))*z(59) + 0.5*(4*e3*z(39)-4*e2*z(42)-z(1)*z(215)-z(2)*z(217)-z(3)*z(216))*z(79) + 0.5*z(221)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.25*z(22)*z(2)*(4*z(39)*z(255)+4*z(41)*z(251)-4*e2*z(64)-4*e4*z(69)-2*z(215)*z(67)-2*z(216)*z(52)-z(4)*z(295)-z(6)*z(296)) + 0.25*z(1)*z(20)*(4*z(39)*z(255)+8*e1*z(66)-4*e4*z(69)-4*z(42)*z(253)-2*z(216)*z(52)-2*z(217)*z(70)-z(5)*z(297)-z(6)*z(296)) + 0.25*z(34)*(4*z(41)*z(260)+4*e3*z(69)+4*z(39)*z(264)-4*e2*z(66)-4*z(42)*z(262)-2*z(215)*z(53)-2*z(216)*z(68)-2*z(217)*z(65)-z(1)*z(295)-z(2)*z(297)-z(3)*z(296)) - z(20)*(z(5)*z(42)-z(6)*z(39))*z(260) - z(22)*(z(1)*z(41)+z(3)*z(39))*z(253) - 2*rad*z(210)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(285) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(304) - 0.5*z(22)*(4*e3*z(39)-z(1)*z(215)-z(3)*z(216))*z(70) - 2*z(22)*e1*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) - z(22)*e2*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e2*z(20)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - rad*z(4)*z(210)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(20)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216))*z(67) - 0.5*rad*z(4)*(4*e3*z(39)-4*e2*z(42)-z(2)*z(217)-z(3)*z(216))*z(59) - 0.5*z(219)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(41)+8*e1*z(42)-4*e4*z(39)-z(4)*z(215)-z(5)*z(217)-z(6)*z(216))*z(77) - 0.25*z(22)*z(5)*(4*z(41)*z(260)+4*e3*z(69)+4*z(39)*z(264)-2*z(215)*z(53)-2*z(216)*z(68)-z(1)*z(295)-z(3)*z(296)) - 0.25*z(4)*z(20)*(4*e3*z(69)+4*z(39)*z(264)-4*e2*z(66)-4*z(42)*z(262)-2*z(216)*z(68)-2*z(217)*z(65)-z(2)*z(297)-z(3)*z(296)) - 0.25*z(31)*(4*z(39)*z(255)+4*z(41)*z(251)+8*e1*z(66)-4*e2*z(64)-4*e4*z(69)-4*z(42)*z(253)-2*z(215)*z(67)-2*z(216)*z(52)-2*z(217)*z(70)-z(4)*z(295)-z(5)*z(297)-z(6)*z(296));
z(306) = u1*z(299) + u2*z(302) + u3*z(305) + z(237)*u1*z(80) - z(236)*u1*z(75);
z(307) = (z(235)*(z(44)*z(90)-z(48)*z(86))-z(49)*(z(44)*z(290)+z(222)*z(90)-z(48)*z(306)-z(228)*z(86)))/z(49)^2;
z(308) = m*(2*z(23)*z(218)+2*z(25)*z(224)+2*z(26)*z(220)+4*z(213)*z(209)+2*z(2)*z(14)*z(218)+4*e2*z(14)*z(25)+8*e1*z(14)*z(28)+8*rad*z(5)*z(28)*z(209)-4*z(22)*e1*z(26)-2*z(22)*e2*z(23)-2*z(5)*z(14)*z(220)-z(22)*z(2)*z(224)-8*rad*z(1)*z(28)*z(211)-8*rad*z(2)*z(25)*z(209)-4*z(213)*z(1)*z(5)*z(211));
z(309) = 4*m*(z(25)*z(218)+z(28)*z(220)) - 2*z(134)*(z(2)*z(218)+2*e2*z(25)+4*e1*z(28)-z(5)*z(220));
z(310) = z(308) + z(50)*z(281) + z(132)*z(236) + z(135)*z(237) - z(51)*z(309);
z(311) = (z(235)*(z(43)*z(90)-z(46)*z(86))-z(49)*(z(43)*z(290)+z(233)*z(90)-z(46)*z(306)-z(234)*z(86)))/z(49)^2;
z(312) = z(278) + z(139)*z(311) + z(282)*z(91) - z(138)*z(307) - z(310)*z(92);
z(313) = -z(14)*z(236) - 4*rad*z(50)*z(209);
z(314) = -4*rad*z(209) - z(22)*z(237);
z(315) = m*(z(8)*z(314)+2*e4*z(14)+2*e4*z(96)-4*e1*z(97)-4*e3*z(95)-2*z(7)*z(313)-4*rad*z(8)*z(209)-z(22)*z(9)*z(236));
z(316) = m*(4*rad*z(14)*z(209)+4*rad*z(1)*z(26)*z(211)+4*rad*z(2)*z(23)*z(209)-z(23)*z(224)-4*e1*z(14)*z(26)-2*e2*z(14)*z(23)-z(2)*z(14)*z(224)-4*rad*z(5)*z(26)*z(209)-4*rad*z(1)*z(5)*z(14)*z(211));
z(317) = z(50)*z(279) + z(128)*z(237) - 4*z(316) - z(51)*z(308) - z(127)*z(236);
z(318) = z(317) + z(50)*z(282) + z(139)*z(237) - z(51)*z(310) - z(138)*z(236);
z(319) = (z(318)*z(143)-z(141)*(9.810000000000001*z(315)+z(312)))/z(141)^2;
a12 = z(319);
z(320) = e2*(z(11)+z(9)^2/z(10)^0.5)/z(11)^2;
z(321) = e3*(z(17)+z(8)^2/z(16)^0.5)/z(17)^2;
z(322) = e2*z(9)/(z(10)^0.5*z(11)^2);
z(323) = e3*z(8)/(z(16)^0.5*z(17)^2);
z(324) = 4*rad*z(20)*z(321) + 4*z(22)*e4*z(35) + 4*e3*z(20)*z(35) + 8*e2*z(14)*z(35) + 4*rad*z(8)*z(35)*z(321) + 8*rad*z(9)*z(35)*z(320) + 4*z(213)*z(7)*(z(323)+2*z(322)) + 4*rad*z(9)*z(14)*(z(323)+2*z(322)) - 8*rad*z(14)*z(320) - 4*rad*z(35)*(z(323)+2*z(322)) - 4*rad*z(8)*z(20)*(z(323)+2*z(322));
z(325) = z(22)*z(324)/(z(37)^0.5*z(38)^2);
z(326) = (8*rad*z(38)*z(320)+z(14)*z(324)/z(37)^0.5)/z(38)^2;
z(327) = (4*rad*z(38)*z(321)-z(20)*z(324)/z(37)^0.5)/z(38)^2;
z(328) = 2*e3*z(15) - 4*rad*z(6)*z(322);
z(329) = z(328) - 2*e3*z(21) - 2*rad*z(6)*z(323);
z(330) = 2*e4*z(15) - 4*rad*z(3)*z(322);
z(331) = 2*e4*z(21) + 2*rad*z(3)*z(323) - z(330);
z(332) = 2*z(22)*e1*(z(4)*z(41)+z(6)*z(39)) + 2*e1*z(20)*(z(2)*z(42)-z(3)*z(39)) + 4*e2*z(20)*(z(5)*z(42)-z(6)*z(39)) + 2*rad*z(4)*z(321)*(z(2)*z(42)-z(3)*z(39)) + z(329)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(1)*z(20)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326)) + 0.5*z(22)*z(2)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(6)*z(326)) + 0.5*z(34)*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326)) - 2*rad*z(1)*z(321)*(z(5)*z(42)-z(6)*z(39)) - z(331)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(5)*(4*e4*z(39)-8*e2*z(41)-z(1)*z(325)-z(3)*z(326)) - 0.5*z(4)*z(20)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326)) - 0.5*z(31)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326));
z(333) = 2*e3*z(14) + 2*rad*z(9)*z(321) - 4*e2*z(20) - 4*rad*z(8)*z(320);
z(334) = 2*e1*z(15) - 4*rad*z(4)*z(322);
z(335) = z(334) - 2*e1*z(21) - 2*rad*z(4)*z(323);
z(336) = -4*e2*z(15) - 4*rad*z(1)*z(322);
z(337) = 2*rad*z(1)*z(323) - z(336) - 4*e2*z(21);
z(338) = 2*e3*z(20)*(z(1)*z(41)-z(2)*z(42)) + 4*rad*z(2)*z(320)*(z(4)*z(41)+z(6)*z(39)) + 2*rad*z(6)*z(321)*(z(1)*z(41)-z(2)*z(42)) + z(335)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(5)*z(14)*(4*e4*z(39)-8*e2*z(41)-z(1)*z(325)-z(3)*z(326)) + 0.5*z(32)*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326)) - 2*e1*z(14)*(z(4)*z(41)+z(6)*z(39)) - 4*rad*z(5)*z(320)*(z(1)*z(41)+z(3)*z(39)) - 2*e4*z(20)*(z(4)*z(41)-z(5)*z(42)) - 2*rad*z(3)*z(321)*(z(4)*z(41)-z(5)*z(42)) - z(337)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(3)*z(20)*(4*e1*z(41)-z(4)*z(325)-z(5)*z(327)) - 0.5*z(6)*z(20)*(z(1)*z(325)+z(2)*z(327)+4*e1*z(42)+8*e2*z(41)) - 0.5*z(2)*z(14)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(6)*z(326)) - 0.5*z(29)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326));
z(339) = 2*e4*z(20) - 2*z(22)*e3 - 2*rad*z(7)*z(321);
z(340) = rad*z(5)*(z(323)+2*z(322));
z(341) = 2*e1*z(15) - 4*rad*z(2)*z(322);
z(342) = 2*e1*z(21) + 2*rad*z(2)*z(323) - z(341);
z(343) = 2*z(22)*e3*(z(1)*z(41)-z(2)*z(42)) + 2*e1*z(14)*(z(2)*z(42)-z(3)*z(39)) + 4*e2*z(14)*(z(5)*z(42)-z(6)*z(39)) + 4*rad*z(1)*z(320)*(z(5)*z(42)-z(6)*z(39)) + 0.5*z(1)*z(14)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326)) + 0.5*z(33)*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326)) - 2*z(22)*e4*(z(4)*z(41)-z(5)*z(42)) - 4*rad*z(4)*z(320)*(z(2)*z(42)-z(3)*z(39)) - 2*z(340)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) - z(342)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(3)*(4*e1*z(41)-z(4)*z(325)-z(5)*z(327)) - 0.5*z(22)*z(6)*(z(1)*z(325)+z(2)*z(327)+4*e1*z(42)+8*e2*z(41)) - 0.5*z(4)*z(14)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326)) - 0.5*z(30)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326));
z(344) = 2*e4*z(14) + 4*rad*z(7)*z(320) - 4*z(22)*e2;
z(345) = z(43)*z(339) + z(48)*z(343) - z(44)*z(344) - z(46)*z(332);
z(346) = (z(345)*(z(44)*z(47)-z(45)*z(48))-z(49)*(z(44)*z(333)+z(47)*z(332)-z(45)*z(339)-z(48)*z(338)))/z(49)^2;
z(347) = (z(345)*(z(43)*z(47)-z(45)*z(46))+z(49)*(z(45)*z(344)+z(46)*z(338)-z(43)*z(333)-z(47)*z(343)))/z(49)^2;
z(348) = u1*(i22*z(346)*u3+i33*z(347)*u2-z(346)*z(101)-z(347)*z(100));
z(349) = u1*(i22*z(347)*u2+i33*z(346)*u3-z(346)*z(106)-z(347)*z(105));
z(350) = -4*rad*z(320)*u2 - z(14)*z(346)*u1;
z(351) = u2*z(350) - 4*rad*z(320)*u1^2 - z(346)*u1*z(110);
z(352) = 0.5*u3 + 0.5*e2*z(347)*u1 + 0.5*e3*z(346)*u1;
z(353) = (e1*z(347)+e4*z(346))*u1;
z(354) = 2*e2*z(353) - 4*e1*z(352) - 4*e2p;
z(355) = (z(9)*z(10)*z(354)-4*e2*z(10)*z(54)-4*e2*z(9)^2*z(54))/z(10)^1.5;
z(356) = (4*e2*z(9)*(z(11)*z(54)+2*z(9)*z(55))/z(10)^0.5+z(11)*(4*e2*z(55)-z(9)*z(355)-z(11)*z(354)))/z(11)^3;
z(357) = -u3*z(350) - z(347)*u1*z(110) - rad*u1*z(356);
z(358) = rad*(z(11)*z(355)-8*e2*z(9)*z(55)/z(10)^0.5)/z(11)^3;
z(359) = 0.5*e2*z(346)*u1 - 0.5*u2 - 0.5*e3*z(347)*u1;
z(360) = (1+e1*z(346)-e4*z(347))*u1;
z(361) = 2*e2*z(352) + 2*e3*z(359) + 2*e1p - e1*z(353) - e4*z(360);
z(362) = z(4)*z(358) + z(15)*z(361) + 2*e1*z(60) - 4*rad*z(322)*z(67);
z(363) = 2*e3*z(360) - 4*e1*z(352);
z(364) = z(15)*z(363) + z(5)*z(358) - 4*rad*z(322)*z(70);
z(365) = 2*e3p - 2*e1*z(359) - 2*e4*z(352) - e2*z(360) - e3*z(353);
z(366) = z(6)*z(358) + z(15)*z(365) + 2*e3*z(60) - 4*rad*z(322)*z(52);
z(367) = u1*z(362) + u2*z(364) + u3*z(366) + z(347)*u1*z(76) - z(346)*u1*z(71);
z(368) = rad*(u2*z(356)+z(346)*u1*z(56));
z(369) = 4*rad*z(320)*u1*u3 + z(368) - z(347)*u1*z(111);
z(370) = 2*e2*z(353) + 2*e3*z(360) - 4*e2p;
z(371) = z(1)*z(358) + z(15)*z(370) - 4*e2*z(60) - 4*rad*z(322)*z(53);
z(372) = e4*z(360) + 2*e2*z(352) + 2*e1p - 2*e3*z(359) - e1*z(353);
z(373) = z(2)*z(358) + z(15)*z(372) + 2*e1*z(60) - 4*rad*z(322)*z(65);
z(374) = 2*e2*z(359) + 2*e3*z(352) + 2*e4p - e1*z(360) - e4*z(353);
z(375) = z(3)*z(358) + z(15)*z(374) + 2*e4*z(60) - 4*rad*z(322)*z(68);
z(376) = z(346)*u1*z(73) - u1*z(371) - u2*z(373) - u3*z(375) - z(347)*u1*z(78);
z(377) = (z(22)*z(347)+4*rad*z(320))*u1;
z(378) = u2*z(350) - u1*z(377) - z(346)*u1*z(110);
z(379) = z(22)*z(346)*u1^2 - u3*z(350) - z(347)*u1*z(110) - rad*u1*z(356);
z(380) = u3*z(377) + z(368) - z(22)*z(346)*u1*u2 - z(346)*u1*z(122) - z(347)*u1*z(121);
z(381) = z(348) + z(349) + m*(z(3)*z(23)*z(351)+z(3)*z(334)*z(119)+2*e4*z(23)*z(119)+z(14)*z(357)+2*z(23)*z(367)+2*z(334)*z(116)+z(1)*z(23)*z(369)+z(1)*z(334)*z(118)+z(2)*z(23)*z(357)+z(2)*z(334)*z(120)+2*e1*z(23)*z(120)+2*z(5)*z(14)*z(376)-2*e3*z(26)*z(119)-z(6)*z(26)*z(351)-z(6)*z(336)*z(119)-4*rad*z(320)*z(120)-4*e2*z(23)*z(118)-2*e1*z(26)*z(118)-z(4)*z(26)*z(369)-z(4)*z(336)*z(118)-z(5)*z(26)*z(357)-z(5)*z(336)*z(120)-8*rad*z(5)*z(320)*z(117)) + m*(z(3)*z(23)*z(378)+z(3)*z(334)*z(124)+2*e4*z(23)*z(124)+z(14)*z(379)+z(1)*z(23)*z(380)+z(1)*z(334)*z(123)+z(2)*z(23)*z(379)+z(2)*z(334)*z(125)+2*e1*z(23)*z(125)+2*z(2)*z(14)*z(367)+4*e1*z(14)*z(116)-2*e3*z(26)*z(124)-z(6)*z(26)*z(378)-z(6)*z(336)*z(124)-2*z(26)*z(376)-2*z(336)*z(117)-4*rad*z(320)*z(125)-4*e2*z(23)*z(123)-2*e1*z(26)*z(123)-z(4)*z(26)*z(380)-z(4)*z(336)*z(123)-z(5)*z(26)*z(379)-z(5)*z(336)*z(125)-8*rad*z(2)*z(320)*z(116));
z(382) = z(347)*u1*(i33*u1-z(99));
z(383) = z(347)*u1*(i22*u1-z(99));
z(384) = m*(z(3)*z(24)*z(351)+2*e4*z(24)*z(119)+2*z(24)*z(367)+z(1)*z(24)*z(369)+z(2)*z(24)*z(357)+2*e1*z(24)*z(120)+4*rad*z(320)*z(118)+8*e2*z(14)*z(116)+8*rad*z(1)*z(320)*z(116)-2*e3*z(27)*z(119)-z(6)*z(27)*z(351)-z(6)*z(341)*z(119)-4*rad*z(3)*z(5)*z(322)*z(119)-2*z(27)*z(376)-2*z(341)*z(117)-z(14)*z(369)-4*e2*z(24)*z(118)-2*e1*z(27)*z(118)-2*z(1)*z(14)*z(367)-z(4)*z(27)*z(369)-z(4)*z(341)*z(118)-z(5)*z(27)*z(357)-z(5)*z(341)*z(120)-8*rad*z(5)*z(322)*z(116)-4*rad*z(1)*z(5)*z(322)*z(118)-4*rad*z(2)*z(5)*z(322)*z(120)) + m*(z(22)*z(378)+z(3)*z(24)*z(378)+2*e4*z(24)*z(124)+z(22)*z(3)*z(367)+z(22)*z(6)*z(376)+z(1)*z(24)*z(380)+z(2)*z(24)*z(379)+2*z(22)*e3*z(117)+2*z(22)*e4*z(116)+2*e1*z(24)*z(125)+4*rad*z(320)*z(123)+8*rad*z(4)*z(320)*z(117)-2*e3*z(27)*z(124)-z(6)*z(27)*z(378)-z(6)*z(341)*z(124)-4*rad*z(3)*z(5)*z(322)*z(124)-z(14)*z(380)-4*e1*z(14)*z(117)-4*e2*z(24)*z(123)-2*e1*z(27)*z(123)-2*z(4)*z(14)*z(376)-z(4)*z(27)*z(380)-z(4)*z(341)*z(123)-z(5)*z(27)*z(379)-z(5)*z(341)*z(125)-4*rad*z(1)*z(5)*z(322)*z(123)-4*rad*z(2)*z(5)*z(322)*z(125)) - z(382) - z(383);
z(385) = z(346)*u1*(i22*u1-z(99));
z(386) = z(346)*u1*(i33*u1-z(99));
z(387) = m*(z(3)*z(25)*z(351)+z(3)*z(328)*z(119)+2*e4*z(25)*z(119)+2*z(25)*z(367)+2*z(328)*z(116)+z(1)*z(25)*z(369)+z(1)*z(328)*z(118)+z(2)*z(25)*z(357)+z(2)*z(328)*z(120)+2*e1*z(25)*z(120)-2*e3*z(28)*z(119)-z(6)*z(28)*z(351)-z(6)*z(330)*z(119)-2*z(28)*z(376)-2*z(330)*z(117)-4*e2*z(25)*z(118)-2*e1*z(28)*z(118)-z(4)*z(28)*z(369)-z(4)*z(330)*z(118)-z(5)*z(28)*z(357)-z(5)*z(330)*z(120)) + m*(z(3)*z(25)*z(378)+z(3)*z(328)*z(124)+2*e4*z(25)*z(124)+z(1)*z(25)*z(380)+z(1)*z(328)*z(123)+z(2)*z(25)*z(379)+z(2)*z(328)*z(125)+2*e1*z(25)*z(125)-2*e3*z(28)*z(124)-z(6)*z(28)*z(378)-z(6)*z(330)*z(124)-z(22)*z(379)-4*e2*z(25)*z(123)-2*z(22)*e1*z(116)-2*e1*z(28)*z(123)-z(22)*z(2)*z(367)-z(22)*z(5)*z(376)-z(4)*z(28)*z(380)-z(4)*z(330)*z(123)-z(5)*z(28)*z(379)-z(5)*z(330)*z(125)) - z(385) - z(386);
z(388) = z(381) + z(50)*z(384) + z(347)*z(136) - z(51)*z(387) - z(346)*z(133);
z(389) = m*(2*z(24)*z(334)+2*z(26)*z(341)+2*z(27)*z(336)+z(22)*z(3)*z(334)+2*z(22)*e4*z(23)+2*z(4)*z(14)*z(336)+4*e1*z(14)*z(24)+4*e1*z(14)*z(26)+8*e2*z(14)*z(23)+8*rad*z(1)*z(23)*z(320)+8*rad*z(5)*z(27)*z(320)-2*z(22)*e3*z(26)-2*z(1)*z(14)*z(334)-2*z(5)*z(14)*z(341)-z(22)*z(6)*z(336)-8*rad*z(2)*z(24)*z(320)-8*rad*z(4)*z(26)*z(320)-8*rad*z(5)*z(23)*z(322)-8*rad*z(2)*z(5)*z(14)*z(322));
z(390) = m*(z(22)*z(6)*z(341)+2*z(22)*e3*z(27)+8*rad*z(14)*z(320)+4*z(213)*z(3)*z(5)*z(322)+8*rad*z(4)*z(27)*z(320)+8*rad*z(5)*z(24)*z(322)-2*z(27)*z(341)-8*e2*z(14)*z(24)-4*e1*z(14)*z(27)-2*z(22)*e4*z(24)-2*z(4)*z(14)*z(341)-8*rad*z(1)*z(24)*z(320)-8*rad*z(1)*z(5)*z(14)*z(322));
z(391) = m*(z(22)*z(6)*z(330)+2*z(22)*e1*z(24)+2*z(22)*e3*z(28)+2*z(1)*z(14)*z(328)+8*rad*z(4)*z(28)*z(320)+8*rad*z(5)*z(25)*z(322)-2*z(24)*z(328)-2*z(27)*z(330)-2*z(28)*z(341)-8*e2*z(14)*z(25)-4*e1*z(14)*z(28)-2*z(22)*e4*z(25)-2*z(4)*z(14)*z(330)-z(22)*z(3)*z(328)-z(22)*z(5)*z(341)-8*rad*z(1)*z(25)*z(320)-4*z(213)*z(2)*z(5)*z(322));
z(392) = z(389) + z(51)*z(391) - 2*z(50)*z(390) - z(131)*z(346) - z(132)*z(347);
z(393) = 2*e1*z(359) + 2*e4*z(352) + 2*e3p - e2*z(360) - e3*z(353);
z(394) = (z(8)*z(16)*z(393)+2*e3*z(16)*z(57)+2*e3*z(8)^2*z(57))/z(16)^1.5;
z(395) = (z(17)*(z(8)*z(394)+z(17)*z(393)+2*e3*z(58))+2*e3*z(8)*(z(17)*z(57)+2*z(8)*z(58))/z(16)^0.5)/z(17)^3;
z(396) = z(14)*z(393) + z(20)*z(354) + rad*z(9)*z(395) + 2*rad*e3*z(56) + 2*rad*z(321)*z(54) - 4*rad*e2*z(59) - 4*rad*z(320)*z(57) - rad*z(8)*z(356);
z(397) = e4*z(353) + 2*e3*z(352) - 2*e2*z(359) - e1*z(360) - 2*e4p;
z(398) = z(22)*z(354) + rad*z(7)*z(356) + 2*rad*e4*z(56) + 4*rad*z(320)*z(62) - z(14)*z(397);
z(399) = 2*rad*e4*z(59) - z(22)*z(393) - z(20)*z(397) - 2*rad*z(321)*z(62) - rad*z(7)*z(395);
z(400) = u1*z(396) + u2*z(398) + u3*z(399) + z(347)*u1*z(89) - z(346)*u1*z(87);
z(401) = rad*(z(17)*z(394)+4*e3*z(8)*z(58)/z(16)^0.5)/z(17)^3;
z(402) = z(362) - 2*e1*z(61) - z(4)*z(401) - z(21)*z(361) - 2*rad*z(323)*z(67);
z(403) = 2*rad*z(20)*z(395) + 2*z(20)*z(35)*z(393) + 4*z(293)*z(321)*z(59) + 2*rad*z(8)*z(35)*z(395) + 2*rad*z(9)*z(35)*z(356) + 4*rad*e3*z(35)*z(59) + 4*rad*z(35)*z(321)*z(57) + 8*rad*e2*z(35)*z(56) + 8*rad*z(35)*z(320)*z(54) + 2*z(35)*(z(358)-z(401)) + 2*z(8)*z(20)*(z(358)-z(401)) + 4*z(22)*e4*(z(60)-z(61)) + 4*z(213)*(z(323)+2*z(322))*z(62) + 4*e3*z(20)*(z(60)-z(61)) + 8*e2*z(14)*(z(60)-z(61)) + 4*rad*z(8)*z(321)*(z(60)-z(61)) + 4*rad*z(14)*(z(323)+2*z(322))*z(54) + 4*z(293)*z(9)*(z(323)+2*z(322))*z(56) + 8*rad*z(9)*z(320)*(z(60)-z(61)) - 8*z(293)*z(320)*z(56) - 2*rad*z(14)*z(356) - 2*z(22)*z(35)*z(397) - 2*z(14)*z(35)*z(354) - 2*z(22)*z(7)*(z(358)-z(401)) - 2*z(9)*z(14)*(z(358)-z(401)) - 4*rad*z(20)*(z(323)+2*z(322))*z(57) - 4*z(293)*z(8)*(z(323)+2*z(322))*z(59) - 4*rad*(z(323)+2*z(322))*(z(60)-z(61));
z(404) = z(22)*(2*z(403)/(z(37)^0.5*z(38)^2)-z(324)*(2+z(38)/z(37)^0.5)*z(63)/(z(37)*z(38)^3));
z(405) = (2*z(324)*(rad*z(38)*z(56)-z(14)*z(63)/z(37)^0.5)/z(37)^0.5+z(38)*(4*rad*z(38)*z(356)+(2*z(14)*z(37)*z(403)-z(14)*z(324)*z(63)-8*rad*z(37)*z(320)*z(63))/z(37)^1.5))/z(38)^3;
z(406) = (2*z(20)*z(324)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(395)+(z(20)*z(324)*z(63)-2*z(20)*z(37)*z(403)-4*rad*z(37)*z(321)*z(63)-2*rad*z(37)*z(324)*z(59))/z(37)^1.5))/z(38)^3;
z(407) = z(1)*z(401) + z(21)*z(370) + 2*rad*z(323)*z(53) - z(371) - 4*e2*z(61);
z(408) = z(14)*(z(1)*z(41)+z(3)*z(39))*z(363) + rad*z(2)*(z(4)*z(41)+z(6)*z(39))*z(356) + z(20)*(z(1)*z(41)-z(2)*z(42))*z(365) + 4*rad*z(320)*(z(4)*z(41)+z(6)*z(39))*z(65) + rad*z(6)*(z(1)*z(41)-z(2)*z(42))*z(395) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(402) + 2*rad*e3*(z(1)*z(41)-z(2)*z(42))*z(59) + 2*rad*z(321)*(z(1)*z(41)-z(2)*z(42))*z(52) + e1*z(14)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) + e3*z(20)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*z(14)*(4*e4*z(39)-8*e2*z(41)-z(1)*z(325)-z(3)*z(326))*z(70) + 2*rad*z(5)*z(320)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) + rad*z(6)*z(321)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*rad*z(5)*(4*e4*z(39)-8*e2*z(41)-z(1)*z(325)-z(3)*z(326))*z(56) + 0.5*z(337)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326))*z(84) + 0.25*z(5)*z(14)*(4*e4*z(69)+4*z(39)*z(374)+4*z(41)*z(370)+8*e2*z(64)-2*z(325)*z(53)-2*z(326)*z(68)-z(1)*z(404)-z(3)*z(405)) + 0.25*z(6)*z(20)*(4*z(41)*z(370)+8*e2*z(64)-4*e1*z(66)-4*z(42)*z(372)-2*z(325)*z(53)-2*z(327)*z(65)-z(1)*z(404)-z(2)*z(406)) + 0.25*z(32)*(4*e4*z(69)+4*z(39)*z(374)+4*z(41)*z(370)+8*e2*z(64)-4*e1*z(66)-4*z(42)*z(372)-2*z(325)*z(53)-2*z(326)*z(68)-2*z(327)*z(65)-z(1)*z(404)-z(2)*z(406)-z(3)*z(405)) - z(14)*(z(4)*z(41)+z(6)*z(39))*z(372) - 4*rad*z(320)*(z(1)*z(41)+z(3)*z(39))*z(70) - 2*rad*e1*(z(4)*z(41)+z(6)*z(39))*z(56) - rad*z(5)*(z(1)*z(41)+z(3)*z(39))*z(356) - z(20)*(z(4)*z(41)-z(5)*z(42))*z(374) - 2*rad*e4*(z(4)*z(41)-z(5)*z(42))*z(59) - 2*rad*z(321)*(z(4)*z(41)-z(5)*z(42))*z(68) - rad*z(3)*(z(4)*z(41)-z(5)*z(42))*z(395) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(407) - 0.5*z(20)*(4*e1*z(41)-z(4)*z(325)-z(5)*z(327))*z(68) - 0.5*rad*z(3)*(4*e1*z(41)-z(4)*z(325)-z(5)*z(327))*z(59) - 0.5*z(20)*(z(1)*z(325)+z(2)*z(327)+4*e1*z(42)+8*e2*z(41))*z(52) - 0.5*rad*z(6)*(z(1)*z(325)+z(2)*z(327)+4*e1*z(42)+8*e2*z(41))*z(59) - 2*rad*z(2)*z(320)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e4*z(20)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(14)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(6)*z(326))*z(65) - rad*z(3)*z(321)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*rad*z(2)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(6)*z(326))*z(56) - 0.5*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326))*z(82) - 0.5*z(335)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.25*z(3)*z(20)*(4*z(41)*z(361)-4*z(42)*z(363)-4*e1*z(64)-2*z(325)*z(67)-2*z(327)*z(70)-z(4)*z(404)-z(5)*z(406)) - 0.25*z(2)*z(14)*(4*e3*z(69)+4*z(39)*z(365)+4*z(41)*z(361)-4*e1*z(64)-2*z(325)*z(67)-2*z(326)*z(52)-z(4)*z(404)-z(6)*z(405)) - 0.25*z(29)*(4*e3*z(69)+4*z(39)*z(365)+4*z(41)*z(361)-4*z(42)*z(363)-4*e1*z(64)-2*z(325)*z(67)-2*z(326)*z(52)-2*z(327)*z(70)-z(4)*z(404)-z(5)*z(406)-z(6)*z(405));
z(409) = z(364) - z(21)*z(363) - z(5)*z(401) - 2*rad*z(323)*z(70);
z(410) = z(2)*z(401) + z(21)*z(372) + 2*e1*z(61) + 2*rad*z(323)*z(65) - z(373);
z(411) = z(22)*(z(1)*z(41)-z(2)*z(42))*z(365) + z(14)*(z(2)*z(42)-z(3)*z(39))*z(361) + rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(356) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(409) + 2*rad*e1*(z(2)*z(42)-z(3)*z(39))*z(56) + 4*rad*e2*(z(5)*z(42)-z(6)*z(39))*z(56) + 4*rad*z(320)*(z(5)*z(42)-z(6)*z(39))*z(53) + 0.5*z(14)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326))*z(53) + 0.5*rad*z(1)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326))*z(56) + z(22)*e3*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 2*rad*z(4)*z(320)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) + z(340)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) + 0.5*z(342)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326))*z(74) + 0.25*z(1)*z(14)*(4*e3*z(69)+4*z(39)*z(365)-4*z(42)*z(363)-2*z(326)*z(52)-2*z(327)*z(70)-z(5)*z(406)-z(6)*z(405)) + 0.25*z(22)*z(6)*(4*z(41)*z(370)+8*e2*z(64)-4*e1*z(66)-4*z(42)*z(372)-2*z(325)*z(53)-2*z(327)*z(65)-z(1)*z(404)-z(2)*z(406)) + 0.25*z(33)*(4*e4*z(69)+4*z(39)*z(374)+4*z(41)*z(370)+8*e2*z(64)-4*e1*z(66)-4*z(42)*z(372)-2*z(325)*z(53)-2*z(326)*z(68)-2*z(327)*z(65)-z(1)*z(404)-z(2)*z(406)-z(3)*z(405)) - z(22)*(z(4)*z(41)-z(5)*z(42))*z(374) - z(14)*(z(5)*z(42)-z(6)*z(39))*z(370) - 4*rad*z(320)*(z(2)*z(42)-z(3)*z(39))*z(67) - rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(356) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(410) - 0.5*z(22)*(4*e1*z(41)-z(4)*z(325)-z(5)*z(327))*z(68) - 0.5*z(22)*(z(1)*z(325)+z(2)*z(327)+4*e1*z(42)+8*e2*z(41))*z(52) - 2*e2*z(14)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - e1*z(14)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 2*rad*z(1)*z(320)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - z(22)*e4*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(14)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326))*z(67) - 0.5*rad*z(4)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326))*z(56) - 0.5*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326))*z(72) - 0.25*z(22)*z(3)*(4*z(41)*z(361)-4*z(42)*z(363)-4*e1*z(64)-2*z(325)*z(67)-2*z(327)*z(70)-z(4)*z(404)-z(5)*z(406)) - 0.25*z(4)*z(14)*(4*e4*z(69)+4*z(39)*z(374)-4*e1*z(66)-4*z(42)*z(372)-2*z(326)*z(68)-2*z(327)*z(65)-z(2)*z(406)-z(3)*z(405)) - 0.25*z(30)*(4*e3*z(69)+4*z(39)*z(365)+4*z(41)*z(361)-4*z(42)*z(363)-4*e1*z(64)-2*z(325)*z(67)-2*z(326)*z(52)-2*z(327)*z(70)-z(4)*z(404)-z(5)*z(406)-z(6)*z(405));
z(412) = z(366) - 2*e3*z(61) - z(6)*z(401) - z(21)*z(365) - 2*rad*z(323)*z(52);
z(413) = z(3)*z(401) + z(21)*z(374) + 2*e4*z(61) + 2*rad*z(323)*z(68) - z(375);
z(414) = z(22)*(z(4)*z(41)+z(6)*z(39))*z(372) + z(20)*(z(2)*z(42)-z(3)*z(39))*z(361) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(395) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(412) + 2*rad*e1*(z(2)*z(42)-z(3)*z(39))*z(59) + 2*rad*z(321)*(z(2)*z(42)-z(3)*z(39))*z(67) + 4*rad*e2*(z(5)*z(42)-z(6)*z(39))*z(59) + 0.5*z(20)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326))*z(53) + 0.5*rad*z(1)*(4*e3*z(39)-z(5)*z(327)-z(6)*z(326))*z(59) + rad*z(1)*z(321)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(22)*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(6)*z(326))*z(65) + 0.5*z(331)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e4*z(39)-8*e2*z(41)-4*e1*z(42)-z(1)*z(325)-z(2)*z(327)-z(3)*z(326))*z(79) + 0.25*z(1)*z(20)*(4*e3*z(69)+4*z(39)*z(365)-4*z(42)*z(363)-2*z(326)*z(52)-2*z(327)*z(70)-z(5)*z(406)-z(6)*z(405)) + 0.25*z(22)*z(2)*(4*e3*z(69)+4*z(39)*z(365)+4*z(41)*z(361)-4*e1*z(64)-2*z(325)*z(67)-2*z(326)*z(52)-z(4)*z(404)-z(6)*z(405)) + 0.25*z(34)*(4*e4*z(69)+4*z(39)*z(374)+4*z(41)*z(370)+8*e2*z(64)-4*e1*z(66)-4*z(42)*z(372)-2*z(325)*z(53)-2*z(326)*z(68)-2*z(327)*z(65)-z(1)*z(404)-z(2)*z(406)-z(3)*z(405)) - z(22)*(z(1)*z(41)+z(3)*z(39))*z(363) - z(20)*(z(5)*z(42)-z(6)*z(39))*z(370) - 2*rad*z(321)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(395) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(413) - 2*e2*z(20)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - z(22)*e1*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e1*z(20)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - rad*z(4)*z(321)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(22)*(4*e4*z(39)-8*e2*z(41)-z(1)*z(325)-z(3)*z(326))*z(70) - 0.5*z(20)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326))*z(67) - 0.5*rad*z(4)*(4*e4*z(39)-4*e1*z(42)-z(2)*z(327)-z(3)*z(326))*z(59) - 0.5*(4*e1*z(41)+4*e3*z(39)-z(4)*z(325)-z(5)*z(327)-z(6)*z(326))*z(77) - 0.5*z(329)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.25*z(22)*z(5)*(4*e4*z(69)+4*z(39)*z(374)+4*z(41)*z(370)+8*e2*z(64)-2*z(325)*z(53)-2*z(326)*z(68)-z(1)*z(404)-z(3)*z(405)) - 0.25*z(4)*z(20)*(4*e4*z(69)+4*z(39)*z(374)-4*e1*z(66)-4*z(42)*z(372)-2*z(326)*z(68)-2*z(327)*z(65)-z(2)*z(406)-z(3)*z(405)) - 0.25*z(31)*(4*e3*z(69)+4*z(39)*z(365)+4*z(41)*z(361)-4*z(42)*z(363)-4*e1*z(64)-2*z(325)*z(67)-2*z(326)*z(52)-2*z(327)*z(70)-z(4)*z(404)-z(5)*z(406)-z(6)*z(405));
z(415) = u1*z(408) + u2*z(411) + u3*z(414) + z(347)*u1*z(80) - z(346)*u1*z(75);
z(416) = (z(345)*(z(44)*z(90)-z(48)*z(86))-z(49)*(z(44)*z(400)+z(332)*z(90)-z(48)*z(415)-z(339)*z(86)))/z(49)^2;
z(417) = m*(z(22)*z(2)*z(334)+2*z(22)*e1*z(23)+2*z(5)*z(14)*z(330)+8*rad*z(2)*z(25)*z(320)-4*z(213)*z(320)-2*z(23)*z(328)-2*z(25)*z(334)-2*z(26)*z(330)-2*z(28)*z(336)-4*e1*z(14)*z(25)-2*z(2)*z(14)*z(328)-z(22)*z(5)*z(336)-8*rad*z(5)*z(28)*z(320));
z(418) = 4*m*(z(25)*z(328)+z(28)*z(330)) - 2*z(134)*(z(2)*z(328)+2*e1*z(25)-z(5)*z(330));
z(419) = z(132)*z(346) + z(135)*z(347) - z(417) - z(50)*z(391) - z(51)*z(418);
z(420) = (z(345)*(z(43)*z(90)-z(46)*z(86))-z(49)*(z(43)*z(400)+z(343)*z(90)-z(46)*z(415)-z(344)*z(86)))/z(49)^2;
z(421) = z(388) + z(139)*z(420) + z(392)*z(91) - z(138)*z(416) - z(419)*z(92);
z(422) = -z(14)*z(346) - 4*rad*z(50)*z(320);
z(423) = -4*rad*z(320) - z(22)*z(347);
z(424) = m*(z(8)*z(423)+2*e3*z(14)+2*e3*z(96)+4*e4*z(95)-4*e2*z(97)-2*z(7)*z(422)-4*rad*z(8)*z(320)-z(22)*z(9)*z(346));
z(425) = m*(z(5)*z(14)*z(336)+4*rad*z(14)*z(320)+4*rad*z(2)*z(23)*z(320)-z(23)*z(334)-z(26)*z(336)-2*e1*z(14)*z(23)-z(2)*z(14)*z(334)-4*rad*z(5)*z(26)*z(320));
z(426) = z(50)*z(389) + z(51)*z(417) + z(128)*z(347) - 4*z(425) - z(127)*z(346);
z(427) = z(426) + z(50)*z(392) + z(139)*z(347) - z(51)*z(419) - z(138)*z(346);
z(428) = (z(427)*z(143)-z(141)*(9.810000000000001*z(424)+z(421)))/z(141)^2;
a13 = z(428);
z(429) = e2*(z(17)+z(8)^2/z(16)^0.5)/z(17)^2;
z(430) = e2*z(8)/(z(16)^0.5*z(17)^2);
z(431) = 4*rad*z(20)*z(429) + 4*z(213)*z(7)*z(430) + 4*e2*z(20)*z(35) + 4*rad*z(8)*z(35)*z(429) + 4*rad*z(9)*z(14)*z(430) - 4*rad*z(35)*z(430) - 4*z(22)*e1*z(35) - 4*rad*z(8)*z(20)*z(430);
z(432) = z(22)*z(431)/(z(37)^0.5*z(38)^2);
z(433) = z(14)*z(431)/(z(37)^0.5*z(38)^2);
z(434) = (4*rad*z(38)*z(429)-z(20)*z(431)/z(37)^0.5)/z(38)^2;
z(435) = 2*e2*z(15) - 2*e2*z(21) - 2*rad*z(6)*z(430);
z(436) = 2*e1*z(21) + 2*rad*z(3)*z(430) - 2*e1*z(15);
z(437) = 4*z(22)*e3*(z(1)*z(41)+z(3)*z(39)) + 2*e4*z(20)*(z(2)*z(42)-z(3)*z(39)) + 4*e3*z(20)*(z(5)*z(42)-z(6)*z(39)) + 2*rad*z(4)*z(429)*(z(2)*z(42)-z(3)*z(39)) + z(435)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(22)*z(2)*(4*e2*z(39)+4*e4*z(41)-z(4)*z(432)-z(6)*z(433)) + 0.5*z(1)*z(20)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433)) + 0.5*z(34)*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433)) - 2*z(22)*e4*(z(4)*z(41)+z(6)*z(39)) - 2*rad*z(1)*z(429)*(z(5)*z(42)-z(6)*z(39)) - z(436)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(5)*(4*e1*z(39)-8*e3*z(41)-z(1)*z(432)-z(3)*z(433)) - 0.5*z(4)*z(20)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433)) - 0.5*z(31)*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433));
z(438) = 2*e2*z(14) + 2*rad*z(9)*z(429);
z(439) = 2*e4*z(15) - 2*e4*z(21) - 2*rad*z(4)*z(430);
z(440) = 4*e3*z(15) + 2*rad*z(1)*z(430) - 4*e3*z(21);
z(441) = 2*e4*z(14)*(z(4)*z(41)+z(6)*z(39)) + 2*e2*z(20)*(z(1)*z(41)-z(2)*z(42)) + 2*rad*z(6)*z(429)*(z(1)*z(41)-z(2)*z(42)) + z(439)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(5)*z(14)*(4*e1*z(39)-8*e3*z(41)-z(1)*z(432)-z(3)*z(433)) + 0.5*z(6)*z(20)*(4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)) + 0.5*z(32)*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433)) - 4*e3*z(14)*(z(1)*z(41)+z(3)*z(39)) - 2*e1*z(20)*(z(4)*z(41)-z(5)*z(42)) - 2*rad*z(3)*z(429)*(z(4)*z(41)-z(5)*z(42)) - z(440)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(2)*z(14)*(4*e2*z(39)+4*e4*z(41)-z(4)*z(432)-z(6)*z(433)) - 0.5*z(3)*z(20)*(4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)) - 0.5*z(29)*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433));
z(442) = -2*z(22)*e2 - 2*e1*z(20) - 2*rad*z(7)*z(429);
z(443) = 4*e3*z(21) - 4*e3*z(15) - 2*rad*z(5)*z(430);
z(444) = 2*e4*z(15) + 2*rad*z(2)*z(430) - 2*e4*z(21);
z(445) = 2*z(22)*e2*(z(1)*z(41)-z(2)*z(42)) + 2*e4*z(14)*(z(2)*z(42)-z(3)*z(39)) + 4*e3*z(14)*(z(5)*z(42)-z(6)*z(39)) + z(443)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(22)*z(6)*(4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)) + 0.5*z(1)*z(14)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433)) + 0.5*z(33)*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433)) - 2*z(22)*e1*(z(4)*z(41)-z(5)*z(42)) - z(444)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(22)*z(3)*(4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)) - 0.5*z(4)*z(14)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433)) - 0.5*z(30)*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433));
z(446) = z(43)*z(442) + z(48)*z(445) + 2*e1*z(14)*z(44) - z(46)*z(437);
z(447) = (z(446)*(z(44)*z(47)-z(45)*z(48))-z(49)*(z(44)*z(438)+z(47)*z(437)-z(45)*z(442)-z(48)*z(441)))/z(49)^2;
z(448) = (z(446)*(z(43)*z(47)-z(45)*z(46))+z(49)*(z(46)*z(441)-z(43)*z(438)-z(47)*z(445)-2*e1*z(14)*z(45)))/z(49)^2;
z(449) = u1*(i22*z(447)*u3+i33*z(448)*u2-z(447)*z(101)-z(448)*z(100));
z(450) = u1*(i22*z(448)*u2+i33*z(447)*u3-z(447)*z(106)-z(448)*z(105));
z(451) = z(447)*u1*(z(14)*u2+z(110));
z(452) = 0.5*e2*z(448)*u1 + 0.5*e3*z(447)*u1 - 0.5*u2;
z(453) = (1-e1*z(448)-e4*z(447))*u1;
z(454) = -4*e1*z(452) - 2*e2*z(453);
z(455) = (z(11)+z(9)^2/z(10)^0.5)*z(454)/z(11)^2;
z(456) = u1*(z(14)*z(447)*u3+rad*z(455)-z(448)*z(110));
z(457) = rad*z(9)*z(454)/(z(10)^0.5*z(11)^2);
z(458) = 0.5*e2*z(447)*u1 - 0.5*u3 - 0.5*e3*z(448)*u1;
z(459) = (e1*z(447)-e4*z(448))*u1;
z(460) = e1*z(453) + 2*e2*z(452) + 2*e3*z(458) + 2*e4p - e4*z(459);
z(461) = z(4)*z(457) + z(15)*z(460) + 2*e4*z(60);
z(462) = 2*e3*z(459) - 4*e1*z(452) - 4*e3p;
z(463) = z(5)*z(457) + z(15)*z(462) - 4*e3*z(60);
z(464) = e3*z(453) + 2*e2p - 2*e1*z(458) - 2*e4*z(452) - e2*z(459);
z(465) = z(6)*z(457) + z(15)*z(464) + 2*e2*z(60);
z(466) = u1*z(461) + u2*z(463) + u3*z(465) + z(448)*u1*z(76) - z(447)*u1*z(71);
z(467) = rad*(u2*z(455)-z(447)*u1*z(56));
z(468) = -z(448)*u1*z(111) - z(467);
z(469) = 2*e3*z(459) - 2*e2*z(453) - 4*e3p;
z(470) = z(1)*z(457) + z(15)*z(469) - 4*e3*z(60);
z(471) = e1*z(453) + e4*z(459) + 2*e2*z(452) - 2*e3*z(458) - 2*e4p;
z(472) = z(2)*z(457) + z(15)*z(471) - 2*e4*z(60);
z(473) = e4*z(453) + 2*e2*z(458) + 2*e3*z(452) + 2*e1p - e1*z(459);
z(474) = z(3)*z(457) + z(15)*z(473) + 2*e1*z(60);
z(475) = z(447)*u1*z(73) - u1*z(470) - u2*z(472) - u3*z(474) - z(448)*u1*z(78);
z(476) = u1*(z(22)*z(448)*u1+z(14)*z(447)*u2+z(447)*z(110));
z(477) = u1*(z(22)*z(447)*u1+z(14)*z(447)*u3+rad*z(455)-z(448)*z(110));
z(478) = z(22)*z(448)*u1*u3 - z(22)*z(447)*u1*u2 - z(447)*u1*z(122) - z(448)*u1*z(121) - z(467);
z(479) = z(449) + z(450) - m*(z(3)*z(23)*z(451)+z(5)*z(26)*z(456)+2*e2*z(26)*z(119)+z(4)*z(26)*z(468)+2*e4*z(23)*z(120)+2*e4*z(26)*z(118)+4*e3*z(23)*z(118)+8*e3*z(14)*z(117)-z(14)*z(456)-2*e1*z(23)*z(119)-z(2)*z(23)*z(456)-z(6)*z(26)*z(451)-4*e3*z(6)*z(15)*z(119)-2*e4*z(3)*z(15)*z(119)-2*z(23)*z(466)-4*e3*z(26)*z(120)-4*e4*z(15)*z(116)-2*z(5)*z(14)*z(475)-z(1)*z(23)*z(468)-4*e3*z(4)*z(15)*z(118)-4*e3*z(5)*z(15)*z(120)-2*e4*z(1)*z(15)*z(118)-2*e4*z(2)*z(15)*z(120)) - m*(z(3)*z(23)*z(476)+z(5)*z(26)*z(477)+2*e2*z(26)*z(124)+2*z(26)*z(475)+z(4)*z(26)*z(478)+2*e4*z(23)*z(125)+2*e4*z(26)*z(123)+4*e3*z(23)*z(123)+4*e4*z(14)*z(116)-z(14)*z(477)-2*e1*z(23)*z(124)-z(2)*z(23)*z(477)-z(6)*z(26)*z(476)-4*e3*z(6)*z(15)*z(124)-2*e4*z(3)*z(15)*z(124)-8*e3*z(15)*z(117)-4*e3*z(26)*z(125)-2*z(2)*z(14)*z(466)-z(1)*z(23)*z(478)-4*e3*z(4)*z(15)*z(123)-4*e3*z(5)*z(15)*z(125)-2*e4*z(1)*z(15)*z(123)-2*e4*z(2)*z(15)*z(125));
z(480) = z(448)*u1*(i33*u1-z(99));
z(481) = z(448)*u1*(i22*u1-z(99));
z(482) = -z(480) - z(481) - m*(z(3)*z(24)*z(451)+z(5)*z(27)*z(456)+2*e2*z(27)*z(119)+4*e3*z(3)*z(15)*z(119)+z(14)*z(468)+2*z(27)*z(475)+z(4)*z(27)*z(468)+2*e4*z(24)*z(120)+2*e4*z(27)*z(118)+2*z(1)*z(14)*z(466)+4*e3*z(24)*z(118)+8*e3*z(15)*z(116)+4*e3*z(1)*z(15)*z(118)+4*e3*z(2)*z(15)*z(120)-2*e1*z(24)*z(119)-z(2)*z(24)*z(456)-z(6)*z(27)*z(451)-2*e4*z(6)*z(15)*z(119)-2*z(24)*z(466)-8*e3*z(14)*z(116)-4*e3*z(27)*z(120)-4*e4*z(15)*z(117)-z(1)*z(24)*z(468)-2*e4*z(4)*z(15)*z(118)-2*e4*z(5)*z(15)*z(120)) - m*(z(22)*z(476)+z(3)*z(24)*z(476)+z(5)*z(27)*z(477)+2*e2*z(27)*z(124)+4*e3*z(3)*z(15)*z(124)+z(14)*z(478)+z(4)*z(27)*z(478)+2*e4*z(24)*z(125)+2*e4*z(27)*z(123)+2*z(4)*z(14)*z(475)+4*e3*z(24)*z(123)+4*e4*z(14)*z(117)+4*e3*z(1)*z(15)*z(123)+4*e3*z(2)*z(15)*z(125)-2*e1*z(24)*z(124)-z(2)*z(24)*z(477)-z(6)*z(27)*z(476)-2*e4*z(6)*z(15)*z(124)-4*e3*z(27)*z(125)-2*z(22)*e1*z(116)-2*z(22)*e2*z(117)-z(22)*z(3)*z(466)-z(22)*z(6)*z(475)-z(1)*z(24)*z(478)-2*e4*z(4)*z(15)*z(123)-2*e4*z(5)*z(15)*z(125));
z(483) = z(447)*u1*(i22*u1-z(99));
z(484) = z(447)*u1*(i33*u1-z(99));
z(485) = -z(483) - z(484) - m*(z(3)*z(25)*z(451)+z(5)*z(28)*z(456)+2*e2*z(28)*z(119)+2*e1*z(6)*z(15)*z(119)+2*z(28)*z(475)+z(4)*z(28)*z(468)+2*e4*z(25)*z(120)+2*e4*z(28)*z(118)+4*e1*z(15)*z(117)+4*e3*z(25)*z(118)+2*e1*z(4)*z(15)*z(118)+2*e1*z(5)*z(15)*z(120)-2*e1*z(25)*z(119)-z(2)*z(25)*z(456)-z(6)*z(28)*z(451)-2*e2*z(3)*z(15)*z(119)-2*z(25)*z(466)-4*e2*z(15)*z(116)-4*e3*z(28)*z(120)-z(1)*z(25)*z(468)-2*e2*z(1)*z(15)*z(118)-2*e2*z(2)*z(15)*z(120)) - m*(z(22)*z(477)+z(3)*z(25)*z(476)+z(5)*z(28)*z(477)+2*e2*z(28)*z(124)+2*e1*z(6)*z(15)*z(124)+z(22)*z(2)*z(466)+z(22)*z(5)*z(475)+z(4)*z(28)*z(478)+2*e4*z(25)*z(125)+2*e4*z(28)*z(123)+4*e3*z(25)*z(123)+2*e1*z(4)*z(15)*z(123)+2*e1*z(5)*z(15)*z(125)-2*e1*z(25)*z(124)-z(2)*z(25)*z(477)-z(6)*z(28)*z(476)-2*e2*z(3)*z(15)*z(124)-4*z(22)*e3*z(117)-4*e3*z(28)*z(125)-2*z(22)*e4*z(116)-z(1)*z(25)*z(478)-2*e2*z(1)*z(15)*z(123)-2*e2*z(2)*z(15)*z(125));
z(486) = z(479) + z(50)*z(482) + z(448)*z(136) - z(51)*z(485) - z(447)*z(133);
z(487) = m*(z(22)*e2*z(26)+2*e4*z(14)*z(24)+2*e4*z(15)*z(26)+4*e3*z(15)*z(23)+4*e3*z(15)*z(27)+2*e4*z(1)*z(14)*z(15)+4*e3*z(2)*z(14)*z(15)+4*e3*z(4)*z(14)*z(15)-4*e3*z(14)*z(23)-4*e3*z(14)*z(27)-2*e4*z(14)*z(26)-2*e4*z(15)*z(24)-z(22)*e1*z(23)-2*z(22)*e3*z(6)*z(15)-2*e4*z(5)*z(14)*z(15)-z(22)*e4*z(3)*z(15));
z(488) = m*(z(22)*e1*z(24)+2*e4*z(14)*z(27)+4*e3*z(14)*z(24)+z(22)*e4*z(6)*z(15)+4*e3*z(1)*z(14)*z(15)-4*e3*z(15)*z(24)-2*e4*z(15)*z(27)-z(22)*e2*z(27)-2*z(22)*e3*z(3)*z(15)-2*e4*z(4)*z(14)*z(15));
z(489) = m*(z(22)*e1*z(25)+z(22)*e4*z(24)+2*e1*z(15)*z(27)+2*e2*z(15)*z(24)+2*e4*z(14)*z(28)+4*e3*z(14)*z(25)+z(22)*e2*z(3)*z(15)+2*z(22)*e3*z(2)*z(15)+2*e1*z(4)*z(14)*z(15)-4*e3*z(15)*z(25)-2*z(22)*e3*z(27)-2*e4*z(15)*z(28)-z(22)*e2*z(28)-2*e2*z(1)*z(14)*z(15)-z(22)*e1*z(6)*z(15)-z(22)*e4*z(5)*z(15));
z(490) = 4*z(50)*z(488) - 2*z(487) - 2*z(51)*z(489) - z(131)*z(447) - z(132)*z(448);
z(491) = e3*z(453) + 2*e1*z(458) + 2*e4*z(452) + 2*e2p - e2*z(459);
z(492) = (z(8)*z(16)*z(491)+2*e2*z(16)*z(57)+2*e2*z(8)^2*z(57))/z(16)^1.5;
z(493) = (z(17)*(z(8)*z(492)+z(17)*z(491)+2*e2*z(58))+2*e2*z(8)*(z(17)*z(57)+2*z(8)*z(58))/z(16)^0.5)/z(17)^3;
z(494) = z(20)*z(454) + rad*z(8)*z(455) + z(14)*z(491) + rad*z(9)*z(493) + 2*rad*e2*z(56) + 2*rad*z(429)*z(54);
z(495) = 2*e3*z(452) + 2*e1p - 2*e2*z(458) - e1*z(459) - e4*z(453);
z(496) = z(22)*z(454) - rad*z(7)*z(455) - z(14)*z(495) - 2*rad*e1*z(56);
z(497) = -z(22)*z(491) - z(20)*z(495) - 2*rad*e1*z(59) - 2*rad*z(429)*z(62) - rad*z(7)*z(493);
z(498) = u1*z(494) + u2*z(496) + u3*z(497) + z(448)*u1*z(89) - z(447)*u1*z(87);
z(499) = rad*(z(17)*z(492)+4*e2*z(8)*z(58)/z(16)^0.5)/z(17)^3;
z(500) = z(461) - 2*e4*z(61) - z(4)*z(499) - z(21)*z(460) - 2*rad*z(430)*z(67);
z(501) = 2*rad*z(14)*z(455) + 2*rad*z(20)*z(493) + 2*z(20)*z(35)*z(491) + 4*z(213)*z(430)*z(62) + 4*z(293)*z(429)*z(59) + 2*rad*z(8)*z(35)*z(493) + 4*rad*e2*z(35)*z(59) + 4*rad*z(14)*z(430)*z(54) + 4*rad*z(35)*z(429)*z(57) + 4*z(293)*z(9)*z(430)*z(56) + 2*z(35)*(z(457)-z(499)) + 2*z(8)*z(20)*(z(457)-z(499)) + 4*e2*z(20)*(z(60)-z(61)) + 4*rad*z(8)*z(429)*(z(60)-z(61)) - 2*z(14)*z(35)*z(454) - 2*rad*z(9)*z(35)*z(455) - 2*z(22)*z(35)*z(495) - 4*rad*z(20)*z(430)*z(57) - 4*z(293)*z(8)*z(430)*z(59) - 4*rad*z(430)*(z(60)-z(61)) - 4*z(22)*e1*(z(60)-z(61)) - 2*z(22)*z(7)*(z(457)-z(499)) - 2*z(9)*z(14)*(z(457)-z(499));
z(502) = z(22)*(2*z(501)/(z(37)^0.5*z(38)^2)-z(431)*(2+z(38)/z(37)^0.5)*z(63)/(z(37)*z(38)^3));
z(503) = (2*z(14)*z(431)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(455)+(z(14)*z(431)*z(63)-2*z(14)*z(37)*z(501)-2*rad*z(37)*z(431)*z(56))/z(37)^1.5))/z(38)^3;
z(504) = (2*z(20)*z(431)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(493)+(z(20)*z(431)*z(63)-2*z(20)*z(37)*z(501)-4*rad*z(37)*z(429)*z(63)-2*rad*z(37)*z(431)*z(59))/z(37)^1.5))/z(38)^3;
z(505) = z(1)*z(499) + z(21)*z(469) + 2*rad*z(430)*z(53) - z(470) - 4*e3*z(61);
z(506) = rad*z(5)*(z(1)*z(41)+z(3)*z(39))*z(455) + z(14)*(z(1)*z(41)+z(3)*z(39))*z(462) + z(20)*(z(1)*z(41)-z(2)*z(42))*z(464) + 2*rad*e4*(z(4)*z(41)+z(6)*z(39))*z(56) + rad*z(6)*(z(1)*z(41)-z(2)*z(42))*z(493) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(500) + 2*rad*e2*(z(1)*z(41)-z(2)*z(42))*z(59) + 2*rad*z(429)*(z(1)*z(41)-z(2)*z(42))*z(52) + 2*e3*z(14)*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) + e2*z(20)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*z(14)*(4*e1*z(39)-8*e3*z(41)-z(1)*z(432)-z(3)*z(433))*z(70) + 0.5*z(20)*(4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434))*z(52) + rad*z(6)*z(429)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*rad*z(5)*(4*e1*z(39)-8*e3*z(41)-z(1)*z(432)-z(3)*z(433))*z(56) + 0.5*rad*z(6)*(4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434))*z(59) + 0.5*z(440)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433))*z(84) + 0.25*z(5)*z(14)*(z(3)*z(503)+4*e1*z(69)+4*z(39)*z(473)+4*z(41)*z(469)+8*e3*z(64)-2*z(432)*z(53)-2*z(433)*z(68)-z(1)*z(502)) + 0.25*z(6)*z(20)*(4*e4*z(66)+4*z(41)*z(469)+8*e3*z(64)-4*z(42)*z(471)-2*z(432)*z(53)-2*z(434)*z(65)-z(1)*z(502)-z(2)*z(504)) + 0.25*z(32)*(z(3)*z(503)+4*e1*z(69)+4*e4*z(66)+4*z(39)*z(473)+4*z(41)*z(469)+8*e3*z(64)-4*z(42)*z(471)-2*z(432)*z(53)-2*z(433)*z(68)-2*z(434)*z(65)-z(1)*z(502)-z(2)*z(504)) - rad*z(2)*(z(4)*z(41)+z(6)*z(39))*z(455) - z(14)*(z(4)*z(41)+z(6)*z(39))*z(471) - 4*rad*e3*(z(1)*z(41)+z(3)*z(39))*z(56) - z(20)*(z(4)*z(41)-z(5)*z(42))*z(473) - 2*rad*e1*(z(4)*z(41)-z(5)*z(42))*z(59) - 2*rad*z(429)*(z(4)*z(41)-z(5)*z(42))*z(68) - rad*z(3)*(z(4)*z(41)-z(5)*z(42))*z(493) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(505) - e4*z(14)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e1*z(20)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(14)*(4*e2*z(39)+4*e4*z(41)-z(4)*z(432)-z(6)*z(433))*z(65) - 0.5*z(20)*(4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434))*z(68) - rad*z(3)*z(429)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*rad*z(2)*(4*e2*z(39)+4*e4*z(41)-z(4)*z(432)-z(6)*z(433))*z(56) - 0.5*rad*z(3)*(4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434))*z(59) - 0.5*z(439)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433))*z(82) - 0.25*z(2)*z(14)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+4*z(41)*z(460)-4*e4*z(64)-2*z(432)*z(67)-2*z(433)*z(52)-z(4)*z(502)) - 0.25*z(3)*z(20)*(4*z(41)*z(460)+8*e3*z(66)-4*e4*z(64)-4*z(42)*z(462)-2*z(432)*z(67)-2*z(434)*z(70)-z(4)*z(502)-z(5)*z(504)) - 0.25*z(29)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+4*z(41)*z(460)+8*e3*z(66)-4*e4*z(64)-4*z(42)*z(462)-2*z(432)*z(67)-2*z(433)*z(52)-2*z(434)*z(70)-z(4)*z(502)-z(5)*z(504));
z(507) = z(463) + 4*e3*z(61) - z(5)*z(499) - z(21)*z(462) - 2*rad*z(430)*z(70);
z(508) = z(2)*z(499) + z(21)*z(471) + 2*rad*z(430)*z(65) - z(472) - 2*e4*z(61);
z(509) = rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(455) + z(22)*(z(1)*z(41)-z(2)*z(42))*z(464) + z(14)*(z(2)*z(42)-z(3)*z(39))*z(460) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(507) + 2*rad*e4*(z(2)*z(42)-z(3)*z(39))*z(56) + 4*rad*e3*(z(5)*z(42)-z(6)*z(39))*z(56) + z(22)*e2*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*z(22)*(4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434))*z(52) + 0.5*z(14)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433))*z(53) + 0.5*rad*z(1)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433))*z(56) + 0.5*z(444)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433))*z(74) + 0.25*z(1)*z(14)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+8*e3*z(66)-4*z(42)*z(462)-2*z(433)*z(52)-2*z(434)*z(70)-z(5)*z(504)) + 0.25*z(22)*z(6)*(4*e4*z(66)+4*z(41)*z(469)+8*e3*z(64)-4*z(42)*z(471)-2*z(432)*z(53)-2*z(434)*z(65)-z(1)*z(502)-z(2)*z(504)) + 0.25*z(33)*(z(3)*z(503)+4*e1*z(69)+4*e4*z(66)+4*z(39)*z(473)+4*z(41)*z(469)+8*e3*z(64)-4*z(42)*z(471)-2*z(432)*z(53)-2*z(433)*z(68)-2*z(434)*z(65)-z(1)*z(502)-z(2)*z(504)) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(455) - z(22)*(z(4)*z(41)-z(5)*z(42))*z(473) - z(14)*(z(5)*z(42)-z(6)*z(39))*z(469) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(508) - 2*e3*z(14)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - e4*z(14)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - z(22)*e1*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(22)*(4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434))*z(68) - 0.5*z(14)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433))*z(67) - 0.5*rad*z(4)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433))*z(56) - 0.5*z(443)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433))*z(72) - 0.25*z(4)*z(14)*(z(3)*z(503)+4*e1*z(69)+4*e4*z(66)+4*z(39)*z(473)-4*z(42)*z(471)-2*z(433)*z(68)-2*z(434)*z(65)-z(2)*z(504)) - 0.25*z(22)*z(3)*(4*z(41)*z(460)+8*e3*z(66)-4*e4*z(64)-4*z(42)*z(462)-2*z(432)*z(67)-2*z(434)*z(70)-z(4)*z(502)-z(5)*z(504)) - 0.25*z(30)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+4*z(41)*z(460)+8*e3*z(66)-4*e4*z(64)-4*z(42)*z(462)-2*z(432)*z(67)-2*z(433)*z(52)-2*z(434)*z(70)-z(4)*z(502)-z(5)*z(504));
z(510) = z(465) - 2*e2*z(61) - z(6)*z(499) - z(21)*z(464) - 2*rad*z(430)*z(52);
z(511) = z(3)*z(499) + z(21)*z(473) + 2*e1*z(61) + 2*rad*z(430)*z(68) - z(474);
z(512) = z(22)*(z(4)*z(41)+z(6)*z(39))*z(471) + z(20)*(z(2)*z(42)-z(3)*z(39))*z(460) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(493) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(510) + 2*rad*e4*(z(2)*z(42)-z(3)*z(39))*z(59) + 2*rad*z(429)*(z(2)*z(42)-z(3)*z(39))*z(67) + 4*rad*e3*(z(5)*z(42)-z(6)*z(39))*z(59) + z(22)*e4*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) + rad*z(1)*z(429)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(22)*(4*e2*z(39)+4*e4*z(41)-z(4)*z(432)-z(6)*z(433))*z(65) + 0.5*z(20)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433))*z(53) + 0.5*rad*z(1)*(4*e2*z(39)+8*e3*z(42)-z(5)*z(434)-z(6)*z(433))*z(59) + 0.5*z(436)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e1*z(39)+4*e4*z(42)-8*e3*z(41)-z(1)*z(432)-z(2)*z(434)-z(3)*z(433))*z(79) + 0.25*z(22)*z(2)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+4*z(41)*z(460)-4*e4*z(64)-2*z(432)*z(67)-2*z(433)*z(52)-z(4)*z(502)) + 0.25*z(1)*z(20)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+8*e3*z(66)-4*z(42)*z(462)-2*z(433)*z(52)-2*z(434)*z(70)-z(5)*z(504)) + 0.25*z(34)*(z(3)*z(503)+4*e1*z(69)+4*e4*z(66)+4*z(39)*z(473)+4*z(41)*z(469)+8*e3*z(64)-4*z(42)*z(471)-2*z(432)*z(53)-2*z(433)*z(68)-2*z(434)*z(65)-z(1)*z(502)-z(2)*z(504)) - z(22)*(z(1)*z(41)+z(3)*z(39))*z(462) - z(20)*(z(5)*z(42)-z(6)*z(39))*z(469) - 2*rad*z(429)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(493) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(511) - 2*z(22)*e3*(z(1)*z(64)-2*z(39)*z(68)-2*z(41)*z(53)-z(3)*z(69)) - 2*e3*z(20)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) - e4*z(20)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - rad*z(4)*z(429)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(22)*(4*e1*z(39)-8*e3*z(41)-z(1)*z(432)-z(3)*z(433))*z(70) - 0.5*z(20)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433))*z(67) - 0.5*rad*z(4)*(4*e1*z(39)+4*e4*z(42)-z(2)*z(434)-z(3)*z(433))*z(59) - 0.5*z(435)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*(4*e2*z(39)+4*e4*z(41)+8*e3*z(42)-z(4)*z(432)-z(5)*z(434)-z(6)*z(433))*z(77) - 0.25*z(22)*z(5)*(z(3)*z(503)+4*e1*z(69)+4*z(39)*z(473)+4*z(41)*z(469)+8*e3*z(64)-2*z(432)*z(53)-2*z(433)*z(68)-z(1)*z(502)) - 0.25*z(4)*z(20)*(z(3)*z(503)+4*e1*z(69)+4*e4*z(66)+4*z(39)*z(473)-4*z(42)*z(471)-2*z(433)*z(68)-2*z(434)*z(65)-z(2)*z(504)) - 0.25*z(31)*(z(6)*z(503)+4*e2*z(69)+4*z(39)*z(464)+4*z(41)*z(460)+8*e3*z(66)-4*e4*z(64)-4*z(42)*z(462)-2*z(432)*z(67)-2*z(433)*z(52)-2*z(434)*z(70)-z(4)*z(502)-z(5)*z(504));
z(513) = u1*z(506) + u2*z(509) + u3*z(512) + z(448)*u1*z(80) - z(447)*u1*z(75);
z(514) = (z(446)*(z(44)*z(90)-z(48)*z(86))-z(49)*(z(44)*z(498)+z(437)*z(90)-z(48)*z(513)-z(442)*z(86)))/z(49)^2;
z(515) = m*(2*z(22)*e3*z(26)+2*e4*z(14)*z(25)+4*e3*z(15)*z(28)+z(22)*e4*z(2)*z(15)+2*z(22)*e3*z(5)*z(15)+2*e1*z(5)*z(14)*z(15)-4*e3*z(14)*z(28)-2*e1*z(15)*z(26)-2*e2*z(15)*z(23)-2*e4*z(15)*z(25)-z(22)*e4*z(23)-2*e2*z(2)*z(14)*z(15));
z(516) = 8*m*z(15)*(e1*z(28)+e2*z(25)) + 4*z(134)*(e4*z(25)+e1*z(5)*z(15)-2*e3*z(28)-e2*z(2)*z(15));
z(517) = z(132)*z(447) + z(135)*z(448) + 2*z(50)*z(489) - 2*z(515) - z(51)*z(516);
z(518) = (z(446)*(z(43)*z(90)-z(46)*z(86))+z(49)*(z(46)*z(513)-z(43)*z(498)-z(445)*z(90)-2*e1*z(14)*z(86)))/z(49)^2;
z(519) = z(486) + z(139)*z(518) + z(490)*z(91) - z(138)*z(514) - z(517)*z(92);
z(520) = m*(4*e1*z(95)+z(22)*z(8)*z(448)+z(22)*z(9)*z(447)-2*e2*z(14)-2*e2*z(96)-2*z(7)*z(14)*z(447));
z(521) = m*(e4*z(14)*z(23)+2*e3*z(15)*z(26)-2*e3*z(14)*z(26)-e4*z(15)*z(23)-2*e3*z(5)*z(14)*z(15)-e4*z(2)*z(14)*z(15));
z(522) = z(128)*z(448) + 2*z(51)*z(515) - 8*z(521) - 2*z(50)*z(487) - z(127)*z(447);
z(523) = z(522) + z(50)*z(490) + z(139)*z(448) - z(51)*z(517) - z(138)*z(447);
z(524) = (z(523)*z(143)+z(141)*(9.810000000000001*z(520)-z(519)))/z(141)^2;
a14 = z(524);
z(525) = e1*(z(17)+z(8)^2/z(16)^0.5)/z(17)^2;
z(526) = e1*z(8)/(z(16)^0.5*z(17)^2);
z(527) = 4*rad*z(20)*z(525) + 4*z(22)*e2*z(35) + 4*z(213)*z(7)*z(526) + 4*e1*z(20)*z(35) + 4*rad*z(8)*z(35)*z(525) + 4*rad*z(9)*z(14)*z(526) - 4*rad*z(35)*z(526) - 4*rad*z(8)*z(20)*z(526);
z(528) = z(22)*z(527)/(z(37)^0.5*z(38)^2);
z(529) = z(14)*z(527)/(z(37)^0.5*z(38)^2);
z(530) = (4*rad*z(38)*z(525)-z(20)*z(527)/z(37)^0.5)/z(38)^2;
z(531) = 2*e1*z(21) - 2*e1*z(15) - 2*rad*z(6)*z(526);
z(532) = 2*e2*z(21) + 2*rad*z(3)*z(526) - 2*e2*z(15);
z(533) = 2*e3*z(20)*(z(2)*z(42)-z(3)*z(39)) + 2*rad*z(4)*z(525)*(z(2)*z(42)-z(3)*z(39)) + z(531)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(22)*z(2)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(6)*z(529)) + 0.5*z(34)*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529)) - 2*z(22)*e3*(z(4)*z(41)+z(6)*z(39)) - 2*rad*z(1)*z(525)*(z(5)*z(42)-z(6)*z(39)) - z(532)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(1)*z(20)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39)) - 0.5*z(22)*z(5)*(4*e2*z(39)-z(1)*z(528)-z(3)*z(529)) - 0.5*z(4)*z(20)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529)) - 0.5*z(31)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529));
z(534) = 2*e1*z(14) + 2*rad*z(9)*z(525);
z(535) = 2*e3*z(15) - 2*e3*z(21) - 2*rad*z(4)*z(526);
z(536) = 2*e3*z(14)*(z(4)*z(41)+z(6)*z(39)) + 2*rad*z(6)*z(525)*(z(1)*z(41)-z(2)*z(42)) + z(535)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) + 0.5*z(5)*z(14)*(4*e2*z(39)-z(1)*z(528)-z(3)*z(529)) + 0.5*z(6)*z(20)*(4*e3*z(42)-z(1)*z(528)-z(2)*z(530)) + 0.5*z(32)*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529)) - 2*e1*z(20)*(z(1)*z(41)-z(2)*z(42)) - 2*e2*z(20)*(z(4)*z(41)-z(5)*z(42)) - 2*rad*z(3)*z(525)*(z(4)*z(41)-z(5)*z(42)) - 2*rad*z(1)*z(526)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(3)*z(20)*(4*e3*z(41)-z(4)*z(528)-z(5)*z(530)) - 0.5*z(2)*z(14)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(6)*z(529)) - 0.5*z(29)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529));
z(537) = 2*e2*z(20) - 2*z(22)*e1 - 2*rad*z(7)*z(525);
z(538) = 2*e3*z(15) + 2*rad*z(2)*z(526) - 2*e3*z(21);
z(539) = 2*e3*z(14)*(z(2)*z(42)-z(3)*z(39)) + 0.5*z(22)*z(6)*(4*e3*z(42)-z(1)*z(528)-z(2)*z(530)) + 0.5*z(33)*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529)) - 2*z(22)*e1*(z(1)*z(41)-z(2)*z(42)) - 2*z(22)*e2*(z(4)*z(41)-z(5)*z(42)) - z(538)*(z(2)*z(42)-z(1)*z(41)-z(3)*z(39)) - 0.5*z(1)*z(14)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39)) - 2*rad*z(5)*z(526)*(z(5)*z(42)-z(4)*z(41)-z(6)*z(39)) - 0.5*z(22)*z(3)*(4*e3*z(41)-z(4)*z(528)-z(5)*z(530)) - 0.5*z(4)*z(14)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529)) - 0.5*z(30)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529));
z(540) = z(43)*z(537) + z(48)*z(539) - z(46)*z(533) - 2*e2*z(14)*z(44);
z(541) = (z(540)*(z(44)*z(47)-z(45)*z(48))-z(49)*(z(44)*z(534)+z(47)*z(533)-z(45)*z(537)-z(48)*z(536)))/z(49)^2;
z(542) = (z(540)*(z(43)*z(47)-z(45)*z(46))-z(49)*(z(43)*z(534)+z(47)*z(539)-z(46)*z(536)-2*e2*z(14)*z(45)))/z(49)^2;
z(543) = u1*(i22*z(541)*u3+i33*z(542)*u2-z(541)*z(101)-z(542)*z(100));
z(544) = u1*(i22*z(542)*u2+i33*z(541)*u3-z(541)*z(106)-z(542)*z(105));
z(545) = z(541)*u1*(z(14)*u2+z(110));
z(546) = (1+e2*z(542)+e3*z(541))*u1;
z(547) = 0.5*u2 - 0.5*e1*z(542)*u1 - 0.5*e4*z(541)*u1;
z(548) = -4*e2*z(547) - 2*e1*z(546);
z(549) = (z(11)+z(9)^2/z(10)^0.5)*z(548)/z(11)^2;
z(550) = u1*(z(14)*z(541)*u3+rad*z(549)-z(542)*z(110));
z(551) = rad*z(9)*z(548)/(z(10)^0.5*z(11)^2);
z(552) = (e2*z(541)-e3*z(542))*u1;
z(553) = 0.5*u3 + 0.5*e4*z(542)*u1 - 0.5*e1*z(541)*u1;
z(554) = e2*z(546) + e3*z(552) + 2*e1*z(547) + 2*e4*z(553) + 2*e3p;
z(555) = z(4)*z(551) + z(15)*z(554) + 2*e3*z(60);
z(556) = -4*e3*z(553) - 2*e1*z(546);
z(557) = z(5)*z(551) + z(15)*z(556);
z(558) = 2*e2*z(553) + 2*e3*z(547) - e1*z(552) - e4*z(546) - 2*e1p;
z(559) = z(6)*z(551) + z(15)*z(558) - 2*e1*z(60);
z(560) = u2*z(557) + u1*z(555) + u3*z(559) + z(542)*u1*z(76) - z(541)*u1*z(71);
z(561) = rad*(u2*z(549)-z(541)*u1*z(56));
z(562) = -z(542)*u1*z(111) - z(561);
z(563) = -4*e2*z(547) - 4*e3*z(553);
z(564) = z(1)*z(551) + z(15)*z(563);
z(565) = e2*z(546) + 2*e1*z(547) - 2*e4*z(553) - e3*z(552) - 2*e3p;
z(566) = z(2)*z(551) + z(15)*z(565) - 2*e3*z(60);
z(567) = e2*z(552) + e3*z(546) + 2*e1*z(553) + 2*e4*z(547) + 2*e2p;
z(568) = z(3)*z(551) + z(15)*z(567) + 2*e2*z(60);
z(569) = z(541)*u1*z(73) - u1*z(564) - u2*z(566) - u3*z(568) - z(542)*u1*z(78);
z(570) = u1*(z(22)*z(542)*u1+z(14)*z(541)*u2+z(541)*z(110));
z(571) = u1*(z(22)*z(541)*u1+z(14)*z(541)*u3+rad*z(549)-z(542)*z(110));
z(572) = z(22)*z(542)*u1*u3 - z(22)*z(541)*u1*u2 - z(541)*u1*z(122) - z(542)*u1*z(121) - z(561);
z(573) = z(543) + z(544) - m*(z(3)*z(23)*z(545)+z(5)*z(26)*z(550)+z(4)*z(26)*z(562)+2*e3*z(23)*z(120)+2*e3*z(26)*z(118)-z(14)*z(550)-2*e1*z(26)*z(119)-2*e2*z(23)*z(119)-z(2)*z(23)*z(550)-z(6)*z(26)*z(545)-2*e3*z(3)*z(15)*z(119)-2*z(23)*z(560)-4*e3*z(15)*z(116)-2*z(5)*z(14)*z(569)-z(1)*z(23)*z(562)-2*e3*z(1)*z(15)*z(118)-2*e3*z(2)*z(15)*z(120)) - m*(z(3)*z(23)*z(570)+z(5)*z(26)*z(571)+2*z(26)*z(569)+z(4)*z(26)*z(572)+2*e3*z(23)*z(125)+2*e3*z(26)*z(123)+4*e3*z(14)*z(116)-z(14)*z(571)-2*e1*z(26)*z(124)-2*e2*z(23)*z(124)-z(2)*z(23)*z(571)-z(6)*z(26)*z(570)-2*e3*z(3)*z(15)*z(124)-2*z(2)*z(14)*z(560)-z(1)*z(23)*z(572)-2*e3*z(1)*z(15)*z(123)-2*e3*z(2)*z(15)*z(125));
z(574) = z(542)*u1*(i33*u1-z(99));
z(575) = z(542)*u1*(i22*u1-z(99));
z(576) = -z(574) - z(575) - m*(z(3)*z(24)*z(545)+z(5)*z(27)*z(550)+z(14)*z(562)+2*z(27)*z(569)+z(4)*z(27)*z(562)+2*e3*z(24)*z(120)+2*e3*z(27)*z(118)+2*z(1)*z(14)*z(560)-2*e1*z(27)*z(119)-2*e2*z(24)*z(119)-z(2)*z(24)*z(550)-z(6)*z(27)*z(545)-2*e3*z(6)*z(15)*z(119)-2*z(24)*z(560)-4*e3*z(15)*z(117)-z(1)*z(24)*z(562)-2*e3*z(4)*z(15)*z(118)-2*e3*z(5)*z(15)*z(120)) - m*(z(22)*z(570)+z(3)*z(24)*z(570)+z(5)*z(27)*z(571)+z(14)*z(572)+z(4)*z(27)*z(572)+2*z(22)*e1*z(117)+2*e3*z(24)*z(125)+2*e3*z(27)*z(123)+2*z(4)*z(14)*z(569)+4*e3*z(14)*z(117)-2*e1*z(27)*z(124)-2*e2*z(24)*z(124)-z(2)*z(24)*z(571)-z(6)*z(27)*z(570)-2*e3*z(6)*z(15)*z(124)-2*z(22)*e2*z(116)-z(22)*z(3)*z(560)-z(22)*z(6)*z(569)-z(1)*z(24)*z(572)-2*e3*z(4)*z(15)*z(123)-2*e3*z(5)*z(15)*z(125));
z(577) = z(541)*u1*(i22*u1-z(99));
z(578) = z(541)*u1*(i33*u1-z(99));
z(579) = m*(z(2)*z(25)*z(550)+z(6)*z(28)*z(545)+2*e1*z(28)*z(119)+2*e2*z(25)*z(119)+2*z(25)*z(560)+z(1)*z(25)*z(562)-z(3)*z(25)*z(545)-z(5)*z(28)*z(550)-2*e1*z(3)*z(15)*z(119)-2*e2*z(6)*z(15)*z(119)-2*z(28)*z(569)-4*e1*z(15)*z(116)-4*e2*z(15)*z(117)-2*e3*z(25)*z(120)-2*e3*z(28)*z(118)-z(4)*z(28)*z(562)-2*e1*z(1)*z(15)*z(118)-2*e1*z(2)*z(15)*z(120)-2*e2*z(4)*z(15)*z(118)-2*e2*z(5)*z(15)*z(120)) + m*(z(2)*z(25)*z(571)+z(6)*z(28)*z(570)+2*e1*z(28)*z(124)+2*e2*z(25)*z(124)+z(1)*z(25)*z(572)+2*z(22)*e3*z(116)-z(22)*z(571)-z(3)*z(25)*z(570)-z(5)*z(28)*z(571)-2*e1*z(3)*z(15)*z(124)-2*e2*z(6)*z(15)*z(124)-2*e3*z(25)*z(125)-2*e3*z(28)*z(123)-z(22)*z(2)*z(560)-z(22)*z(5)*z(569)-z(4)*z(28)*z(572)-2*e1*z(1)*z(15)*z(123)-2*e1*z(2)*z(15)*z(125)-2*e2*z(4)*z(15)*z(123)-2*e2*z(5)*z(15)*z(125)) - z(577) - z(578);
z(580) = z(573) + z(50)*z(576) + z(542)*z(136) - z(51)*z(579) - z(541)*z(133);
z(581) = m*(2*e3*z(14)*z(24)+2*e3*z(15)*z(26)+2*e3*z(1)*z(14)*z(15)-2*e3*z(14)*z(26)-2*e3*z(15)*z(24)-z(22)*e1*z(26)-z(22)*e2*z(23)-2*e3*z(5)*z(14)*z(15)-z(22)*e3*z(3)*z(15));
z(582) = m*(z(22)*e1*z(27)+z(22)*e2*z(24)+2*e3*z(14)*z(27)+z(22)*e3*z(6)*z(15)-2*e3*z(15)*z(27)-2*e3*z(4)*z(14)*z(15));
z(583) = m*(2*e1*z(15)*z(24)+2*e3*z(15)*z(28)+z(22)*e1*z(3)*z(15)+z(22)*e2*z(6)*z(15)+z(22)*e3*z(5)*z(15)-2*e2*z(15)*z(27)-2*e3*z(14)*z(28)-z(22)*e1*z(28)-z(22)*e2*z(25)-z(22)*e3*z(24)-2*e1*z(1)*z(14)*z(15)-2*e2*z(4)*z(14)*z(15));
z(584) = 2*z(51)*z(583) + 4*z(50)*z(582) - 2*z(581) - z(131)*z(541) - z(132)*z(542);
z(585) = e1*z(552) + e4*z(546) + 2*e2*z(553) + 2*e3*z(547) + 2*e1p;
z(586) = (z(8)*z(16)*z(585)+2*e1*z(16)*z(57)+2*e1*z(8)^2*z(57))/z(16)^1.5;
z(587) = (z(17)*(z(8)*z(586)+z(17)*z(585)+2*e1*z(58))+2*e1*z(8)*(z(17)*z(57)+2*z(8)*z(58))/z(16)^0.5)/z(17)^3;
z(588) = z(20)*z(548) + rad*z(8)*z(549) + z(14)*z(585) + rad*z(9)*z(587) + 2*rad*e1*z(56) + 2*rad*z(525)*z(54);
z(589) = e3*z(546) + 2*e1*z(553) - 2*e4*z(547) - e2*z(552) - 2*e2p;
z(590) = z(22)*z(548) + 2*rad*e2*z(56) - rad*z(7)*z(549) - z(14)*z(589);
z(591) = 2*rad*e2*z(59) - z(22)*z(585) - z(20)*z(589) - 2*rad*z(525)*z(62) - rad*z(7)*z(587);
z(592) = u1*z(588) + u2*z(590) + u3*z(591) + z(542)*u1*z(89) - z(541)*u1*z(87);
z(593) = rad*(z(17)*z(586)+4*e1*z(8)*z(58)/z(16)^0.5)/z(17)^3;
z(594) = z(555) - 2*e3*z(61) - z(4)*z(593) - z(21)*z(554) - 2*rad*z(526)*z(67);
z(595) = 2*rad*z(14)*z(549) + 2*rad*z(20)*z(587) + 2*z(20)*z(35)*z(585) + 4*z(213)*z(526)*z(62) + 4*z(293)*z(525)*z(59) + 2*rad*z(8)*z(35)*z(587) + 4*rad*e1*z(35)*z(59) + 4*rad*z(14)*z(526)*z(54) + 4*rad*z(35)*z(525)*z(57) + 4*z(293)*z(9)*z(526)*z(56) + 2*z(35)*(z(551)-z(593)) + 2*z(8)*z(20)*(z(551)-z(593)) + 4*z(22)*e2*(z(60)-z(61)) + 4*e1*z(20)*(z(60)-z(61)) + 4*rad*z(8)*z(525)*(z(60)-z(61)) - 2*z(14)*z(35)*z(548) - 2*rad*z(9)*z(35)*z(549) - 2*z(22)*z(35)*z(589) - 4*rad*z(20)*z(526)*z(57) - 4*z(293)*z(8)*z(526)*z(59) - 4*rad*z(526)*(z(60)-z(61)) - 2*z(22)*z(7)*(z(551)-z(593)) - 2*z(9)*z(14)*(z(551)-z(593));
z(596) = z(22)*(2*z(595)/(z(37)^0.5*z(38)^2)-z(527)*(2+z(38)/z(37)^0.5)*z(63)/(z(37)*z(38)^3));
z(597) = (2*z(14)*z(527)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(549)+(z(14)*z(527)*z(63)-2*z(14)*z(37)*z(595)-2*rad*z(37)*z(527)*z(56))/z(37)^1.5))/z(38)^3;
z(598) = (2*z(20)*z(527)*z(63)/z(37)+z(38)*(4*rad*z(38)*z(587)+(z(20)*z(527)*z(63)-2*z(20)*z(37)*z(595)-4*rad*z(37)*z(525)*z(63)-2*rad*z(37)*z(527)*z(59))/z(37)^1.5))/z(38)^3;
z(599) = z(21)*z(563) + z(1)*z(593) + 2*rad*z(526)*z(53) - z(564);
z(600) = z(14)*(z(1)*z(41)+z(3)*z(39))*z(556) + rad*z(5)*(z(1)*z(41)+z(3)*z(39))*z(549) + z(20)*(z(1)*z(41)-z(2)*z(42))*z(558) + 2*rad*e3*(z(4)*z(41)+z(6)*z(39))*z(56) + rad*z(6)*(z(1)*z(41)-z(2)*z(42))*z(587) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(594) + 2*rad*z(525)*(z(1)*z(41)-z(2)*z(42))*z(52) + 0.5*z(14)*(4*e2*z(39)-z(1)*z(528)-z(3)*z(529))*z(70) + 0.5*z(20)*(4*e3*z(42)-z(1)*z(528)-z(2)*z(530))*z(52) + 0.5*rad*z(5)*(4*e2*z(39)-z(1)*z(528)-z(3)*z(529))*z(56) + 0.5*rad*z(6)*(4*e3*z(42)-z(1)*z(528)-z(2)*z(530))*z(59) + rad*z(6)*z(525)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.5*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529))*z(84) + rad*z(1)*z(526)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.25*z(5)*z(14)*(4*z(41)*z(563)+z(3)*z(597)+4*e2*z(69)+4*z(39)*z(567)-2*z(528)*z(53)-2*z(529)*z(68)-z(1)*z(596)) + 0.25*z(6)*z(20)*(4*z(41)*z(563)+4*e3*z(66)-4*z(42)*z(565)-2*z(528)*z(53)-2*z(530)*z(65)-z(1)*z(596)-z(2)*z(598)) + 0.25*z(32)*(4*z(41)*z(563)+z(3)*z(597)+4*e2*z(69)+4*e3*z(66)+4*z(39)*z(567)-4*z(42)*z(565)-2*z(528)*z(53)-2*z(529)*z(68)-2*z(530)*z(65)-z(1)*z(596)-z(2)*z(598)) - rad*z(2)*(z(4)*z(41)+z(6)*z(39))*z(549) - z(14)*(z(4)*z(41)+z(6)*z(39))*z(565) - z(20)*(z(4)*z(41)-z(5)*z(42))*z(567) - 2*rad*e1*(z(1)*z(41)-z(2)*z(42))*z(59) - 2*rad*e2*(z(4)*z(41)-z(5)*z(42))*z(59) - 2*rad*z(525)*(z(4)*z(41)-z(5)*z(42))*z(68) - rad*z(3)*(z(4)*z(41)-z(5)*z(42))*z(587) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(599) - 0.5*z(20)*(4*e3*z(41)-z(4)*z(528)-z(5)*z(530))*z(68) - 0.5*rad*z(3)*(4*e3*z(41)-z(4)*z(528)-z(5)*z(530))*z(59) - e3*z(14)*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) - e1*z(20)*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - e2*z(20)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(14)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(6)*z(529))*z(65) - rad*z(3)*z(525)*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*rad*z(2)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(6)*z(529))*z(56) - 0.5*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529))*z(82) - 0.5*z(535)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.25*z(3)*z(20)*(4*z(41)*z(554)-4*z(42)*z(556)-4*e3*z(64)-2*z(528)*z(67)-2*z(530)*z(70)-z(4)*z(596)-z(5)*z(598)) - 0.25*z(2)*z(14)*(z(6)*z(597)+4*z(39)*z(558)+4*z(41)*z(554)-4*e1*z(69)-4*e3*z(64)-2*z(528)*z(67)-2*z(529)*z(52)-z(4)*z(596)) - 0.25*z(29)*(z(6)*z(597)+4*z(39)*z(558)+4*z(41)*z(554)-4*z(42)*z(556)-4*e1*z(69)-4*e3*z(64)-2*z(528)*z(67)-2*z(529)*z(52)-2*z(530)*z(70)-z(4)*z(596)-z(5)*z(598));
z(601) = z(557) - z(21)*z(556) - z(5)*z(593) - 2*rad*z(526)*z(70);
z(602) = z(2)*z(593) + z(21)*z(565) + 2*rad*z(526)*z(65) - z(566) - 2*e3*z(61);
z(603) = rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(549) + z(22)*(z(1)*z(41)-z(2)*z(42))*z(558) + z(14)*(z(2)*z(42)-z(3)*z(39))*z(554) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(601) + 2*rad*e3*(z(2)*z(42)-z(3)*z(39))*z(56) + 0.5*z(22)*(4*e3*z(42)-z(1)*z(528)-z(2)*z(530))*z(52) + 0.5*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529))*z(74) + 0.5*z(538)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + rad*z(5)*z(526)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) + 0.25*z(1)*z(14)*(z(6)*z(597)+4*z(39)*z(558)-4*z(42)*z(556)-4*e1*z(69)-2*z(529)*z(52)-2*z(530)*z(70)-z(5)*z(598)) + 0.25*z(22)*z(6)*(4*z(41)*z(563)+4*e3*z(66)-4*z(42)*z(565)-2*z(528)*z(53)-2*z(530)*z(65)-z(1)*z(596)-z(2)*z(598)) + 0.25*z(33)*(4*z(41)*z(563)+z(3)*z(597)+4*e2*z(69)+4*e3*z(66)+4*z(39)*z(567)-4*z(42)*z(565)-2*z(528)*z(53)-2*z(529)*z(68)-2*z(530)*z(65)-z(1)*z(596)-z(2)*z(598)) - z(14)*(z(5)*z(42)-z(6)*z(39))*z(563) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(549) - z(22)*(z(4)*z(41)-z(5)*z(42))*z(567) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(602) - 0.5*z(14)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39))*z(53) - 0.5*rad*z(1)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39))*z(56) - 0.5*z(22)*(4*e3*z(41)-z(4)*z(528)-z(5)*z(530))*z(68) - e3*z(14)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - z(22)*e1*(2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) - z(22)*e2*(2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.5*z(14)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529))*z(67) - 0.5*rad*z(4)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529))*z(56) - 0.5*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529))*z(72) - 0.25*z(22)*z(3)*(4*z(41)*z(554)-4*z(42)*z(556)-4*e3*z(64)-2*z(528)*z(67)-2*z(530)*z(70)-z(4)*z(596)-z(5)*z(598)) - 0.25*z(4)*z(14)*(z(3)*z(597)+4*e2*z(69)+4*e3*z(66)+4*z(39)*z(567)-4*z(42)*z(565)-2*z(529)*z(68)-2*z(530)*z(65)-z(2)*z(598)) - 0.25*z(30)*(z(6)*z(597)+4*z(39)*z(558)+4*z(41)*z(554)-4*z(42)*z(556)-4*e1*z(69)-4*e3*z(64)-2*z(528)*z(67)-2*z(529)*z(52)-2*z(530)*z(70)-z(4)*z(596)-z(5)*z(598));
z(604) = z(559) + 2*e1*z(61) - z(6)*z(593) - z(21)*z(558) - 2*rad*z(526)*z(52);
z(605) = z(3)*z(593) + z(21)*z(567) + 2*e2*z(61) + 2*rad*z(526)*z(68) - z(568);
z(606) = z(22)*(z(4)*z(41)+z(6)*z(39))*z(565) + z(20)*(z(2)*z(42)-z(3)*z(39))*z(554) + rad*z(4)*(z(2)*z(42)-z(3)*z(39))*z(587) + (z(5)*z(42)-z(4)*z(41)-z(6)*z(39))*z(604) + 2*rad*e3*(z(2)*z(42)-z(3)*z(39))*z(59) + 2*rad*z(525)*(z(2)*z(42)-z(3)*z(39))*z(67) + z(22)*e3*(z(4)*z(64)-2*z(39)*z(52)-2*z(41)*z(67)-z(6)*z(69)) + rad*z(1)*z(525)*(z(6)*z(69)+2*z(39)*z(52)-2*z(42)*z(70)-z(5)*z(66)) + 0.5*z(22)*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(6)*z(529))*z(65) + 0.5*(4*e2*z(39)+4*e3*z(42)-z(1)*z(528)-z(2)*z(530)-z(3)*z(529))*z(79) + 0.5*z(532)*(z(3)*z(69)+2*z(39)*z(68)+2*z(41)*z(53)-2*z(42)*z(65)-z(1)*z(64)-z(2)*z(66)) + 0.25*z(1)*z(20)*(z(6)*z(597)+4*z(39)*z(558)-4*z(42)*z(556)-4*e1*z(69)-2*z(529)*z(52)-2*z(530)*z(70)-z(5)*z(598)) + 0.25*z(22)*z(2)*(z(6)*z(597)+4*z(39)*z(558)+4*z(41)*z(554)-4*e1*z(69)-4*e3*z(64)-2*z(528)*z(67)-2*z(529)*z(52)-z(4)*z(596)) + 0.25*z(34)*(4*z(41)*z(563)+z(3)*z(597)+4*e2*z(69)+4*e3*z(66)+4*z(39)*z(567)-4*z(42)*z(565)-2*z(528)*z(53)-2*z(529)*z(68)-2*z(530)*z(65)-z(1)*z(596)-z(2)*z(598)) - z(22)*(z(1)*z(41)+z(3)*z(39))*z(556) - z(20)*(z(5)*z(42)-z(6)*z(39))*z(563) - 2*rad*z(525)*(z(5)*z(42)-z(6)*z(39))*z(53) - rad*z(1)*(z(5)*z(42)-z(6)*z(39))*z(587) - (z(2)*z(42)-z(1)*z(41)-z(3)*z(39))*z(605) - 0.5*z(20)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39))*z(53) - 0.5*rad*z(1)*(z(5)*z(530)+z(6)*z(529)+4*e1*z(39))*z(59) - 0.5*z(22)*(4*e2*z(39)-z(1)*z(528)-z(3)*z(529))*z(70) - e3*z(20)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - rad*z(4)*z(525)*(z(3)*z(69)+2*z(39)*z(68)-2*z(42)*z(65)-z(2)*z(66)) - 0.5*z(20)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529))*z(67) - 0.5*rad*z(4)*(4*e2*z(39)+4*e3*z(42)-z(2)*z(530)-z(3)*z(529))*z(59) - 0.5*(4*e3*z(41)-4*e1*z(39)-z(4)*z(528)-z(5)*z(530)-z(6)*z(529))*z(77) - 0.5*z(531)*(z(6)*z(69)+2*z(39)*z(52)+2*z(41)*z(67)-2*z(42)*z(70)-z(4)*z(64)-z(5)*z(66)) - 0.25*z(22)*z(5)*(4*z(41)*z(563)+z(3)*z(597)+4*e2*z(69)+4*z(39)*z(567)-2*z(528)*z(53)-2*z(529)*z(68)-z(1)*z(596)) - 0.25*z(4)*z(20)*(z(3)*z(597)+4*e2*z(69)+4*e3*z(66)+4*z(39)*z(567)-4*z(42)*z(565)-2*z(529)*z(68)-2*z(530)*z(65)-z(2)*z(598)) - 0.25*z(31)*(z(6)*z(597)+4*z(39)*z(558)+4*z(41)*z(554)-4*z(42)*z(556)-4*e1*z(69)-4*e3*z(64)-2*z(528)*z(67)-2*z(529)*z(52)-2*z(530)*z(70)-z(4)*z(596)-z(5)*z(598));
z(607) = u1*z(600) + u2*z(603) + u3*z(606) + z(542)*u1*z(80) - z(541)*u1*z(75);
z(608) = (z(540)*(z(44)*z(90)-z(48)*z(86))-z(49)*(z(44)*z(592)+z(533)*z(90)-z(48)*z(607)-z(537)*z(86)))/z(49)^2;
z(609) = m*(z(22)*e3*z(23)+2*e2*z(15)*z(26)+2*e3*z(15)*z(25)-2*e1*z(15)*z(23)-2*e3*z(14)*z(25)-2*e1*z(2)*z(14)*z(15)-2*e2*z(5)*z(14)*z(15)-z(22)*e3*z(2)*z(15));
z(610) = 4*z(134)*(e3*z(25)+e1*z(2)*z(15)+e2*z(5)*z(15)) - 8*m*z(15)*(e1*z(25)-e2*z(28));
z(611) = 2*z(609) + z(132)*z(541) + z(135)*z(542) - 2*z(50)*z(583) - z(51)*z(610);
z(612) = (z(540)*(z(43)*z(90)-z(46)*z(86))-z(49)*(z(43)*z(592)+z(539)*z(90)-z(46)*z(607)-2*e2*z(14)*z(86)))/z(49)^2;
z(613) = z(580) + z(139)*z(612) + z(584)*z(91) - z(138)*z(608) - z(611)*z(92);
z(614) = m*(z(22)*z(8)*z(542)+z(22)*z(9)*z(541)-4*e2*z(95)-2*e1*z(14)-2*e1*z(96)-2*z(7)*z(14)*z(541));
z(615) = m*e3*(z(14)*z(23)-z(15)*z(23)-z(2)*z(14)*z(15));
z(616) = z(128)*z(542) - 8*z(615) - 2*z(50)*z(581) - 2*z(51)*z(609) - z(127)*z(541);
z(617) = z(616) + z(50)*z(584) + z(139)*z(542) - z(51)*z(611) - z(138)*z(541);
z(618) = (z(617)*z(143)+z(141)*(9.810000000000001*z(614)-z(613)))/z(141)^2;
a15 = z(618);
a21 = z(148);
a22 = 0.5*z(242);
a23 = z(352);
a24 = z(452);
a25 = 0.5*z(546);
a31 = z(149);
a32 = z(243);
a33 = -0.5*z(353);
a34 = 0.5*z(453);
a35 = z(547);
a41 = z(156);
a42 = z(250);
a43 = -0.5*z(360);
a44 = -0.5*z(459);
a45 = z(553);
a51 = z(155);
a52 = 0.5*z(249);
a53 = z(359);
a54 = z(458);
a55 = 0.5*z(552);
Amatr = [a11,a12,a13,a14,a15;a21,a22,a23,a24,a25;a31,a32,a33,a34,a35;a41,a42,a43,a44,a45;a51,a52,a53,a54,a55];
% eig(Amatr)
Output(1)=T; Output(2)=ke; Output(3)=pe; Output(4)=ke+pe; Output(5)=(-(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*z(14)*u2*z(4)-(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*z(14)*u2*z(1)+(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*(z(14)*u1-z(22)*u3)*z(5)+(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*(z(14)*u1-z(22)*u3)*z(2)+(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*z(22)*u2*z(6)+(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*z(22)*u2*z(3)-(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*z(20)*u3*z(4)-(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*z(20)*u3*z(1)+(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*z(20)*u1*z(6)+(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*z(20)*u1*z(3)+(z(5)*z(20)-z(22)*z(4)-z(6)*z(14))*(z(29)*u1+z(30)*u2+z(31)*u3)+(z(22)*z(1)+z(3)*z(14)-z(2)*z(20))*(z(32)*u1+z(33)*u2+z(34)*u3)); Output(6)=(-z(14)*u2*z(7)+(z(14)*u1-z(22)*u3)*z(8)+z(22)*u2*z(9)-z(20)*u3*z(7)+z(20)*u1*z(9)); Output(7)=e1; Output(8)=e2; Output(9)=e3; Output(10)=e4; Output(11)=u1; Output(12)=cahat1; Output(13)=cahat2; Output(14)=cahat3; Output(15)=cbhat1; Output(16)=cbhat2; Output(17)=cbhat3; Output(18)=u1; Output(19)=u2; Output(20)=u3;
FileIdentifier = fopen('all');
WriteOutput( 1, Output(1:20));
WriteOutput( FileIdentifier(1), Output(1:20));
%===========================================================================
function WriteOutput( fileIdentifier, Output )
numberOfOutputQuantities = length( Output );
if numberOfOutputQuantities > 0,
for i=1:numberOfOutputQuantities,
fprintf( fileIdentifier, ' %- 14.6E', Output(i) );
end
fprintf( fileIdentifier, '\n' );
end
%===========================================================================
% mdlTerminate: Perform end of simulation tasks and set sys=[]
%===========================================================================
function sys = mdlTerminate(T,VAR,u)
FileIdentifier = fopen('all');
fclose( FileIdentifier(1) );
fprintf( 1, '\n Output is in the file TwinDisks.1\n' );
fprintf( 1, '\n To load and plot columns 1 and 2 with a solid line and columns 1 and 3 with a dashed line, enter:\n' );
fprintf( 1, ' someName = load( ''TwinDisks.1'' );\n' );
fprintf( 1, ' plot( someName(:,1), someName(:,2), ''-'', someName(:,1), someName(:,3), ''--'' )\n\n' );
sys = [];
%===========================================================================
% Sfunction: System/Simulink function from standard template
%===========================================================================
function [sys,x0,str,ts] = Sfunction(t,x,u,flag)
switch flag,
case 0, [sys,x0,str,ts] = mdlInitializeSizes; % Initialization of sys, initial state x0, state ordering string str, and sample times ts
case 1, sys = mdlDerivatives(t,x,u); % Calculate the derivatives of continuous states and store them in sys
case 2, sys = mdlUpdate(t,x,u); % Update discrete states x(n+1) in sys
case 3, sys = mdlOutputs(t,x,u); % Calculate outputs in sys
case 4, sys = mdlGetTimeOfNextVarHit(t,x,u); % Return next sample time for variable-step in sys
case 9, sys = mdlTerminate(t,x,u); % Perform end of simulation tasks and set sys=[]
otherwise error(['Unhandled flag = ',num2str(flag)]);
end
%===========================================================================
% mdlInitializeSizes: Return the sizes, initial state VAR, and sample times ts
%===========================================================================
function [sys,VAR,stateOrderingStrings,timeSampling] = mdlInitializeSizes
sizes = simsizes; % Call simsizes to create a sizes structure
sizes.NumContStates = 17; % sys(1) is the number of continuous states
sizes.NumDiscStates = 0; % sys(2) is the number of discrete states
sizes.NumOutputs = 45; % sys(3) is the number of outputs
sizes.NumInputs = 0; % sys(4) is the number of inputs
sizes.DirFeedthrough = 1; % sys(6) is 1, and allows for the output to be a function of the input
sizes.NumSampleTimes = 1; % sys(7) is the number of samples times (the number of rows in ts)
sys = simsizes(sizes); % Convert it to a sizes array
stateOrderingStrings = [];
timeSampling = [0 0]; % m-by-2 matrix containing the sample times
OpenOutputFilesAndWriteHeadings
VAR = ReadUserInput