-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathTwinDisksEig.m
932 lines (891 loc) · 115 KB
/
TwinDisksEig.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
923
924
925
926
927
928
929
930
931
932
function TwinDisksEig(uu,tfina)
SolveOrdinaryDifferentialEquations(uu,tfina)
% File TwinDisksEig.m created by Autolev 4.1 on Sun Sep 26 14:52:12 2010
%===========================================================================
function VAR = ReadUserInput(uu,tfina)
global l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 q1 q2 q3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p q1p q2p q3p 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
%-------------------------------|--------------------------|-------------------|-----------------
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.1; % UNITS Initial Value
cbhat2 = 0.1414213562373095; % 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; % UNITS Initial Value
posa2 = 0.07071067811865475; % UNITS Initial Value
posa3 = 0.07071067811865477; % UNITS Initial Value
posb1 = -0.1; % UNITS Initial Value
posb2 = 0.07071067811865475; % UNITS Initial Value
posb3 = 0.07071067811865477; % UNITS Initial Value
q1 = 0.0; % UNITS Initial Value
q2 = 0.0; % UNITS Initial Value
q3 = 0.0; % UNITS Initial Value
u1 = uu; % UNITS Initial Value
TINITIAL = 0.0; % UNITS Initial Time
TFINAL = tfina; % 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(616,1);
% Evaluate constants
z(23) = l*rad;
z(1) = m*rad^2;
cahat3p = 0;
cbhat3p = 0;
z(37) = z(23)^2;
z(132) = m*z(23);
z(211) = rad*z(23);
z(291) = 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) = q1;
VAR(18) = q2;
VAR(19) = q3;
VAR(20) = u1;
%===========================================================================
function OpenOutputFilesAndWriteHeadings
FileIdentifier = fopen('TwinDisksEig.1', 'wt'); if( FileIdentifier == -1 ) error('Error: unable to open file TwinDisksEig.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 q1 q2 q3\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) (UNITS) (UNITS) (UNITS)\n\n' );
fprintf(FileIdentifier, '%% FILE: TwinDisksEig.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 q1 q2 q3\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) (UNITS) (UNITS) (UNITS)\n\n' );
%===========================================================================
% Main driver loop for numerical integration of differential equations
%===========================================================================
function SolveOrdinaryDifferentialEquations(uu,tfina)
global l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 q1 q2 q3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p q1p q2p q3p 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(uu,tfina);
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('TwinDisksEig.1');
% size(data)
% 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))
% figure(4), plot(data(:,1),data(:,46),data(:,1),data(:,47),data(:,1),data(:,48))
%===========================================================================
% mdlDerivatives: Calculates and returns the derivatives of the continuous states
%===========================================================================
function sys = mdlDerivatives(T,VAR,u)
global l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 q1 q2 q3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p q1p q2p q3p 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);
q1 = VAR(17);
q2 = VAR(18);
q3 = VAR(19);
u1 = VAR(20);
z(3) = 2*e1*e2 - 2*e3*e4;
z(5) = 2*e1*e2 + 2*e3*e4;
z(10) = 1 - 2*e1^2 - 2*e2^2;
z(11) = 1 - z(10)^2;
z(12) = z(11)^0.5;
z(13) = z(10)/z(12);
z(15) = rad*z(13);
z(9) = 2*e1*e4 + 2*e2*e3;
z(17) = 1 - z(9)^2;
z(18) = z(17)^0.5;
z(19) = z(9)/z(18);
z(21) = rad*z(19);
z(14) = 1/z(12);
z(16) = rad*z(14);
z(20) = 1/z(18);
z(22) = rad*z(20);
z(36) = z(16) - z(22);
z(8) = 2*e1*e3 - 2*e2*e4;
z(38) = z(37) + z(15)^2 + z(21)^2 + z(36)^2 + 2*z(9)*z(21)*z(36) - 2*z(23)*z(8)*z(36) - 2*z(10)*z(15)*z(36);
z(39) = z(38)^0.5;
z(42) = z(23)/z(39);
z(7) = 2*e2*e3 - 2*e1*e4;
z(40) = z(15)/z(39);
z(43) = z(21)/z(39);
z(4) = 2*e1*e3 + 2*e2*e4;
z(26) = z(7)*z(16);
z(32) = z(26) - z(7)*z(22);
z(6) = 1 - 2*e1^2 - 2*e3^2;
z(2) = 1 - 2*e2^2 - 2*e3^2;
z(29) = z(4)*z(16);
z(35) = z(4)*z(22) - z(29);
z(45) = z(23)*z(3)*(z(5)*z(42)+z(7)*z(40)) + z(5)*z(21)*(z(3)*z(43)-z(4)*z(40)) + z(32)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) - z(23)*z(6)*(z(2)*z(42)+z(4)*z(40)) - z(2)*z(21)*(z(6)*z(43)-z(7)*z(40)) - z(35)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40));
z(48) = z(9)*z(15) + z(10)*z(21);
z(24) = z(5)*z(16);
z(30) = z(24) - z(5)*z(22);
z(27) = z(2)*z(16);
z(33) = z(2)*z(22) - z(27);
z(46) = z(6)*z(15)*(z(2)*z(42)+z(4)*z(40)) + z(7)*z(21)*(z(2)*z(42)-z(3)*z(43)) + z(30)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) - z(3)*z(15)*(z(5)*z(42)+z(7)*z(40)) - z(4)*z(21)*(z(5)*z(42)-z(6)*z(43)) - z(33)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40));
z(49) = -z(23)*z(9) - z(8)*z(21);
z(25) = z(6)*z(16);
z(31) = z(25) - z(6)*z(22);
z(28) = z(3)*z(16);
z(34) = z(3)*z(22) - z(28);
z(44) = z(23)*z(7)*(z(2)*z(42)-z(3)*z(43)) + z(5)*z(15)*(z(3)*z(43)-z(4)*z(40)) + z(31)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) - z(23)*z(4)*(z(5)*z(42)-z(6)*z(43)) - z(2)*z(15)*(z(6)*z(43)-z(7)*z(40)) - z(34)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40));
z(47) = z(23)*z(10) - z(8)*z(15);
z(50) = z(44)*z(49) - z(45)*z(47);
z(51) = (z(45)*z(48)-z(46)*z(49))/z(50);
u2 = z(51)*u1;
z(52) = (z(44)*z(48)-z(46)*z(47))/z(50);
u3 = -z(52)*u1;
z(102) = z(1)*u3;
z(101) = z(1)*u2;
z(105) = 0.5*u2*z(102) - 0.25*u3*z(101);
z(108) = 0.25*u2*z(102) - 0.5*u3*z(101);
z(110) = z(15)*u1;
z(109) = z(15)*u2;
z(118) = u1*z(110) + u2*z(109);
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(55) = -4*e1*e1p - 4*e2*e2p;
z(56) = z(10)*z(55)/z(11)^0.5;
z(57) = (z(10)*z(56)+z(12)*z(55))/z(12)^2;
z(114) = rad*u1*z(57);
z(119) = z(114) - u3*z(109);
z(61) = rad*z(56)/z(12)^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(68) = 2*e1*e2p + 2*e2*e1p + 2*e3*e4p + 2*e4*e3p;
z(82) = z(5)*z(61) + z(16)*z(68);
z(71) = -4*e1*e1p - 4*e3*e3p;
z(72) = z(6)*z(61) + z(16)*z(71);
z(53) = 2*e2*e3p + 2*e3*e2p - 2*e1*e4p - 2*e4*e1p;
z(77) = z(7)*z(61) + z(16)*z(53);
z(115) = u1*z(82) + u2*z(72) + u3*z(77);
z(113) = rad*u2*z(57);
z(117) = -u3*z(110) - z(113);
z(54) = -4*e2*e2p - 4*e3*e3p;
z(84) = z(2)*z(61) + z(16)*z(54);
z(66) = 2*e1*e2p + 2*e2*e1p - 2*e3*e4p - 2*e4*e3p;
z(74) = z(3)*z(61) + z(16)*z(66);
z(69) = 2*e1*e3p + 2*e2*e4p + 2*e3*e1p + 2*e4*e2p;
z(79) = z(4)*z(61) + z(16)*z(69);
z(116) = -u1*z(84) - u2*z(74) - u3*z(79);
z(120) = z(15)*u1 - z(23)*u3;
z(123) = u1*z(120) + u2*z(109);
z(121) = z(23)*u2;
z(124) = z(114) - u1*z(121) - u3*z(109);
z(122) = u2*z(121) - u3*z(120) - z(113);
z(128) = z(105) + z(108) + m*(z(4)*z(24)*z(118)+z(15)*z(119)+2*z(24)*z(115)+z(2)*z(24)*z(117)+z(3)*z(24)*z(119)+2*z(6)*z(15)*z(116)-z(7)*z(27)*z(118)-z(5)*z(27)*z(117)-z(6)*z(27)*z(119)) + m*(z(4)*z(24)*z(123)+z(15)*z(124)+z(2)*z(24)*z(122)+z(3)*z(24)*z(124)+2*z(3)*z(15)*z(115)-z(7)*z(27)*z(123)-2*z(27)*z(116)-z(5)*z(27)*z(122)-z(6)*z(27)*z(124));
z(100) = z(1)*u1;
z(104) = 0.25*u3*z(100) - 0.5*u1*z(102);
z(107) = 0.25*u1*z(102) - 0.25*u3*z(100);
z(131) = z(104) + m*(z(4)*z(25)*z(118)+2*z(25)*z(115)+z(2)*z(25)*z(117)+z(3)*z(25)*z(119)-z(7)*z(28)*z(118)-2*z(28)*z(116)-z(15)*z(117)-2*z(2)*z(15)*z(115)-z(5)*z(28)*z(117)-z(6)*z(28)*z(119)) + m*(z(23)*z(123)+z(4)*z(25)*z(123)+z(23)*z(4)*z(115)+z(23)*z(7)*z(116)+z(2)*z(25)*z(122)+z(3)*z(25)*z(124)-z(7)*z(28)*z(123)-z(15)*z(122)-2*z(5)*z(15)*z(116)-z(5)*z(28)*z(122)-z(6)*z(28)*z(124)) - z(107);
z(103) = 0.25*u1*z(101) - 0.25*u2*z(100);
z(106) = 0.5*u1*z(101) - 0.25*u2*z(100);
z(134) = z(103) + z(106) + m*(z(4)*z(26)*z(118)+2*z(26)*z(115)+z(2)*z(26)*z(117)+z(3)*z(26)*z(119)-z(7)*z(29)*z(118)-2*z(29)*z(116)-z(5)*z(29)*z(117)-z(6)*z(29)*z(119)) + m*(z(4)*z(26)*z(123)+z(2)*z(26)*z(122)+z(3)*z(26)*z(124)-z(7)*z(29)*z(123)-z(23)*z(124)-z(23)*z(3)*z(115)-z(23)*z(6)*z(116)-z(5)*z(29)*z(122)-z(6)*z(29)*z(124));
z(138) = z(128) + z(51)*z(131) - z(52)*z(134);
z(126) = m*(2*z(24)*z(25)+2*z(27)*z(28)+z(23)*z(4)*z(24)+2*z(3)*z(15)*z(25)+2*z(5)*z(15)*z(27)-2*z(2)*z(15)*z(24)-2*z(6)*z(15)*z(28)-z(23)*z(7)*z(27));
z(129) = 0.75*z(1) + 2*m*(z(15)^2+z(25)^2+z(28)^2-2*z(2)*z(15)*z(25)) - m*(2*z(23)*z(7)*z(28)-z(37)-4*z(5)*z(15)*z(28)-2*z(23)*z(4)*z(25));
z(130) = m*(z(23)*z(3)*z(25)+z(23)*z(7)*z(29)+2*z(2)*z(15)*z(26)-2*z(25)*z(26)-2*z(28)*z(29)-2*z(5)*z(15)*z(29)-z(23)*z(4)*z(26)-z(23)*z(6)*z(28));
z(136) = z(126) + z(51)*z(129) + z(52)*z(130);
z(58) = 2*e1*e4p + 2*e2*e3p + 2*e3*e2p + 2*e4*e1p;
z(59) = z(9)*z(58)/z(17)^0.5;
z(60) = (z(9)*z(59)+z(18)*z(58))/z(18)^2;
z(89) = z(15)*z(58) + z(21)*z(55) + rad*z(9)*z(57) + rad*z(10)*z(60);
z(63) = 2*e1*e3p + 2*e3*e1p - 2*e2*e4p - 2*e4*e2p;
z(88) = z(23)*z(55) - z(15)*z(63) - rad*z(8)*z(57);
z(90) = -z(23)*z(58) - z(21)*z(63) - rad*z(8)*z(60);
z(91) = u1*z(89) + u2*z(88) + u3*z(90);
z(62) = rad*z(59)/z(18)^2;
z(83) = z(82) - z(5)*z(62) - z(22)*z(68);
z(64) = 2*rad*z(15)*z(57) + 2*rad*z(21)*z(60) + 2*z(21)*z(36)*z(58) + 2*rad*z(9)*z(36)*z(60) + 2*z(36)*(z(61)-z(62)) + 2*z(9)*z(21)*(z(61)-z(62)) - 2*z(23)*z(36)*z(63) - 2*z(15)*z(36)*z(55) - 2*rad*z(10)*z(36)*z(57) - 2*z(23)*z(8)*(z(61)-z(62)) - 2*z(10)*z(15)*(z(61)-z(62));
z(65) = z(23)*z(64)/(z(38)^0.5*z(39)^2);
z(70) = (2*rad*z(39)*z(57)-z(15)*z(64)/z(38)^0.5)/z(39)^2;
z(67) = (2*rad*z(39)*z(60)-z(21)*z(64)/z(38)^0.5)/z(39)^2;
z(85) = z(2)*z(62) + z(22)*z(54) - z(84);
z(86) = z(15)*(z(2)*z(42)+z(4)*z(40))*z(71) + rad*z(6)*(z(2)*z(42)+z(4)*z(40))*z(57) + z(21)*(z(2)*z(42)-z(3)*z(43))*z(53) + rad*z(7)*(z(2)*z(42)-z(3)*z(43))*z(60) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(83) + 0.5*z(3)*z(15)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) + 0.5*z(7)*z(21)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*z(33)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - z(15)*(z(5)*z(42)+z(7)*z(40))*z(66) - rad*z(3)*(z(5)*z(42)+z(7)*z(40))*z(57) - z(21)*(z(5)*z(42)-z(6)*z(43))*z(69) - rad*z(4)*(z(5)*z(42)-z(6)*z(43))*z(60) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(85) - 0.5*z(6)*z(15)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) - 0.5*z(4)*z(21)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(30)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67));
z(73) = z(72) - z(6)*z(62) - z(22)*z(71);
z(75) = z(3)*z(62) + z(22)*z(66) - z(74);
z(76) = z(23)*(z(2)*z(42)-z(3)*z(43))*z(53) + z(15)*(z(3)*z(43)-z(4)*z(40))*z(68) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(57) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(73) + 0.5*z(2)*z(15)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(23)*z(7)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*z(34)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - z(23)*(z(5)*z(42)-z(6)*z(43))*z(69) - z(15)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(57) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(75) - 0.5*z(5)*z(15)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(23)*z(4)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(31)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67));
z(78) = z(77) - z(7)*z(62) - z(22)*z(53);
z(80) = z(4)*z(62) + z(22)*z(69) - z(79);
z(81) = z(23)*(z(5)*z(42)+z(7)*z(40))*z(66) + z(21)*(z(3)*z(43)-z(4)*z(40))*z(68) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(60) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(78) + 0.5*z(23)*z(6)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) + 0.5*z(2)*z(21)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(35)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - z(23)*(z(2)*z(42)+z(4)*z(40))*z(71) - z(21)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(60) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(80) - 0.5*z(23)*z(3)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - 0.5*z(5)*z(21)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(32)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67));
z(87) = u1*z(86) + u2*z(76) + u3*z(81);
z(92) = (z(45)*z(91)-z(49)*z(87))/z(50);
z(127) = m*(2*z(24)*z(26)+2*z(27)*z(29)+z(23)*z(6)*z(27)+2*z(3)*z(15)*z(26)-z(23)*z(15)-2*z(6)*z(15)*z(29)-z(23)*z(3)*z(24));
z(133) = 0.75*z(1) + 2*m*(z(26)^2+z(29)^2) - z(132)*(2*z(3)*z(26)-z(23)-2*z(6)*z(29));
z(137) = z(127) - z(51)*z(130) - z(52)*z(133);
z(93) = (z(44)*z(91)-z(47)*z(87))/z(50);
z(140) = z(138) + z(136)*z(92) - z(137)*z(93);
z(96) = z(15)*z(51);
z(97) = z(15) + z(23)*z(52);
z(98) = z(23)*z(51);
z(99) = m*(2*z(8)*z(96)-z(9)*z(15)-z(9)*z(97)-z(10)*z(98));
z(141) = z(140) - 9.810000000000001*z(99);
z(125) = 0.5*z(1) + 2*m*(z(15)^2+z(24)^2+z(27)^2+2*z(3)*z(15)*z(24)-2*z(6)*z(15)*z(27));
z(135) = z(125) + z(51)*z(126) - z(52)*z(127);
z(139) = z(135) + z(51)*z(136) - z(52)*z(137);
z(142) = z(141)/z(139);
u1p = -z(142);
posa1p = z(24)*u1 + z(25)*u2 + z(26)*u3 + z(3)*z(15)*u1 - z(2)*z(15)*u2;
posa2p = z(6)*z(15)*u1 - z(27)*u1 - z(28)*u2 - z(29)*u3 - z(5)*z(15)*u2;
posa3p = -z(15)*(z(8)*u2-z(9)*u1);
posb1p = z(24)*u1 + z(25)*u2 + z(26)*u3 + z(23)*z(4)*u2 - z(2)*z(15)*u2 - z(3)*(z(23)*u3-z(15)*u1);
posb2p = z(23)*z(7)*u2 - z(27)*u1 - z(28)*u2 - z(29)*u3 - z(5)*z(15)*u2 - z(6)*(z(23)*u3-z(15)*u1);
posb3p = z(23)*z(10)*u2 - z(8)*z(15)*u2 - z(9)*(z(23)*u3-z(15)*u1);
cahat1p = z(24)*u1 + z(25)*u2 + z(26)*u3 - z(5)*z(16)*u1 - z(6)*z(16)*u2;
cahat2p = z(2)*z(16)*u1 + z(3)*z(16)*u2 - z(27)*u1 - z(28)*u2 - z(29)*u3;
cbhat1p = z(24)*u1 + z(25)*u2 + z(26)*u3 + z(23)*z(4)*u2 + z(4)*z(21)*u1 - z(2)*z(15)*u2 - z(2)*z(21)*u3 - z(5)*z(22)*u1 - z(7)*z(22)*u3 - z(3)*(z(23)*u3-z(15)*u1);
cbhat2p = z(23)*z(7)*u2 + z(2)*z(22)*u1 + z(4)*z(22)*u3 + z(7)*z(21)*u1 - z(27)*u1 - z(28)*u2 - z(29)*u3 - z(5)*z(15)*u2 - z(5)*z(21)*u3 - z(6)*(z(23)*u3-z(15)*u1);
q1p = u1;
q2p = u2;
q3p = u3;
% 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) = q1p;
VARp(18) = q2p;
VARp(19) = q3p;
VARp(20) = u1p;
sys = VARp';
%===========================================================================
% mdlOutputs: Calculates and return the outputs
%===========================================================================
function Output = mdlOutputs(T,VAR,u)
global l m rad;
global cahat1 cahat2 cahat3 cbhat1 cbhat2 cbhat3 e1 e2 e3 e4 posa1 posa2 posa3 posb1 posb2 posb3 q1 q2 q3 u1;
global ke pe u2 u3 cahat1p cahat2p cahat3p cbhat1p cbhat2p cbhat3p e1p e2p e3p e4p posa1p posa2p posa3p posb1p posb2p posb3p q1p q2p q3p 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 = 0.375*z(1)*u3^2 + 0.25*z(1)*(u1^2+1.5*u2^2) + 0.5*m*(z(23)^2*u2^2+(z(23)*u3-z(15)*u1)^2+2*z(23)*z(4)*u2*(z(24)*u1+z(25)*u2+z(26)*u3)+2*z(23)*z(6)*u3*(z(27)*u1+z(28)*u2+z(29)*u3)-2*z(23)*z(3)*u3*(z(24)*u1+z(25)*u2+z(26)*u3)-2*z(23)*z(7)*u2*(z(27)*u1+z(28)*u2+z(29)*u3)) + 0.5*m*(z(15)^2*u1^2+2*z(15)^2*u2^2+2*(z(24)*u1+z(25)*u2+z(26)*u3)^2+2*(z(27)*u1+z(28)*u2+z(29)*u3)^2+4*z(3)*z(15)*u1*(z(24)*u1+z(25)*u2+z(26)*u3)+4*z(5)*z(15)*u2*(z(27)*u1+z(28)*u2+z(29)*u3)-4*z(2)*z(15)*u2*(z(24)*u1+z(25)*u2+z(26)*u3)-4*z(6)*z(15)*u1*(z(27)*u1+z(28)*u2+z(29)*u3));
pe = 9.810000000000001*m*(posa3+posb3);
z(143) = 0.25*z(52)*z(101) + 0.5*z(51)*z(102) - 0.5*z(1)*z(52)*u2 - 0.25*z(1)*z(51)*u3;
z(144) = 0.25*z(51)*z(102) + 0.5*z(52)*z(101) - 0.5*z(1)*z(51)*u3 - 0.25*z(1)*z(52)*u2;
z(145) = z(15)*u1 + z(15)*z(51)*u2 + z(110) + z(51)*z(109);
z(146) = 0.5*e4 - 0.5*e2*z(52) - 0.5*e3*z(51);
z(147) = 0.5*e3 + 0.5*e1*z(52) + 0.5*e4*z(51);
z(148) = -4*e1*z(146) - 4*e2*z(147);
z(149) = z(148)*(z(12)+z(10)^2/z(11)^0.5)/z(12)^2;
z(150) = rad*(z(149)*u1+z(57));
z(151) = z(52)*z(109) + z(150) - z(15)*z(51)*u3;
z(152) = rad*z(10)*z(148)/(z(11)^0.5*z(12)^2);
z(153) = 0.5*e3*z(52) - 0.5*e1 - 0.5*e2*z(51);
z(154) = 0.5*e1*z(51) - 0.5*e2 - 0.5*e4*z(52);
z(155) = 2*e1*z(147) + 2*e2*z(146) + 2*e3*z(153) + 2*e4*z(154);
z(156) = z(5)*z(152) + z(16)*z(155);
z(157) = -4*e1*z(146) - 4*e3*z(154);
z(158) = z(6)*z(152) + z(16)*z(157);
z(159) = 2*e2*z(154) + 2*e3*z(147) - 2*e1*z(153) - 2*e4*z(146);
z(160) = z(7)*z(152) + z(16)*z(159);
z(161) = z(156)*u1 + z(158)*u2 + z(160)*u3 + z(82) + z(51)*z(72) - z(52)*z(77);
z(162) = rad*(z(149)*u2+z(51)*z(57));
z(163) = z(52)*z(110) - z(15)*u3 - z(162);
z(164) = -4*e2*z(147) - 4*e3*z(154);
z(165) = z(2)*z(152) + z(16)*z(164);
z(166) = 2*e1*z(147) + 2*e2*z(146) - 2*e3*z(153) - 2*e4*z(154);
z(167) = z(3)*z(152) + z(16)*z(166);
z(168) = 2*e1*z(154) + 2*e2*z(153) + 2*e3*z(146) + 2*e4*z(147);
z(169) = z(4)*z(152) + z(16)*z(168);
z(170) = z(52)*z(79) - z(165)*u1 - z(167)*u2 - z(169)*u3 - z(84) - z(51)*z(74);
z(171) = z(15)*z(51)*u2 + (z(15)+z(23)*z(52))*u1 + z(120) + z(51)*z(109);
z(172) = z(52)*z(109) + z(150) - z(23)*z(51)*u1 - z(15)*z(51)*u3 - z(121);
z(173) = z(23)*z(51)*u2 + z(51)*z(121) + z(52)*z(120) - (z(15)+z(23)*z(52))*u3 - z(162);
z(174) = z(143) + z(144) + m*(z(4)*z(24)*z(145)+z(15)*z(151)+2*z(24)*z(161)+z(2)*z(24)*z(163)+z(3)*z(24)*z(151)+2*z(6)*z(15)*z(170)-z(7)*z(27)*z(145)-z(5)*z(27)*z(163)-z(6)*z(27)*z(151)) + m*(z(4)*z(24)*z(171)+z(15)*z(172)+z(2)*z(24)*z(173)+z(3)*z(24)*z(172)+2*z(3)*z(15)*z(161)-z(7)*z(27)*z(171)-2*z(27)*z(170)-z(5)*z(27)*z(173)-z(6)*z(27)*z(172));
z(175) = 0.25*z(1)*u3 + 0.5*z(1)*z(52)*u1 - 0.5*z(102) - 0.25*z(52)*z(100);
z(176) = 0.25*z(102) + 0.25*z(52)*z(100) - 0.25*z(1)*u3 - 0.25*z(1)*z(52)*u1;
z(177) = z(175) + m*(z(4)*z(25)*z(145)+2*z(25)*z(161)+z(2)*z(25)*z(163)+z(3)*z(25)*z(151)-z(7)*z(28)*z(145)-2*z(28)*z(170)-z(15)*z(163)-2*z(2)*z(15)*z(161)-z(5)*z(28)*z(163)-z(6)*z(28)*z(151)) + m*(z(23)*z(171)+z(4)*z(25)*z(171)+z(23)*z(4)*z(161)+z(23)*z(7)*z(170)+z(2)*z(25)*z(173)+z(3)*z(25)*z(172)-z(7)*z(28)*z(171)-z(15)*z(173)-2*z(5)*z(15)*z(170)-z(5)*z(28)*z(173)-z(6)*z(28)*z(172)) - z(176);
z(178) = 0.25*z(1)*z(51)*u1 + 0.25*z(101) - 0.25*z(1)*u2 - 0.25*z(51)*z(100);
z(179) = 0.5*z(1)*z(51)*u1 + 0.5*z(101) - 0.25*z(1)*u2 - 0.25*z(51)*z(100);
z(180) = z(178) + z(179) + m*(z(4)*z(26)*z(145)+2*z(26)*z(161)+z(2)*z(26)*z(163)+z(3)*z(26)*z(151)-z(7)*z(29)*z(145)-2*z(29)*z(170)-z(5)*z(29)*z(163)-z(6)*z(29)*z(151)) + m*(z(4)*z(26)*z(171)+z(2)*z(26)*z(173)+z(3)*z(26)*z(172)-z(7)*z(29)*z(171)-z(23)*z(172)-z(23)*z(3)*z(161)-z(23)*z(6)*z(170)-z(5)*z(29)*z(173)-z(6)*z(29)*z(172));
z(181) = z(174) + z(51)*z(177) - z(52)*z(180);
z(182) = 2*e1*z(153) + 2*e2*z(154) + 2*e3*z(147) + 2*e4*z(146);
z(183) = z(182)*(z(18)+z(9)^2/z(17)^0.5)/z(18)^2;
z(184) = z(15)*z(182) + z(21)*z(148) + rad*z(9)*z(149) + rad*z(10)*z(183);
z(185) = 2*e1*z(154) + 2*e3*z(146) - 2*e2*z(153) - 2*e4*z(147);
z(186) = z(23)*z(148) - z(15)*z(185) - rad*z(8)*z(149);
z(187) = -z(23)*z(182) - z(21)*z(185) - rad*z(8)*z(183);
z(188) = z(184)*u1 + z(186)*u2 + z(187)*u3 + z(89) + z(51)*z(88) - z(52)*z(90);
z(189) = rad*z(9)*z(182)/(z(17)^0.5*z(18)^2);
z(190) = z(156) - z(5)*z(189) - z(22)*z(155);
z(191) = 2*rad*z(15)*z(149) + 2*rad*z(21)*z(183) + 2*z(21)*z(36)*z(182) + 2*rad*z(9)*z(36)*z(183) + 2*z(36)*(z(152)-z(189)) + 2*z(9)*z(21)*(z(152)-z(189)) - 2*z(23)*z(36)*z(185) - 2*z(15)*z(36)*z(148) - 2*rad*z(10)*z(36)*z(149) - 2*z(23)*z(8)*(z(152)-z(189)) - 2*z(10)*z(15)*(z(152)-z(189));
z(192) = z(23)*z(191)/(z(38)^0.5*z(39)^2);
z(193) = (2*rad*z(39)*z(149)-z(15)*z(191)/z(38)^0.5)/z(39)^2;
z(194) = (2*rad*z(39)*z(183)-z(21)*z(191)/z(38)^0.5)/z(39)^2;
z(195) = z(2)*z(189) + z(22)*z(164) - z(165);
z(196) = z(15)*z(157)*(z(2)*z(42)+z(4)*z(40)) + rad*z(6)*z(149)*(z(2)*z(42)+z(4)*z(40)) + z(21)*z(159)*(z(2)*z(42)-z(3)*z(43)) + rad*z(7)*z(183)*(z(2)*z(42)-z(3)*z(43)) + z(190)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(3)*z(15)*(z(5)*z(192)-2*z(40)*z(159)-2*z(42)*z(155)-z(7)*z(193)) + 0.5*z(7)*z(21)*(2*z(42)*z(164)-2*z(43)*z(166)-z(2)*z(192)-z(3)*z(194)) + 0.5*z(33)*(z(4)*z(193)+2*z(40)*z(168)+2*z(42)*z(164)-2*z(43)*z(166)-z(2)*z(192)-z(3)*z(194)) - z(15)*z(166)*(z(5)*z(42)+z(7)*z(40)) - rad*z(3)*z(149)*(z(5)*z(42)+z(7)*z(40)) - z(21)*z(168)*(z(5)*z(42)-z(6)*z(43)) - rad*z(4)*z(183)*(z(5)*z(42)-z(6)*z(43)) - z(195)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(6)*z(15)*(z(2)*z(192)-2*z(40)*z(168)-2*z(42)*z(164)-z(4)*z(193)) - 0.5*z(4)*z(21)*(2*z(42)*z(155)-2*z(43)*z(157)-z(5)*z(192)-z(6)*z(194)) - 0.5*z(30)*(z(7)*z(193)+2*z(40)*z(159)+2*z(42)*z(155)-2*z(43)*z(157)-z(5)*z(192)-z(6)*z(194));
z(197) = z(158) - z(6)*z(189) - z(22)*z(157);
z(198) = z(3)*z(189) + z(22)*z(166) - z(167);
z(199) = z(23)*z(159)*(z(2)*z(42)-z(3)*z(43)) + z(15)*z(155)*(z(3)*z(43)-z(4)*z(40)) + rad*z(5)*z(149)*(z(3)*z(43)-z(4)*z(40)) + z(197)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(2)*z(15)*(z(7)*z(193)+2*z(40)*z(159)-2*z(43)*z(157)-z(6)*z(194)) + 0.5*z(23)*z(7)*(2*z(42)*z(164)-2*z(43)*z(166)-z(2)*z(192)-z(3)*z(194)) + 0.5*z(34)*(z(4)*z(193)+2*z(40)*z(168)+2*z(42)*z(164)-2*z(43)*z(166)-z(2)*z(192)-z(3)*z(194)) - z(23)*z(168)*(z(5)*z(42)-z(6)*z(43)) - z(15)*z(164)*(z(6)*z(43)-z(7)*z(40)) - rad*z(2)*z(149)*(z(6)*z(43)-z(7)*z(40)) - z(198)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(5)*z(15)*(z(4)*z(193)+2*z(40)*z(168)-2*z(43)*z(166)-z(3)*z(194)) - 0.5*z(23)*z(4)*(2*z(42)*z(155)-2*z(43)*z(157)-z(5)*z(192)-z(6)*z(194)) - 0.5*z(31)*(z(7)*z(193)+2*z(40)*z(159)+2*z(42)*z(155)-2*z(43)*z(157)-z(5)*z(192)-z(6)*z(194));
z(200) = z(160) - z(7)*z(189) - z(22)*z(159);
z(201) = z(4)*z(189) + z(22)*z(168) - z(169);
z(202) = z(23)*z(166)*(z(5)*z(42)+z(7)*z(40)) + z(21)*z(155)*(z(3)*z(43)-z(4)*z(40)) + rad*z(5)*z(183)*(z(3)*z(43)-z(4)*z(40)) + z(200)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(23)*z(6)*(z(2)*z(192)-2*z(40)*z(168)-2*z(42)*z(164)-z(4)*z(193)) + 0.5*z(2)*z(21)*(z(7)*z(193)+2*z(40)*z(159)-2*z(43)*z(157)-z(6)*z(194)) + 0.5*z(35)*(z(4)*z(193)+2*z(40)*z(168)+2*z(42)*z(164)-2*z(43)*z(166)-z(2)*z(192)-z(3)*z(194)) - z(23)*z(157)*(z(2)*z(42)+z(4)*z(40)) - z(21)*z(164)*(z(6)*z(43)-z(7)*z(40)) - rad*z(2)*z(183)*(z(6)*z(43)-z(7)*z(40)) - z(201)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(3)*(z(5)*z(192)-2*z(40)*z(159)-2*z(42)*z(155)-z(7)*z(193)) - 0.5*z(5)*z(21)*(z(4)*z(193)+2*z(40)*z(168)-2*z(43)*z(166)-z(3)*z(194)) - 0.5*z(32)*(z(7)*z(193)+2*z(40)*z(159)+2*z(42)*z(155)-2*z(43)*z(157)-z(5)*z(192)-z(6)*z(194));
z(203) = z(196)*u1 + z(199)*u2 + z(202)*u3 + z(86) + z(51)*z(76) - z(52)*z(81);
z(204) = (z(45)*z(188)-z(49)*z(203))/z(50);
z(205) = (z(44)*z(188)-z(47)*z(203))/z(50);
z(206) = z(181) + z(136)*z(204) - z(137)*z(205);
a11 = -z(206)/z(139);
z(207) = e1*(z(12)+z(10)^2/z(11)^0.5)/z(12)^2;
z(208) = e4*(z(18)+z(9)^2/z(17)^0.5)/z(18)^2;
z(209) = e1*z(10)/(z(11)^0.5*z(12)^2);
z(210) = e4*z(9)/(z(17)^0.5*z(18)^2);
z(212) = 4*rad*z(21)*z(208) + 4*e4*z(21)*z(36) + 8*e1*z(15)*z(36) + 4*rad*z(9)*z(36)*z(208) + 8*rad*z(10)*z(36)*z(207) + 4*z(211)*z(8)*(z(210)+2*z(209)) + 4*rad*z(10)*z(15)*(z(210)+2*z(209)) - 8*rad*z(15)*z(207) - 4*z(23)*e3*z(36) - 4*rad*z(36)*(z(210)+2*z(209)) - 4*rad*z(9)*z(21)*(z(210)+2*z(209));
z(213) = z(23)*z(212)/(z(38)^0.5*z(39)^2);
z(214) = (8*rad*z(39)*z(207)+z(15)*z(212)/z(38)^0.5)/z(39)^2;
z(215) = (4*rad*z(39)*z(208)-z(21)*z(212)/z(38)^0.5)/z(39)^2;
z(216) = -2*e4*z(16) - 4*rad*z(7)*z(209);
z(217) = z(216) + 2*e4*z(22) - 2*rad*z(7)*z(210);
z(218) = 2*e3*z(16) - 4*rad*z(4)*z(209);
z(219) = 2*e3*z(22) + 2*rad*z(4)*z(210) - z(218);
z(220) = 2*z(23)*e2*(z(5)*z(42)+z(7)*z(40)) + 4*z(23)*e1*(z(2)*z(42)+z(4)*z(40)) + 2*e2*z(21)*(z(3)*z(43)-z(4)*z(40)) + 2*rad*z(5)*z(208)*(z(3)*z(43)-z(4)*z(40)) + z(217)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(23)*z(3)*(4*e2*z(42)-4*e4*z(40)-z(5)*z(213)-z(7)*z(214)) + 0.5*z(2)*z(21)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214)) + 0.5*z(35)*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214)) - 2*rad*z(2)*z(208)*(z(6)*z(43)-z(7)*z(40)) - z(219)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(6)*(4*e3*z(40)-z(2)*z(213)-z(4)*z(214)) - 0.5*z(5)*z(21)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214)) - 0.5*z(32)*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214));
z(221) = 2*e4*z(15) + 2*rad*z(10)*z(208) - 4*e1*z(21) - 4*rad*z(9)*z(207);
z(222) = 2*e2*z(16) - 4*rad*z(5)*z(209);
z(223) = z(222) - 2*e2*z(22) - 2*rad*z(5)*z(210);
z(224) = rad*z(2)*(z(210)+2*z(209));
z(225) = 4*rad*z(3)*z(207)*(z(5)*z(42)+z(7)*z(40)) + 2*rad*z(7)*z(208)*(z(2)*z(42)-z(3)*z(43)) + z(223)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(6)*z(15)*(4*e3*z(40)-z(2)*z(213)-z(4)*z(214)) + 0.5*z(33)*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214)) - 4*e1*z(15)*(z(2)*z(42)+z(4)*z(40)) - 2*e2*z(15)*(z(5)*z(42)+z(7)*z(40)) - 4*rad*z(6)*z(207)*(z(2)*z(42)+z(4)*z(40)) - 2*e3*z(21)*(z(5)*z(42)-z(6)*z(43)) - 2*e4*z(21)*(z(2)*z(42)-z(3)*z(43)) - 2*rad*z(4)*z(208)*(z(5)*z(42)-z(6)*z(43)) - 2*z(224)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(7)*z(21)*(z(2)*z(213)+z(3)*z(215)+4*e2*z(43)) - 0.5*z(3)*z(15)*(4*e2*z(42)-4*e4*z(40)-z(5)*z(213)-z(7)*z(214)) - 0.5*z(4)*z(21)*(4*e2*z(42)+8*e1*z(43)-z(5)*z(213)-z(6)*z(215)) - 0.5*z(30)*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214));
z(226) = -2*z(23)*e4 - 2*e3*z(21) - 2*rad*z(8)*z(208);
z(227) = -4*e1*z(16) - 4*rad*z(6)*z(209);
z(228) = z(227) + 4*e1*z(22) - 2*rad*z(6)*z(210);
z(229) = 2*e2*z(16) - 4*rad*z(3)*z(209);
z(230) = 2*e2*z(22) + 2*rad*z(3)*z(210) - z(229);
z(231) = 2*e2*z(15)*(z(3)*z(43)-z(4)*z(40)) + 4*rad*z(2)*z(207)*(z(6)*z(43)-z(7)*z(40)) + z(228)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(2)*z(15)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214)) + 0.5*z(34)*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214)) - 2*z(23)*e3*(z(5)*z(42)-z(6)*z(43)) - 2*z(23)*e4*(z(2)*z(42)-z(3)*z(43)) - 4*rad*z(5)*z(207)*(z(3)*z(43)-z(4)*z(40)) - z(230)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(7)*(z(2)*z(213)+z(3)*z(215)+4*e2*z(43)) - 0.5*z(23)*z(4)*(4*e2*z(42)+8*e1*z(43)-z(5)*z(213)-z(6)*z(215)) - 0.5*z(5)*z(15)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214)) - 0.5*z(31)*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214));
z(232) = 4*rad*z(8)*z(207) - 4*z(23)*e1 - 2*e3*z(15);
z(233) = z(44)*z(226) + z(49)*z(231) - z(45)*z(232) - z(47)*z(220);
z(234) = (z(233)*(z(45)*z(48)-z(46)*z(49))-z(50)*(z(45)*z(221)+z(48)*z(220)-z(46)*z(226)-z(49)*z(225)))/z(50)^2;
z(235) = (z(233)*(z(44)*z(48)-z(46)*z(47))+z(50)*(z(46)*z(232)+z(47)*z(225)-z(44)*z(221)-z(48)*z(231)))/z(50)^2;
z(236) = u1*(z(1)*z(234)*u3+2*z(1)*z(235)*u2-2*z(234)*z(102)-z(235)*z(101));
z(237) = u1*(z(1)*z(235)*u2+2*z(1)*z(234)*u3-2*z(235)*z(101)-z(234)*z(102));
z(238) = -4*rad*z(207)*u2 - z(15)*z(234)*u1;
z(239) = u2*z(238) - 4*rad*z(207)*u1^2 - z(234)*u1*z(109);
z(240) = (e2*z(235)+e3*z(234))*u1;
z(241) = -0.5*u3 - 0.5*e1*z(235)*u1 - 0.5*e4*z(234)*u1;
z(242) = -4*e2*z(241) - 2*e1*z(240) - 4*e1p;
z(243) = (z(10)*z(11)*z(242)-4*e1*z(11)*z(55)-4*e1*z(10)^2*z(55))/z(11)^1.5;
z(244) = (4*e1*z(10)*(z(12)*z(55)+2*z(10)*z(56))/z(11)^0.5+z(12)*(4*e1*z(56)-z(10)*z(243)-z(12)*z(242)))/z(12)^3;
z(245) = -u3*z(238) - z(235)*u1*z(109) - rad*u1*z(244);
z(246) = rad*(z(12)*z(243)-8*e1*z(10)*z(56)/z(11)^0.5)/z(12)^3;
z(247) = (-1+e2*z(234)-e3*z(235))*u1;
z(248) = 0.5*u2 + 0.5*e4*z(235)*u1 - 0.5*e1*z(234)*u1;
z(249) = e2*z(240) + e3*z(247) + 2*e1*z(241) + 2*e4*z(248) + 2*e2p;
z(250) = z(5)*z(246) + z(16)*z(249) + 2*e2*z(61) - 4*rad*z(209)*z(68);
z(251) = -4*e3*z(248) - 2*e1*z(240) - 4*e1p;
z(252) = z(6)*z(246) + z(16)*z(251) - 4*e1*z(61) - 4*rad*z(209)*z(71);
z(253) = 2*e2*z(248) + 2*e3*z(241) - e1*z(247) - e4*z(240) - 2*e4p;
z(254) = z(7)*z(246) + z(16)*z(253) - 2*e4*z(61) - 4*rad*z(209)*z(53);
z(255) = u1*z(250) + u2*z(252) + u3*z(254) + z(235)*u1*z(77) - z(234)*u1*z(72);
z(256) = rad*(u2*z(244)+z(234)*u1*z(57));
z(257) = 4*rad*z(207)*u1*u3 + z(256) - z(235)*u1*z(110);
z(258) = -4*e2*z(241) - 4*e3*z(248);
z(259) = z(16)*z(258) + z(2)*z(246) - 4*rad*z(209)*z(54);
z(260) = e2*z(240) + 2*e1*z(241) + 2*e2p - 2*e4*z(248) - e3*z(247);
z(261) = z(3)*z(246) + z(16)*z(260) + 2*e2*z(61) - 4*rad*z(209)*z(66);
z(262) = e2*z(247) + e3*z(240) + 2*e1*z(248) + 2*e4*z(241) + 2*e3p;
z(263) = z(4)*z(246) + z(16)*z(262) + 2*e3*z(61) - 4*rad*z(209)*z(69);
z(264) = z(234)*u1*z(74) - u1*z(259) - u2*z(261) - u3*z(263) - z(235)*u1*z(79);
z(265) = (z(23)*z(235)+4*rad*z(207))*u1;
z(266) = u2*z(238) - u1*z(265) - z(234)*u1*z(109);
z(267) = z(23)*z(234)*u1^2 - u3*z(238) - z(235)*u1*z(109) - rad*u1*z(244);
z(268) = u3*z(265) + z(256) - z(23)*z(234)*u1*u2 - z(234)*u1*z(121) - z(235)*u1*z(120);
z(269) = 0.25*z(236) + 0.25*z(237) - m*(z(7)*z(27)*z(239)+z(5)*z(27)*z(257)+z(6)*z(27)*z(245)+2*e2*z(27)*z(117)+4*rad*z(207)*z(119)+8*e1*z(15)*z(116)+8*rad*z(6)*z(207)*z(116)-2*e3*z(24)*z(118)-2*e4*z(27)*z(118)-z(4)*z(24)*z(239)-z(4)*z(222)*z(118)-4*rad*z(2)*z(7)*z(209)*z(118)-2*z(24)*z(255)-2*z(222)*z(115)-z(15)*z(245)-4*e1*z(27)*z(119)-2*e2*z(24)*z(119)-2*z(6)*z(15)*z(264)-z(2)*z(24)*z(257)-z(2)*z(222)*z(117)-z(3)*z(24)*z(245)-z(3)*z(222)*z(119)-4*rad*z(2)*z(5)*z(209)*z(117)-4*rad*z(2)*z(6)*z(209)*z(119)) - m*(z(7)*z(27)*z(266)+2*z(27)*z(264)+z(5)*z(27)*z(268)+z(6)*z(27)*z(267)+2*e2*z(27)*z(122)+4*rad*z(207)*z(124)+8*rad*z(3)*z(207)*z(115)-2*e3*z(24)*z(123)-2*e4*z(27)*z(123)-z(4)*z(24)*z(266)-z(4)*z(222)*z(123)-4*rad*z(2)*z(7)*z(209)*z(123)-z(15)*z(267)-4*e1*z(27)*z(124)-4*e2*z(15)*z(115)-2*e2*z(24)*z(124)-2*z(3)*z(15)*z(255)-z(2)*z(24)*z(268)-z(2)*z(222)*z(122)-z(3)*z(24)*z(267)-z(3)*z(222)*z(124)-8*rad*z(2)*z(209)*z(116)-4*rad*z(2)*z(5)*z(209)*z(122)-4*rad*z(2)*z(6)*z(209)*z(124));
z(270) = z(235)*u1*(2*z(1)*u1-z(100));
z(271) = z(235)*u1*(z(1)*u1-z(100));
z(272) = -0.25*z(270) - 0.25*z(271) - m*(z(7)*z(28)*z(239)+z(7)*z(229)*z(118)+z(15)*z(257)+2*z(28)*z(264)+2*z(229)*z(116)+z(5)*z(28)*z(257)+z(5)*z(229)*z(117)+z(6)*z(28)*z(245)+z(6)*z(229)*z(119)+2*e2*z(28)*z(117)+2*z(2)*z(15)*z(255)-2*e3*z(25)*z(118)-2*e4*z(28)*z(118)-z(4)*z(25)*z(239)-z(4)*z(227)*z(118)-2*z(25)*z(255)-2*z(227)*z(115)-4*rad*z(207)*z(117)-4*e1*z(28)*z(119)-2*e2*z(25)*z(119)-z(2)*z(25)*z(257)-z(2)*z(227)*z(117)-z(3)*z(25)*z(245)-z(3)*z(227)*z(119)-8*rad*z(2)*z(207)*z(115)) - m*(z(7)*z(28)*z(266)+z(7)*z(229)*z(123)+z(15)*z(268)+z(5)*z(28)*z(268)+z(5)*z(229)*z(122)+z(6)*z(28)*z(267)+z(6)*z(229)*z(124)+2*z(23)*e4*z(116)+2*e2*z(28)*z(122)+2*z(5)*z(15)*z(264)+4*e2*z(15)*z(116)-z(23)*z(266)-2*e3*z(25)*z(123)-2*e4*z(28)*z(123)-z(4)*z(25)*z(266)-z(4)*z(227)*z(123)-4*rad*z(207)*z(122)-4*e1*z(28)*z(124)-2*z(23)*e3*z(115)-2*e2*z(25)*z(124)-z(23)*z(4)*z(255)-z(23)*z(7)*z(264)-z(2)*z(25)*z(268)-z(2)*z(227)*z(122)-z(3)*z(25)*z(267)-z(3)*z(227)*z(124)-8*rad*z(5)*z(207)*z(116));
z(273) = z(234)*u1*(z(1)*u1-z(100));
z(274) = z(234)*u1*(2*z(1)*u1-z(100));
z(275) = m*(z(4)*z(26)*z(239)+z(4)*z(216)*z(118)+2*e3*z(26)*z(118)+2*e4*z(29)*z(118)+2*z(26)*z(255)+2*z(216)*z(115)+z(2)*z(26)*z(257)+z(2)*z(216)*z(117)+z(3)*z(26)*z(245)+z(3)*z(216)*z(119)+2*e2*z(26)*z(119)+4*e1*z(29)*z(119)-z(7)*z(29)*z(239)-z(7)*z(218)*z(118)-2*z(29)*z(264)-2*z(218)*z(116)-2*e2*z(29)*z(117)-z(5)*z(29)*z(257)-z(5)*z(218)*z(117)-z(6)*z(29)*z(245)-z(6)*z(218)*z(119)) + m*(z(4)*z(26)*z(266)+z(4)*z(216)*z(123)+2*e3*z(26)*z(123)+2*e4*z(29)*z(123)+z(2)*z(26)*z(268)+z(2)*z(216)*z(122)+z(3)*z(26)*z(267)+z(3)*z(216)*z(124)+2*e2*z(26)*z(124)+4*z(23)*e1*z(116)+4*e1*z(29)*z(124)-z(7)*z(29)*z(266)-z(7)*z(218)*z(123)-z(23)*z(267)-2*z(23)*e2*z(115)-2*e2*z(29)*z(122)-z(23)*z(3)*z(255)-z(23)*z(6)*z(264)-z(5)*z(29)*z(268)-z(5)*z(218)*z(122)-z(6)*z(29)*z(267)-z(6)*z(218)*z(124)) - 0.25*z(273) - 0.25*z(274);
z(276) = z(269) + z(51)*z(272) + z(235)*z(134) - z(52)*z(275) - z(234)*z(131);
z(277) = m*(2*z(24)*z(227)+2*z(25)*z(222)+2*z(27)*z(229)+z(23)*z(4)*z(222)+2*z(23)*e3*z(24)+2*z(23)*e4*z(27)+2*z(3)*z(15)*z(227)+4*e2*z(15)*z(25)+4*e2*z(15)*z(27)+8*e1*z(15)*z(28)+4*z(211)*z(2)*z(7)*z(209)+8*rad*z(2)*z(24)*z(207)+8*rad*z(6)*z(28)*z(207)-2*z(2)*z(15)*z(222)-2*z(6)*z(15)*z(229)-8*rad*z(2)*z(28)*z(209)-8*rad*z(3)*z(25)*z(207)-8*rad*z(5)*z(27)*z(207)-8*rad*z(2)*z(5)*z(15)*z(209));
z(278) = m*(2*z(25)*z(227)+2*z(28)*z(229)+z(23)*z(4)*z(227)+2*z(23)*e3*z(25)+2*z(23)*e4*z(28)+2*z(5)*z(15)*z(229)+4*e2*z(15)*z(28)+8*rad*z(2)*z(25)*z(207)-8*rad*z(15)*z(207)-2*z(2)*z(15)*z(227)-z(23)*z(7)*z(229)-8*rad*z(5)*z(28)*z(207));
z(279) = m*(2*z(25)*z(216)+2*z(26)*z(227)+2*z(28)*z(218)+2*z(29)*z(229)+z(23)*z(4)*z(216)+z(23)*z(6)*z(229)+2*z(23)*e3*z(26)+2*z(23)*e4*z(29)+2*z(5)*z(15)*z(218)+4*e2*z(15)*z(29)+8*rad*z(2)*z(26)*z(207)-4*z(23)*e1*z(28)-2*z(23)*e2*z(25)-2*z(2)*z(15)*z(216)-z(23)*z(3)*z(227)-z(23)*z(7)*z(218)-8*rad*z(5)*z(29)*z(207));
z(280) = z(277) + 2*z(51)*z(278) - z(52)*z(279) - z(129)*z(234) - z(130)*z(235);
z(281) = e1*z(247) + e4*z(240) + 2*e2*z(248) + 2*e3*z(241) + 2*e4p;
z(282) = (z(9)*z(17)*z(281)+2*e4*z(17)*z(58)+2*e4*z(9)^2*z(58))/z(17)^1.5;
z(283) = (z(18)*(z(9)*z(282)+z(18)*z(281)+2*e4*z(59))+2*e4*z(9)*(z(18)*z(58)+2*z(9)*z(59))/z(17)^0.5)/z(18)^3;
z(284) = z(15)*z(281) + z(21)*z(242) + rad*z(10)*z(283) + 2*rad*e4*z(57) + 2*rad*z(208)*z(55) - 4*rad*e1*z(60) - 4*rad*z(207)*z(58) - rad*z(9)*z(244);
z(285) = e3*z(240) + 2*e1*z(248) + 2*e3p - 2*e4*z(241) - e2*z(247);
z(286) = z(23)*z(242) + rad*z(8)*z(244) + 4*rad*z(207)*z(63) - z(15)*z(285) - 2*rad*e3*z(57);
z(287) = -z(23)*z(281) - z(21)*z(285) - 2*rad*e3*z(60) - 2*rad*z(208)*z(63) - rad*z(8)*z(283);
z(288) = u1*z(284) + u2*z(286) + u3*z(287) + z(235)*u1*z(90) - z(234)*u1*z(88);
z(289) = rad*(z(18)*z(282)+4*e4*z(9)*z(59)/z(17)^0.5)/z(18)^3;
z(290) = z(250) - 2*e2*z(62) - z(5)*z(289) - z(22)*z(249) - 2*rad*z(210)*z(68);
z(292) = 2*rad*z(21)*z(283) + 2*z(21)*z(36)*z(281) + 4*z(291)*z(208)*z(60) + 2*rad*z(9)*z(36)*z(283) + 2*rad*z(10)*z(36)*z(244) + 4*rad*e4*z(36)*z(60) + 4*rad*z(36)*z(208)*z(58) + 8*rad*e1*z(36)*z(57) + 8*rad*z(36)*z(207)*z(55) + 2*z(36)*(z(246)-z(289)) + 2*z(9)*z(21)*(z(246)-z(289)) + 4*z(211)*(z(210)+2*z(209))*z(63) + 4*e4*z(21)*(z(61)-z(62)) + 8*e1*z(15)*(z(61)-z(62)) + 4*rad*z(9)*z(208)*(z(61)-z(62)) + 4*rad*z(15)*(z(210)+2*z(209))*z(55) + 4*z(291)*z(10)*(z(210)+2*z(209))*z(57) + 8*rad*z(10)*z(207)*(z(61)-z(62)) - 8*z(291)*z(207)*z(57) - 2*rad*z(15)*z(244) - 2*z(23)*z(36)*z(285) - 2*z(15)*z(36)*z(242) - 4*z(23)*e3*(z(61)-z(62)) - 2*z(23)*z(8)*(z(246)-z(289)) - 2*z(10)*z(15)*(z(246)-z(289)) - 4*rad*z(21)*(z(210)+2*z(209))*z(58) - 4*z(291)*z(9)*(z(210)+2*z(209))*z(60) - 4*rad*(z(210)+2*z(209))*(z(61)-z(62));
z(293) = z(23)*(2*z(292)/(z(38)^0.5*z(39)^2)-z(212)*(2+z(39)/z(38)^0.5)*z(64)/(z(38)*z(39)^3));
z(294) = (2*z(212)*(rad*z(39)*z(57)-z(15)*z(64)/z(38)^0.5)/z(38)^0.5+z(39)*(4*rad*z(39)*z(244)+(2*z(15)*z(38)*z(292)-z(15)*z(212)*z(64)-8*rad*z(38)*z(207)*z(64))/z(38)^1.5))/z(39)^3;
z(295) = (2*z(21)*z(212)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(283)+(z(21)*z(212)*z(64)-2*z(21)*z(38)*z(292)-4*rad*z(38)*z(208)*z(64)-2*rad*z(38)*z(212)*z(60))/z(38)^1.5))/z(39)^3;
z(296) = z(22)*z(258) + z(2)*z(289) + 2*rad*z(210)*z(54) - z(259);
z(297) = z(15)*(z(2)*z(42)+z(4)*z(40))*z(251) + rad*z(3)*(z(5)*z(42)+z(7)*z(40))*z(244) + z(21)*(z(2)*z(42)-z(3)*z(43))*z(253) + 4*rad*z(207)*(z(5)*z(42)+z(7)*z(40))*z(66) + rad*z(7)*(z(2)*z(42)-z(3)*z(43))*z(283) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(290) + 2*rad*z(208)*(z(2)*z(42)-z(3)*z(43))*z(53) + 0.5*z(15)*(4*e3*z(40)-z(2)*z(213)-z(4)*z(214))*z(71) + 0.5*rad*z(6)*(4*e3*z(40)-z(2)*z(213)-z(4)*z(214))*z(57) + e2*z(15)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) + 2*e1*z(15)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) + 2*rad*z(6)*z(207)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) + rad*z(7)*z(208)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214))*z(85) + z(224)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.25*z(6)*z(15)*(4*z(42)*z(258)+4*e3*z(70)+4*z(40)*z(262)-2*z(213)*z(54)-2*z(214)*z(69)-z(2)*z(293)-z(4)*z(294)) + 0.25*z(7)*z(21)*(4*z(42)*z(258)-4*e2*z(67)-4*z(43)*z(260)-2*z(213)*z(54)-2*z(215)*z(66)-z(2)*z(293)-z(3)*z(295)) + 0.25*z(33)*(4*z(42)*z(258)+4*e3*z(70)+4*z(40)*z(262)-4*e2*z(67)-4*z(43)*z(260)-2*z(213)*z(54)-2*z(214)*z(69)-2*z(215)*z(66)-z(2)*z(293)-z(3)*z(295)-z(4)*z(294)) - z(15)*(z(5)*z(42)+z(7)*z(40))*z(260) - 4*rad*e1*(z(2)*z(42)+z(4)*z(40))*z(57) - 4*rad*z(207)*(z(2)*z(42)+z(4)*z(40))*z(71) - 2*rad*e2*(z(5)*z(42)+z(7)*z(40))*z(57) - rad*z(6)*(z(2)*z(42)+z(4)*z(40))*z(244) - z(21)*(z(5)*z(42)-z(6)*z(43))*z(262) - 2*rad*e3*(z(5)*z(42)-z(6)*z(43))*z(60) - 2*rad*e4*(z(2)*z(42)-z(3)*z(43))*z(60) - 2*rad*z(208)*(z(5)*z(42)-z(6)*z(43))*z(69) - rad*z(4)*(z(5)*z(42)-z(6)*z(43))*z(283) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(296) - 0.5*z(21)*(z(2)*z(213)+z(3)*z(215)+4*e2*z(43))*z(53) - 0.5*rad*z(7)*(z(2)*z(213)+z(3)*z(215)+4*e2*z(43))*z(60) - 2*rad*z(3)*z(207)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e3*z(21)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - e4*z(21)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - 0.5*z(15)*(4*e2*z(42)-4*e4*z(40)-z(5)*z(213)-z(7)*z(214))*z(66) - 0.5*z(21)*(4*e2*z(42)+8*e1*z(43)-z(5)*z(213)-z(6)*z(215))*z(69) - rad*z(4)*z(208)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*rad*z(3)*(4*e2*z(42)-4*e4*z(40)-z(5)*z(213)-z(7)*z(214))*z(57) - 0.5*rad*z(4)*(4*e2*z(42)+8*e1*z(43)-z(5)*z(213)-z(6)*z(215))*z(60) - 0.5*z(223)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214))*z(83) - 0.25*z(3)*z(15)*(4*z(40)*z(253)+4*z(42)*z(249)-4*e2*z(65)-4*e4*z(70)-2*z(213)*z(68)-2*z(214)*z(53)-z(5)*z(293)-z(7)*z(294)) - 0.25*z(4)*z(21)*(4*z(42)*z(249)+8*e1*z(67)-4*e2*z(65)-4*z(43)*z(251)-2*z(213)*z(68)-2*z(215)*z(71)-z(5)*z(293)-z(6)*z(295)) - 0.25*z(30)*(4*z(40)*z(253)+4*z(42)*z(249)+8*e1*z(67)-4*e2*z(65)-4*e4*z(70)-4*z(43)*z(251)-2*z(213)*z(68)-2*z(214)*z(53)-2*z(215)*z(71)-z(5)*z(293)-z(6)*z(295)-z(7)*z(294));
z(298) = z(252) + 4*e1*z(62) - z(6)*z(289) - z(22)*z(251) - 2*rad*z(210)*z(71);
z(299) = z(3)*z(289) + z(22)*z(260) + 2*e2*z(62) + 2*rad*z(210)*z(66) - z(261);
z(300) = z(23)*(z(2)*z(42)-z(3)*z(43))*z(253) + z(15)*(z(3)*z(43)-z(4)*z(40))*z(249) + rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(244) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(298) + 2*rad*e2*(z(3)*z(43)-z(4)*z(40))*z(57) + 4*rad*z(207)*(z(6)*z(43)-z(7)*z(40))*z(54) + 0.5*z(15)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214))*z(54) + 2*rad*z(5)*z(207)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) + 0.5*rad*z(2)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214))*z(57) + 0.5*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214))*z(75) + 0.5*z(230)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.25*z(23)*z(7)*(4*z(42)*z(258)-4*e2*z(67)-4*z(43)*z(260)-2*z(213)*z(54)-2*z(215)*z(66)-z(2)*z(293)-z(3)*z(295)) + 0.25*z(2)*z(15)*(4*z(40)*z(253)+8*e1*z(67)-4*e4*z(70)-4*z(43)*z(251)-2*z(214)*z(53)-2*z(215)*z(71)-z(6)*z(295)-z(7)*z(294)) + 0.25*z(34)*(4*z(42)*z(258)+4*e3*z(70)+4*z(40)*z(262)-4*e2*z(67)-4*z(43)*z(260)-2*z(213)*z(54)-2*z(214)*z(69)-2*z(215)*z(66)-z(2)*z(293)-z(3)*z(295)-z(4)*z(294)) - z(15)*(z(6)*z(43)-z(7)*z(40))*z(258) - z(23)*(z(5)*z(42)-z(6)*z(43))*z(262) - 4*rad*z(207)*(z(3)*z(43)-z(4)*z(40))*z(68) - rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(244) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(299) - 0.5*z(23)*(z(2)*z(213)+z(3)*z(215)+4*e2*z(43))*z(53) - e2*z(15)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 2*rad*z(2)*z(207)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - z(23)*e3*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - z(23)*e4*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - 0.5*z(23)*(4*e2*z(42)+8*e1*z(43)-z(5)*z(213)-z(6)*z(215))*z(69) - 0.5*z(15)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214))*z(68) - 0.5*rad*z(5)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214))*z(57) - 0.5*z(228)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214))*z(73) - 0.25*z(23)*z(4)*(4*z(42)*z(249)+8*e1*z(67)-4*e2*z(65)-4*z(43)*z(251)-2*z(213)*z(68)-2*z(215)*z(71)-z(5)*z(293)-z(6)*z(295)) - 0.25*z(5)*z(15)*(4*e3*z(70)+4*z(40)*z(262)-4*e2*z(67)-4*z(43)*z(260)-2*z(214)*z(69)-2*z(215)*z(66)-z(3)*z(295)-z(4)*z(294)) - 0.25*z(31)*(4*z(40)*z(253)+4*z(42)*z(249)+8*e1*z(67)-4*e2*z(65)-4*e4*z(70)-4*z(43)*z(251)-2*z(213)*z(68)-2*z(214)*z(53)-2*z(215)*z(71)-z(5)*z(293)-z(6)*z(295)-z(7)*z(294));
z(301) = z(254) + 2*e4*z(62) - z(7)*z(289) - z(22)*z(253) - 2*rad*z(210)*z(53);
z(302) = z(4)*z(289) + z(22)*z(262) + 2*e3*z(62) + 2*rad*z(210)*z(69) - z(263);
z(303) = z(23)*(z(5)*z(42)+z(7)*z(40))*z(260) + z(21)*(z(3)*z(43)-z(4)*z(40))*z(249) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(283) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(301) + 2*rad*e2*(z(3)*z(43)-z(4)*z(40))*z(60) + 2*rad*z(208)*(z(3)*z(43)-z(4)*z(40))*z(68) + rad*z(2)*z(208)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(23)*(4*e2*z(42)-4*e4*z(40)-z(5)*z(213)-z(7)*z(214))*z(66) + 0.5*z(21)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214))*z(54) + 0.5*rad*z(2)*(8*e1*z(43)-4*e4*z(40)-z(6)*z(215)-z(7)*z(214))*z(60) + 0.5*(4*e3*z(40)-4*e2*z(43)-z(2)*z(213)-z(3)*z(215)-z(4)*z(214))*z(80) + 0.5*z(219)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.25*z(23)*z(3)*(4*z(40)*z(253)+4*z(42)*z(249)-4*e2*z(65)-4*e4*z(70)-2*z(213)*z(68)-2*z(214)*z(53)-z(5)*z(293)-z(7)*z(294)) + 0.25*z(2)*z(21)*(4*z(40)*z(253)+8*e1*z(67)-4*e4*z(70)-4*z(43)*z(251)-2*z(214)*z(53)-2*z(215)*z(71)-z(6)*z(295)-z(7)*z(294)) + 0.25*z(35)*(4*z(42)*z(258)+4*e3*z(70)+4*z(40)*z(262)-4*e2*z(67)-4*z(43)*z(260)-2*z(213)*z(54)-2*z(214)*z(69)-2*z(215)*z(66)-z(2)*z(293)-z(3)*z(295)-z(4)*z(294)) - z(21)*(z(6)*z(43)-z(7)*z(40))*z(258) - z(23)*(z(2)*z(42)+z(4)*z(40))*z(251) - 2*rad*z(208)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(283) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(302) - 0.5*z(23)*(4*e3*z(40)-z(2)*z(213)-z(4)*z(214))*z(71) - 2*z(23)*e1*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) - z(23)*e2*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e2*z(21)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - rad*z(5)*z(208)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(21)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214))*z(68) - 0.5*rad*z(5)*(4*e3*z(40)-4*e2*z(43)-z(3)*z(215)-z(4)*z(214))*z(60) - 0.5*z(217)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(42)+8*e1*z(43)-4*e4*z(40)-z(5)*z(213)-z(6)*z(215)-z(7)*z(214))*z(78) - 0.25*z(23)*z(6)*(4*z(42)*z(258)+4*e3*z(70)+4*z(40)*z(262)-2*z(213)*z(54)-2*z(214)*z(69)-z(2)*z(293)-z(4)*z(294)) - 0.25*z(5)*z(21)*(4*e3*z(70)+4*z(40)*z(262)-4*e2*z(67)-4*z(43)*z(260)-2*z(214)*z(69)-2*z(215)*z(66)-z(3)*z(295)-z(4)*z(294)) - 0.25*z(32)*(4*z(40)*z(253)+4*z(42)*z(249)+8*e1*z(67)-4*e2*z(65)-4*e4*z(70)-4*z(43)*z(251)-2*z(213)*z(68)-2*z(214)*z(53)-2*z(215)*z(71)-z(5)*z(293)-z(6)*z(295)-z(7)*z(294));
z(304) = u1*z(297) + u2*z(300) + u3*z(303) + z(235)*u1*z(81) - z(234)*u1*z(76);
z(305) = (z(233)*(z(45)*z(91)-z(49)*z(87))-z(50)*(z(45)*z(288)+z(220)*z(91)-z(49)*z(304)-z(226)*z(87)))/z(50)^2;
z(306) = m*(2*z(24)*z(216)+2*z(26)*z(222)+2*z(27)*z(218)+4*z(211)*z(207)+2*z(3)*z(15)*z(216)+4*e2*z(15)*z(26)+8*e1*z(15)*z(29)+8*rad*z(6)*z(29)*z(207)-4*z(23)*e1*z(27)-2*z(23)*e2*z(24)-2*z(6)*z(15)*z(218)-z(23)*z(3)*z(222)-8*rad*z(2)*z(29)*z(209)-8*rad*z(3)*z(26)*z(207)-4*z(211)*z(2)*z(6)*z(209));
z(307) = 4*m*(z(26)*z(216)+z(29)*z(218)) - 2*z(132)*(z(3)*z(216)+2*e2*z(26)+4*e1*z(29)-z(6)*z(218));
z(308) = z(306) + z(51)*z(279) + z(130)*z(234) + z(133)*z(235) - z(52)*z(307);
z(309) = (z(233)*(z(44)*z(91)-z(47)*z(87))-z(50)*(z(44)*z(288)+z(231)*z(91)-z(47)*z(304)-z(232)*z(87)))/z(50)^2;
z(310) = z(276) + z(137)*z(309) + z(280)*z(92) - z(136)*z(305) - z(308)*z(93);
z(311) = -z(15)*z(234) - 4*rad*z(51)*z(207);
z(312) = -4*rad*z(207) - z(23)*z(235);
z(313) = m*(z(9)*z(312)+2*e4*z(15)+2*e4*z(97)-4*e1*z(98)-4*e3*z(96)-2*z(8)*z(311)-4*rad*z(9)*z(207)-z(23)*z(10)*z(234));
z(314) = m*(4*rad*z(15)*z(207)+4*rad*z(2)*z(27)*z(209)+4*rad*z(3)*z(24)*z(207)-z(24)*z(222)-4*e1*z(15)*z(27)-2*e2*z(15)*z(24)-z(3)*z(15)*z(222)-4*rad*z(6)*z(27)*z(207)-4*rad*z(2)*z(6)*z(15)*z(209));
z(315) = z(51)*z(277) + z(127)*z(235) - 4*z(314) - z(52)*z(306) - z(126)*z(234);
z(316) = z(315) + z(51)*z(280) + z(137)*z(235) - z(52)*z(308) - z(136)*z(234);
z(317) = (z(316)*z(141)-z(139)*(9.810000000000001*z(313)+z(310)))/z(139)^2;
a12 = z(317);
z(318) = e2*(z(12)+z(10)^2/z(11)^0.5)/z(12)^2;
z(319) = e3*(z(18)+z(9)^2/z(17)^0.5)/z(18)^2;
z(320) = e2*z(10)/(z(11)^0.5*z(12)^2);
z(321) = e3*z(9)/(z(17)^0.5*z(18)^2);
z(322) = 4*rad*z(21)*z(319) + 4*z(23)*e4*z(36) + 4*e3*z(21)*z(36) + 8*e2*z(15)*z(36) + 4*rad*z(9)*z(36)*z(319) + 8*rad*z(10)*z(36)*z(318) + 4*z(211)*z(8)*(z(321)+2*z(320)) + 4*rad*z(10)*z(15)*(z(321)+2*z(320)) - 8*rad*z(15)*z(318) - 4*rad*z(36)*(z(321)+2*z(320)) - 4*rad*z(9)*z(21)*(z(321)+2*z(320));
z(323) = z(23)*z(322)/(z(38)^0.5*z(39)^2);
z(324) = (8*rad*z(39)*z(318)+z(15)*z(322)/z(38)^0.5)/z(39)^2;
z(325) = (4*rad*z(39)*z(319)-z(21)*z(322)/z(38)^0.5)/z(39)^2;
z(326) = 2*e3*z(16) - 4*rad*z(7)*z(320);
z(327) = z(326) - 2*e3*z(22) - 2*rad*z(7)*z(321);
z(328) = 2*e4*z(16) - 4*rad*z(4)*z(320);
z(329) = 2*e4*z(22) + 2*rad*z(4)*z(321) - z(328);
z(330) = 2*z(23)*e1*(z(5)*z(42)+z(7)*z(40)) + 2*e1*z(21)*(z(3)*z(43)-z(4)*z(40)) + 4*e2*z(21)*(z(6)*z(43)-z(7)*z(40)) + 2*rad*z(5)*z(319)*(z(3)*z(43)-z(4)*z(40)) + z(327)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(2)*z(21)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324)) + 0.5*z(23)*z(3)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(7)*z(324)) + 0.5*z(35)*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324)) - 2*rad*z(2)*z(319)*(z(6)*z(43)-z(7)*z(40)) - z(329)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(6)*(4*e4*z(40)-8*e2*z(42)-z(2)*z(323)-z(4)*z(324)) - 0.5*z(5)*z(21)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324)) - 0.5*z(32)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324));
z(331) = 2*e3*z(15) + 2*rad*z(10)*z(319) - 4*e2*z(21) - 4*rad*z(9)*z(318);
z(332) = 2*e1*z(16) - 4*rad*z(5)*z(320);
z(333) = z(332) - 2*e1*z(22) - 2*rad*z(5)*z(321);
z(334) = -4*e2*z(16) - 4*rad*z(2)*z(320);
z(335) = 2*rad*z(2)*z(321) - z(334) - 4*e2*z(22);
z(336) = 2*e3*z(21)*(z(2)*z(42)-z(3)*z(43)) + 4*rad*z(3)*z(318)*(z(5)*z(42)+z(7)*z(40)) + 2*rad*z(7)*z(319)*(z(2)*z(42)-z(3)*z(43)) + z(333)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(6)*z(15)*(4*e4*z(40)-8*e2*z(42)-z(2)*z(323)-z(4)*z(324)) + 0.5*z(33)*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324)) - 2*e1*z(15)*(z(5)*z(42)+z(7)*z(40)) - 4*rad*z(6)*z(318)*(z(2)*z(42)+z(4)*z(40)) - 2*e4*z(21)*(z(5)*z(42)-z(6)*z(43)) - 2*rad*z(4)*z(319)*(z(5)*z(42)-z(6)*z(43)) - z(335)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(4)*z(21)*(4*e1*z(42)-z(5)*z(323)-z(6)*z(325)) - 0.5*z(7)*z(21)*(z(2)*z(323)+z(3)*z(325)+4*e1*z(43)+8*e2*z(42)) - 0.5*z(3)*z(15)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(7)*z(324)) - 0.5*z(30)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324));
z(337) = 2*e4*z(21) - 2*z(23)*e3 - 2*rad*z(8)*z(319);
z(338) = rad*z(6)*(z(321)+2*z(320));
z(339) = 2*e1*z(16) - 4*rad*z(3)*z(320);
z(340) = 2*e1*z(22) + 2*rad*z(3)*z(321) - z(339);
z(341) = 2*z(23)*e3*(z(2)*z(42)-z(3)*z(43)) + 2*e1*z(15)*(z(3)*z(43)-z(4)*z(40)) + 4*e2*z(15)*(z(6)*z(43)-z(7)*z(40)) + 4*rad*z(2)*z(318)*(z(6)*z(43)-z(7)*z(40)) + 0.5*z(2)*z(15)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324)) + 0.5*z(34)*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324)) - 2*z(23)*e4*(z(5)*z(42)-z(6)*z(43)) - 4*rad*z(5)*z(318)*(z(3)*z(43)-z(4)*z(40)) - 2*z(338)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) - z(340)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(4)*(4*e1*z(42)-z(5)*z(323)-z(6)*z(325)) - 0.5*z(23)*z(7)*(z(2)*z(323)+z(3)*z(325)+4*e1*z(43)+8*e2*z(42)) - 0.5*z(5)*z(15)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324)) - 0.5*z(31)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324));
z(342) = 2*e4*z(15) + 4*rad*z(8)*z(318) - 4*z(23)*e2;
z(343) = z(44)*z(337) + z(49)*z(341) - z(45)*z(342) - z(47)*z(330);
z(344) = (z(343)*(z(45)*z(48)-z(46)*z(49))-z(50)*(z(45)*z(331)+z(48)*z(330)-z(46)*z(337)-z(49)*z(336)))/z(50)^2;
z(345) = (z(343)*(z(44)*z(48)-z(46)*z(47))+z(50)*(z(46)*z(342)+z(47)*z(336)-z(44)*z(331)-z(48)*z(341)))/z(50)^2;
z(346) = u1*(z(1)*z(344)*u3+2*z(1)*z(345)*u2-2*z(344)*z(102)-z(345)*z(101));
z(347) = u1*(z(1)*z(345)*u2+2*z(1)*z(344)*u3-2*z(345)*z(101)-z(344)*z(102));
z(348) = -4*rad*z(318)*u2 - z(15)*z(344)*u1;
z(349) = u2*z(348) - 4*rad*z(318)*u1^2 - z(344)*u1*z(109);
z(350) = 0.5*u3 + 0.5*e2*z(345)*u1 + 0.5*e3*z(344)*u1;
z(351) = (e1*z(345)+e4*z(344))*u1;
z(352) = 2*e2*z(351) - 4*e1*z(350) - 4*e2p;
z(353) = (z(10)*z(11)*z(352)-4*e2*z(11)*z(55)-4*e2*z(10)^2*z(55))/z(11)^1.5;
z(354) = (4*e2*z(10)*(z(12)*z(55)+2*z(10)*z(56))/z(11)^0.5+z(12)*(4*e2*z(56)-z(10)*z(353)-z(12)*z(352)))/z(12)^3;
z(355) = -u3*z(348) - z(345)*u1*z(109) - rad*u1*z(354);
z(356) = rad*(z(12)*z(353)-8*e2*z(10)*z(56)/z(11)^0.5)/z(12)^3;
z(357) = 0.5*e2*z(344)*u1 - 0.5*u2 - 0.5*e3*z(345)*u1;
z(358) = (1+e1*z(344)-e4*z(345))*u1;
z(359) = 2*e2*z(350) + 2*e3*z(357) + 2*e1p - e1*z(351) - e4*z(358);
z(360) = z(5)*z(356) + z(16)*z(359) + 2*e1*z(61) - 4*rad*z(320)*z(68);
z(361) = 2*e3*z(358) - 4*e1*z(350);
z(362) = z(16)*z(361) + z(6)*z(356) - 4*rad*z(320)*z(71);
z(363) = 2*e3p - 2*e1*z(357) - 2*e4*z(350) - e2*z(358) - e3*z(351);
z(364) = z(7)*z(356) + z(16)*z(363) + 2*e3*z(61) - 4*rad*z(320)*z(53);
z(365) = u1*z(360) + u2*z(362) + u3*z(364) + z(345)*u1*z(77) - z(344)*u1*z(72);
z(366) = rad*(u2*z(354)+z(344)*u1*z(57));
z(367) = 4*rad*z(318)*u1*u3 + z(366) - z(345)*u1*z(110);
z(368) = 2*e2*z(351) + 2*e3*z(358) - 4*e2p;
z(369) = z(2)*z(356) + z(16)*z(368) - 4*e2*z(61) - 4*rad*z(320)*z(54);
z(370) = e4*z(358) + 2*e2*z(350) + 2*e1p - 2*e3*z(357) - e1*z(351);
z(371) = z(3)*z(356) + z(16)*z(370) + 2*e1*z(61) - 4*rad*z(320)*z(66);
z(372) = 2*e2*z(357) + 2*e3*z(350) + 2*e4p - e1*z(358) - e4*z(351);
z(373) = z(4)*z(356) + z(16)*z(372) + 2*e4*z(61) - 4*rad*z(320)*z(69);
z(374) = z(344)*u1*z(74) - u1*z(369) - u2*z(371) - u3*z(373) - z(345)*u1*z(79);
z(375) = (z(23)*z(345)+4*rad*z(318))*u1;
z(376) = u2*z(348) - u1*z(375) - z(344)*u1*z(109);
z(377) = z(23)*z(344)*u1^2 - u3*z(348) - z(345)*u1*z(109) - rad*u1*z(354);
z(378) = u3*z(375) + z(366) - z(23)*z(344)*u1*u2 - z(344)*u1*z(121) - z(345)*u1*z(120);
z(379) = 0.25*z(346) + 0.25*z(347) + m*(z(4)*z(24)*z(349)+z(4)*z(332)*z(118)+2*e4*z(24)*z(118)+z(15)*z(355)+2*z(24)*z(365)+2*z(332)*z(115)+z(2)*z(24)*z(367)+z(2)*z(332)*z(117)+z(3)*z(24)*z(355)+z(3)*z(332)*z(119)+2*e1*z(24)*z(119)+2*z(6)*z(15)*z(374)-2*e3*z(27)*z(118)-z(7)*z(27)*z(349)-z(7)*z(334)*z(118)-4*rad*z(318)*z(119)-4*e2*z(24)*z(117)-2*e1*z(27)*z(117)-z(5)*z(27)*z(367)-z(5)*z(334)*z(117)-z(6)*z(27)*z(355)-z(6)*z(334)*z(119)-8*rad*z(6)*z(318)*z(116)) + m*(z(4)*z(24)*z(376)+z(4)*z(332)*z(123)+2*e4*z(24)*z(123)+z(15)*z(377)+z(2)*z(24)*z(378)+z(2)*z(332)*z(122)+z(3)*z(24)*z(377)+z(3)*z(332)*z(124)+2*e1*z(24)*z(124)+2*z(3)*z(15)*z(365)+4*e1*z(15)*z(115)-2*e3*z(27)*z(123)-z(7)*z(27)*z(376)-z(7)*z(334)*z(123)-2*z(27)*z(374)-2*z(334)*z(116)-4*rad*z(318)*z(124)-4*e2*z(24)*z(122)-2*e1*z(27)*z(122)-z(5)*z(27)*z(378)-z(5)*z(334)*z(122)-z(6)*z(27)*z(377)-z(6)*z(334)*z(124)-8*rad*z(3)*z(318)*z(115));
z(380) = z(345)*u1*(2*z(1)*u1-z(100));
z(381) = z(345)*u1*(z(1)*u1-z(100));
z(382) = m*(z(4)*z(25)*z(349)+2*e4*z(25)*z(118)+2*z(25)*z(365)+z(2)*z(25)*z(367)+z(3)*z(25)*z(355)+2*e1*z(25)*z(119)+4*rad*z(318)*z(117)+8*e2*z(15)*z(115)+8*rad*z(2)*z(318)*z(115)-2*e3*z(28)*z(118)-z(7)*z(28)*z(349)-z(7)*z(339)*z(118)-4*rad*z(4)*z(6)*z(320)*z(118)-2*z(28)*z(374)-2*z(339)*z(116)-z(15)*z(367)-4*e2*z(25)*z(117)-2*e1*z(28)*z(117)-2*z(2)*z(15)*z(365)-z(5)*z(28)*z(367)-z(5)*z(339)*z(117)-z(6)*z(28)*z(355)-z(6)*z(339)*z(119)-8*rad*z(6)*z(320)*z(115)-4*rad*z(2)*z(6)*z(320)*z(117)-4*rad*z(3)*z(6)*z(320)*z(119)) + m*(z(23)*z(376)+z(4)*z(25)*z(376)+2*e4*z(25)*z(123)+z(23)*z(4)*z(365)+z(23)*z(7)*z(374)+z(2)*z(25)*z(378)+z(3)*z(25)*z(377)+2*z(23)*e3*z(116)+2*z(23)*e4*z(115)+2*e1*z(25)*z(124)+4*rad*z(318)*z(122)+8*rad*z(5)*z(318)*z(116)-2*e3*z(28)*z(123)-z(7)*z(28)*z(376)-z(7)*z(339)*z(123)-4*rad*z(4)*z(6)*z(320)*z(123)-z(15)*z(378)-4*e1*z(15)*z(116)-4*e2*z(25)*z(122)-2*e1*z(28)*z(122)-2*z(5)*z(15)*z(374)-z(5)*z(28)*z(378)-z(5)*z(339)*z(122)-z(6)*z(28)*z(377)-z(6)*z(339)*z(124)-4*rad*z(2)*z(6)*z(320)*z(122)-4*rad*z(3)*z(6)*z(320)*z(124)) - 0.25*z(380) - 0.25*z(381);
z(383) = z(344)*u1*(z(1)*u1-z(100));
z(384) = z(344)*u1*(2*z(1)*u1-z(100));
z(385) = m*(z(4)*z(26)*z(349)+z(4)*z(326)*z(118)+2*e4*z(26)*z(118)+2*z(26)*z(365)+2*z(326)*z(115)+z(2)*z(26)*z(367)+z(2)*z(326)*z(117)+z(3)*z(26)*z(355)+z(3)*z(326)*z(119)+2*e1*z(26)*z(119)-2*e3*z(29)*z(118)-z(7)*z(29)*z(349)-z(7)*z(328)*z(118)-2*z(29)*z(374)-2*z(328)*z(116)-4*e2*z(26)*z(117)-2*e1*z(29)*z(117)-z(5)*z(29)*z(367)-z(5)*z(328)*z(117)-z(6)*z(29)*z(355)-z(6)*z(328)*z(119)) + m*(z(4)*z(26)*z(376)+z(4)*z(326)*z(123)+2*e4*z(26)*z(123)+z(2)*z(26)*z(378)+z(2)*z(326)*z(122)+z(3)*z(26)*z(377)+z(3)*z(326)*z(124)+2*e1*z(26)*z(124)-2*e3*z(29)*z(123)-z(7)*z(29)*z(376)-z(7)*z(328)*z(123)-z(23)*z(377)-4*e2*z(26)*z(122)-2*z(23)*e1*z(115)-2*e1*z(29)*z(122)-z(23)*z(3)*z(365)-z(23)*z(6)*z(374)-z(5)*z(29)*z(378)-z(5)*z(328)*z(122)-z(6)*z(29)*z(377)-z(6)*z(328)*z(124)) - 0.25*z(383) - 0.25*z(384);
z(386) = z(379) + z(51)*z(382) + z(345)*z(134) - z(52)*z(385) - z(344)*z(131);
z(387) = m*(2*z(25)*z(332)+2*z(27)*z(339)+2*z(28)*z(334)+z(23)*z(4)*z(332)+2*z(23)*e4*z(24)+2*z(5)*z(15)*z(334)+4*e1*z(15)*z(25)+4*e1*z(15)*z(27)+8*e2*z(15)*z(24)+8*rad*z(2)*z(24)*z(318)+8*rad*z(6)*z(28)*z(318)-2*z(23)*e3*z(27)-2*z(2)*z(15)*z(332)-2*z(6)*z(15)*z(339)-z(23)*z(7)*z(334)-8*rad*z(3)*z(25)*z(318)-8*rad*z(5)*z(27)*z(318)-8*rad*z(6)*z(24)*z(320)-8*rad*z(3)*z(6)*z(15)*z(320));
z(388) = m*(z(23)*z(7)*z(339)+2*z(23)*e3*z(28)+8*rad*z(15)*z(318)+4*z(211)*z(4)*z(6)*z(320)+8*rad*z(5)*z(28)*z(318)+8*rad*z(6)*z(25)*z(320)-2*z(28)*z(339)-8*e2*z(15)*z(25)-4*e1*z(15)*z(28)-2*z(23)*e4*z(25)-2*z(5)*z(15)*z(339)-8*rad*z(2)*z(25)*z(318)-8*rad*z(2)*z(6)*z(15)*z(320));
z(389) = m*(z(23)*z(7)*z(328)+2*z(23)*e1*z(25)+2*z(23)*e3*z(29)+2*z(2)*z(15)*z(326)+8*rad*z(5)*z(29)*z(318)+8*rad*z(6)*z(26)*z(320)-2*z(25)*z(326)-2*z(28)*z(328)-2*z(29)*z(339)-8*e2*z(15)*z(26)-4*e1*z(15)*z(29)-2*z(23)*e4*z(26)-2*z(5)*z(15)*z(328)-z(23)*z(4)*z(326)-z(23)*z(6)*z(339)-8*rad*z(2)*z(26)*z(318)-4*z(211)*z(3)*z(6)*z(320));
z(390) = z(387) + z(52)*z(389) - 2*z(51)*z(388) - z(129)*z(344) - z(130)*z(345);
z(391) = 2*e1*z(357) + 2*e4*z(350) + 2*e3p - e2*z(358) - e3*z(351);
z(392) = (z(9)*z(17)*z(391)+2*e3*z(17)*z(58)+2*e3*z(9)^2*z(58))/z(17)^1.5;
z(393) = (z(18)*(z(9)*z(392)+z(18)*z(391)+2*e3*z(59))+2*e3*z(9)*(z(18)*z(58)+2*z(9)*z(59))/z(17)^0.5)/z(18)^3;
z(394) = z(15)*z(391) + z(21)*z(352) + rad*z(10)*z(393) + 2*rad*e3*z(57) + 2*rad*z(319)*z(55) - 4*rad*e2*z(60) - 4*rad*z(318)*z(58) - rad*z(9)*z(354);
z(395) = e4*z(351) + 2*e3*z(350) - 2*e2*z(357) - e1*z(358) - 2*e4p;
z(396) = z(23)*z(352) + rad*z(8)*z(354) + 2*rad*e4*z(57) + 4*rad*z(318)*z(63) - z(15)*z(395);
z(397) = 2*rad*e4*z(60) - z(23)*z(391) - z(21)*z(395) - 2*rad*z(319)*z(63) - rad*z(8)*z(393);
z(398) = u1*z(394) + u2*z(396) + u3*z(397) + z(345)*u1*z(90) - z(344)*u1*z(88);
z(399) = rad*(z(18)*z(392)+4*e3*z(9)*z(59)/z(17)^0.5)/z(18)^3;
z(400) = z(360) - 2*e1*z(62) - z(5)*z(399) - z(22)*z(359) - 2*rad*z(321)*z(68);
z(401) = 2*rad*z(21)*z(393) + 2*z(21)*z(36)*z(391) + 4*z(291)*z(319)*z(60) + 2*rad*z(9)*z(36)*z(393) + 2*rad*z(10)*z(36)*z(354) + 4*rad*e3*z(36)*z(60) + 4*rad*z(36)*z(319)*z(58) + 8*rad*e2*z(36)*z(57) + 8*rad*z(36)*z(318)*z(55) + 2*z(36)*(z(356)-z(399)) + 2*z(9)*z(21)*(z(356)-z(399)) + 4*z(23)*e4*(z(61)-z(62)) + 4*z(211)*(z(321)+2*z(320))*z(63) + 4*e3*z(21)*(z(61)-z(62)) + 8*e2*z(15)*(z(61)-z(62)) + 4*rad*z(9)*z(319)*(z(61)-z(62)) + 4*rad*z(15)*(z(321)+2*z(320))*z(55) + 4*z(291)*z(10)*(z(321)+2*z(320))*z(57) + 8*rad*z(10)*z(318)*(z(61)-z(62)) - 8*z(291)*z(318)*z(57) - 2*rad*z(15)*z(354) - 2*z(23)*z(36)*z(395) - 2*z(15)*z(36)*z(352) - 2*z(23)*z(8)*(z(356)-z(399)) - 2*z(10)*z(15)*(z(356)-z(399)) - 4*rad*z(21)*(z(321)+2*z(320))*z(58) - 4*z(291)*z(9)*(z(321)+2*z(320))*z(60) - 4*rad*(z(321)+2*z(320))*(z(61)-z(62));
z(402) = z(23)*(2*z(401)/(z(38)^0.5*z(39)^2)-z(322)*(2+z(39)/z(38)^0.5)*z(64)/(z(38)*z(39)^3));
z(403) = (2*z(322)*(rad*z(39)*z(57)-z(15)*z(64)/z(38)^0.5)/z(38)^0.5+z(39)*(4*rad*z(39)*z(354)+(2*z(15)*z(38)*z(401)-z(15)*z(322)*z(64)-8*rad*z(38)*z(318)*z(64))/z(38)^1.5))/z(39)^3;
z(404) = (2*z(21)*z(322)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(393)+(z(21)*z(322)*z(64)-2*z(21)*z(38)*z(401)-4*rad*z(38)*z(319)*z(64)-2*rad*z(38)*z(322)*z(60))/z(38)^1.5))/z(39)^3;
z(405) = z(2)*z(399) + z(22)*z(368) + 2*rad*z(321)*z(54) - z(369) - 4*e2*z(62);
z(406) = z(15)*(z(2)*z(42)+z(4)*z(40))*z(361) + rad*z(3)*(z(5)*z(42)+z(7)*z(40))*z(354) + z(21)*(z(2)*z(42)-z(3)*z(43))*z(363) + 4*rad*z(318)*(z(5)*z(42)+z(7)*z(40))*z(66) + rad*z(7)*(z(2)*z(42)-z(3)*z(43))*z(393) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(400) + 2*rad*e3*(z(2)*z(42)-z(3)*z(43))*z(60) + 2*rad*z(319)*(z(2)*z(42)-z(3)*z(43))*z(53) + e1*z(15)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) + e3*z(21)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*z(15)*(4*e4*z(40)-8*e2*z(42)-z(2)*z(323)-z(4)*z(324))*z(71) + 2*rad*z(6)*z(318)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) + rad*z(7)*z(319)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*rad*z(6)*(4*e4*z(40)-8*e2*z(42)-z(2)*z(323)-z(4)*z(324))*z(57) + 0.5*z(335)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324))*z(85) + 0.25*z(6)*z(15)*(4*e4*z(70)+4*z(40)*z(372)+4*z(42)*z(368)+8*e2*z(65)-2*z(323)*z(54)-2*z(324)*z(69)-z(2)*z(402)-z(4)*z(403)) + 0.25*z(7)*z(21)*(4*z(42)*z(368)+8*e2*z(65)-4*e1*z(67)-4*z(43)*z(370)-2*z(323)*z(54)-2*z(325)*z(66)-z(2)*z(402)-z(3)*z(404)) + 0.25*z(33)*(4*e4*z(70)+4*z(40)*z(372)+4*z(42)*z(368)+8*e2*z(65)-4*e1*z(67)-4*z(43)*z(370)-2*z(323)*z(54)-2*z(324)*z(69)-2*z(325)*z(66)-z(2)*z(402)-z(3)*z(404)-z(4)*z(403)) - z(15)*(z(5)*z(42)+z(7)*z(40))*z(370) - 4*rad*z(318)*(z(2)*z(42)+z(4)*z(40))*z(71) - 2*rad*e1*(z(5)*z(42)+z(7)*z(40))*z(57) - rad*z(6)*(z(2)*z(42)+z(4)*z(40))*z(354) - z(21)*(z(5)*z(42)-z(6)*z(43))*z(372) - 2*rad*e4*(z(5)*z(42)-z(6)*z(43))*z(60) - 2*rad*z(319)*(z(5)*z(42)-z(6)*z(43))*z(69) - rad*z(4)*(z(5)*z(42)-z(6)*z(43))*z(393) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(405) - 0.5*z(21)*(4*e1*z(42)-z(5)*z(323)-z(6)*z(325))*z(69) - 0.5*rad*z(4)*(4*e1*z(42)-z(5)*z(323)-z(6)*z(325))*z(60) - 0.5*z(21)*(z(2)*z(323)+z(3)*z(325)+4*e1*z(43)+8*e2*z(42))*z(53) - 0.5*rad*z(7)*(z(2)*z(323)+z(3)*z(325)+4*e1*z(43)+8*e2*z(42))*z(60) - 2*rad*z(3)*z(318)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e4*z(21)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(15)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(7)*z(324))*z(66) - rad*z(4)*z(319)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*rad*z(3)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(7)*z(324))*z(57) - 0.5*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324))*z(83) - 0.5*z(333)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.25*z(4)*z(21)*(4*z(42)*z(359)-4*z(43)*z(361)-4*e1*z(65)-2*z(323)*z(68)-2*z(325)*z(71)-z(5)*z(402)-z(6)*z(404)) - 0.25*z(3)*z(15)*(4*e3*z(70)+4*z(40)*z(363)+4*z(42)*z(359)-4*e1*z(65)-2*z(323)*z(68)-2*z(324)*z(53)-z(5)*z(402)-z(7)*z(403)) - 0.25*z(30)*(4*e3*z(70)+4*z(40)*z(363)+4*z(42)*z(359)-4*z(43)*z(361)-4*e1*z(65)-2*z(323)*z(68)-2*z(324)*z(53)-2*z(325)*z(71)-z(5)*z(402)-z(6)*z(404)-z(7)*z(403));
z(407) = z(362) - z(22)*z(361) - z(6)*z(399) - 2*rad*z(321)*z(71);
z(408) = z(3)*z(399) + z(22)*z(370) + 2*e1*z(62) + 2*rad*z(321)*z(66) - z(371);
z(409) = z(23)*(z(2)*z(42)-z(3)*z(43))*z(363) + z(15)*(z(3)*z(43)-z(4)*z(40))*z(359) + rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(354) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(407) + 2*rad*e1*(z(3)*z(43)-z(4)*z(40))*z(57) + 4*rad*e2*(z(6)*z(43)-z(7)*z(40))*z(57) + 4*rad*z(318)*(z(6)*z(43)-z(7)*z(40))*z(54) + 0.5*z(15)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324))*z(54) + 0.5*rad*z(2)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324))*z(57) + z(23)*e3*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 2*rad*z(5)*z(318)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) + z(338)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) + 0.5*z(340)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324))*z(75) + 0.25*z(2)*z(15)*(4*e3*z(70)+4*z(40)*z(363)-4*z(43)*z(361)-2*z(324)*z(53)-2*z(325)*z(71)-z(6)*z(404)-z(7)*z(403)) + 0.25*z(23)*z(7)*(4*z(42)*z(368)+8*e2*z(65)-4*e1*z(67)-4*z(43)*z(370)-2*z(323)*z(54)-2*z(325)*z(66)-z(2)*z(402)-z(3)*z(404)) + 0.25*z(34)*(4*e4*z(70)+4*z(40)*z(372)+4*z(42)*z(368)+8*e2*z(65)-4*e1*z(67)-4*z(43)*z(370)-2*z(323)*z(54)-2*z(324)*z(69)-2*z(325)*z(66)-z(2)*z(402)-z(3)*z(404)-z(4)*z(403)) - z(23)*(z(5)*z(42)-z(6)*z(43))*z(372) - z(15)*(z(6)*z(43)-z(7)*z(40))*z(368) - 4*rad*z(318)*(z(3)*z(43)-z(4)*z(40))*z(68) - rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(354) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(408) - 0.5*z(23)*(4*e1*z(42)-z(5)*z(323)-z(6)*z(325))*z(69) - 0.5*z(23)*(z(2)*z(323)+z(3)*z(325)+4*e1*z(43)+8*e2*z(42))*z(53) - 2*e2*z(15)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - e1*z(15)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 2*rad*z(2)*z(318)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - z(23)*e4*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(15)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324))*z(68) - 0.5*rad*z(5)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324))*z(57) - 0.5*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324))*z(73) - 0.25*z(23)*z(4)*(4*z(42)*z(359)-4*z(43)*z(361)-4*e1*z(65)-2*z(323)*z(68)-2*z(325)*z(71)-z(5)*z(402)-z(6)*z(404)) - 0.25*z(5)*z(15)*(4*e4*z(70)+4*z(40)*z(372)-4*e1*z(67)-4*z(43)*z(370)-2*z(324)*z(69)-2*z(325)*z(66)-z(3)*z(404)-z(4)*z(403)) - 0.25*z(31)*(4*e3*z(70)+4*z(40)*z(363)+4*z(42)*z(359)-4*z(43)*z(361)-4*e1*z(65)-2*z(323)*z(68)-2*z(324)*z(53)-2*z(325)*z(71)-z(5)*z(402)-z(6)*z(404)-z(7)*z(403));
z(410) = z(364) - 2*e3*z(62) - z(7)*z(399) - z(22)*z(363) - 2*rad*z(321)*z(53);
z(411) = z(4)*z(399) + z(22)*z(372) + 2*e4*z(62) + 2*rad*z(321)*z(69) - z(373);
z(412) = z(23)*(z(5)*z(42)+z(7)*z(40))*z(370) + z(21)*(z(3)*z(43)-z(4)*z(40))*z(359) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(393) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(410) + 2*rad*e1*(z(3)*z(43)-z(4)*z(40))*z(60) + 2*rad*z(319)*(z(3)*z(43)-z(4)*z(40))*z(68) + 4*rad*e2*(z(6)*z(43)-z(7)*z(40))*z(60) + 0.5*z(21)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324))*z(54) + 0.5*rad*z(2)*(4*e3*z(40)-z(6)*z(325)-z(7)*z(324))*z(60) + rad*z(2)*z(319)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(23)*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(7)*z(324))*z(66) + 0.5*z(329)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e4*z(40)-8*e2*z(42)-4*e1*z(43)-z(2)*z(323)-z(3)*z(325)-z(4)*z(324))*z(80) + 0.25*z(2)*z(21)*(4*e3*z(70)+4*z(40)*z(363)-4*z(43)*z(361)-2*z(324)*z(53)-2*z(325)*z(71)-z(6)*z(404)-z(7)*z(403)) + 0.25*z(23)*z(3)*(4*e3*z(70)+4*z(40)*z(363)+4*z(42)*z(359)-4*e1*z(65)-2*z(323)*z(68)-2*z(324)*z(53)-z(5)*z(402)-z(7)*z(403)) + 0.25*z(35)*(4*e4*z(70)+4*z(40)*z(372)+4*z(42)*z(368)+8*e2*z(65)-4*e1*z(67)-4*z(43)*z(370)-2*z(323)*z(54)-2*z(324)*z(69)-2*z(325)*z(66)-z(2)*z(402)-z(3)*z(404)-z(4)*z(403)) - z(23)*(z(2)*z(42)+z(4)*z(40))*z(361) - z(21)*(z(6)*z(43)-z(7)*z(40))*z(368) - 2*rad*z(319)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(393) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(411) - 2*e2*z(21)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - z(23)*e1*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e1*z(21)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - rad*z(5)*z(319)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(23)*(4*e4*z(40)-8*e2*z(42)-z(2)*z(323)-z(4)*z(324))*z(71) - 0.5*z(21)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324))*z(68) - 0.5*rad*z(5)*(4*e4*z(40)-4*e1*z(43)-z(3)*z(325)-z(4)*z(324))*z(60) - 0.5*(4*e1*z(42)+4*e3*z(40)-z(5)*z(323)-z(6)*z(325)-z(7)*z(324))*z(78) - 0.5*z(327)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.25*z(23)*z(6)*(4*e4*z(70)+4*z(40)*z(372)+4*z(42)*z(368)+8*e2*z(65)-2*z(323)*z(54)-2*z(324)*z(69)-z(2)*z(402)-z(4)*z(403)) - 0.25*z(5)*z(21)*(4*e4*z(70)+4*z(40)*z(372)-4*e1*z(67)-4*z(43)*z(370)-2*z(324)*z(69)-2*z(325)*z(66)-z(3)*z(404)-z(4)*z(403)) - 0.25*z(32)*(4*e3*z(70)+4*z(40)*z(363)+4*z(42)*z(359)-4*z(43)*z(361)-4*e1*z(65)-2*z(323)*z(68)-2*z(324)*z(53)-2*z(325)*z(71)-z(5)*z(402)-z(6)*z(404)-z(7)*z(403));
z(413) = u1*z(406) + u2*z(409) + u3*z(412) + z(345)*u1*z(81) - z(344)*u1*z(76);
z(414) = (z(343)*(z(45)*z(91)-z(49)*z(87))-z(50)*(z(45)*z(398)+z(330)*z(91)-z(49)*z(413)-z(337)*z(87)))/z(50)^2;
z(415) = m*(z(23)*z(3)*z(332)+2*z(23)*e1*z(24)+2*z(6)*z(15)*z(328)+8*rad*z(3)*z(26)*z(318)-4*z(211)*z(318)-2*z(24)*z(326)-2*z(26)*z(332)-2*z(27)*z(328)-2*z(29)*z(334)-4*e1*z(15)*z(26)-2*z(3)*z(15)*z(326)-z(23)*z(6)*z(334)-8*rad*z(6)*z(29)*z(318));
z(416) = 4*m*(z(26)*z(326)+z(29)*z(328)) - 2*z(132)*(z(3)*z(326)+2*e1*z(26)-z(6)*z(328));
z(417) = z(130)*z(344) + z(133)*z(345) - z(415) - z(51)*z(389) - z(52)*z(416);
z(418) = (z(343)*(z(44)*z(91)-z(47)*z(87))-z(50)*(z(44)*z(398)+z(341)*z(91)-z(47)*z(413)-z(342)*z(87)))/z(50)^2;
z(419) = z(386) + z(137)*z(418) + z(390)*z(92) - z(136)*z(414) - z(417)*z(93);
z(420) = -z(15)*z(344) - 4*rad*z(51)*z(318);
z(421) = -4*rad*z(318) - z(23)*z(345);
z(422) = m*(z(9)*z(421)+2*e3*z(15)+2*e3*z(97)+4*e4*z(96)-4*e2*z(98)-2*z(8)*z(420)-4*rad*z(9)*z(318)-z(23)*z(10)*z(344));
z(423) = m*(z(6)*z(15)*z(334)+4*rad*z(15)*z(318)+4*rad*z(3)*z(24)*z(318)-z(24)*z(332)-z(27)*z(334)-2*e1*z(15)*z(24)-z(3)*z(15)*z(332)-4*rad*z(6)*z(27)*z(318));
z(424) = z(51)*z(387) + z(52)*z(415) + z(127)*z(345) - 4*z(423) - z(126)*z(344);
z(425) = z(424) + z(51)*z(390) + z(137)*z(345) - z(52)*z(417) - z(136)*z(344);
z(426) = (z(425)*z(141)-z(139)*(9.810000000000001*z(422)+z(419)))/z(139)^2;
a13 = z(426);
z(427) = e2*(z(18)+z(9)^2/z(17)^0.5)/z(18)^2;
z(428) = e2*z(9)/(z(17)^0.5*z(18)^2);
z(429) = 4*rad*z(21)*z(427) + 4*z(211)*z(8)*z(428) + 4*e2*z(21)*z(36) + 4*rad*z(9)*z(36)*z(427) + 4*rad*z(10)*z(15)*z(428) - 4*rad*z(36)*z(428) - 4*z(23)*e1*z(36) - 4*rad*z(9)*z(21)*z(428);
z(430) = z(23)*z(429)/(z(38)^0.5*z(39)^2);
z(431) = z(15)*z(429)/(z(38)^0.5*z(39)^2);
z(432) = (4*rad*z(39)*z(427)-z(21)*z(429)/z(38)^0.5)/z(39)^2;
z(433) = 2*e2*z(16) - 2*e2*z(22) - 2*rad*z(7)*z(428);
z(434) = 2*e1*z(22) + 2*rad*z(4)*z(428) - 2*e1*z(16);
z(435) = 4*z(23)*e3*(z(2)*z(42)+z(4)*z(40)) + 2*e4*z(21)*(z(3)*z(43)-z(4)*z(40)) + 4*e3*z(21)*(z(6)*z(43)-z(7)*z(40)) + 2*rad*z(5)*z(427)*(z(3)*z(43)-z(4)*z(40)) + z(433)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(23)*z(3)*(4*e2*z(40)+4*e4*z(42)-z(5)*z(430)-z(7)*z(431)) + 0.5*z(2)*z(21)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431)) + 0.5*z(35)*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431)) - 2*z(23)*e4*(z(5)*z(42)+z(7)*z(40)) - 2*rad*z(2)*z(427)*(z(6)*z(43)-z(7)*z(40)) - z(434)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(6)*(4*e1*z(40)-8*e3*z(42)-z(2)*z(430)-z(4)*z(431)) - 0.5*z(5)*z(21)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431)) - 0.5*z(32)*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431));
z(436) = 2*e2*z(15) + 2*rad*z(10)*z(427);
z(437) = 2*e4*z(16) - 2*e4*z(22) - 2*rad*z(5)*z(428);
z(438) = 4*e3*z(16) + 2*rad*z(2)*z(428) - 4*e3*z(22);
z(439) = 2*e4*z(15)*(z(5)*z(42)+z(7)*z(40)) + 2*e2*z(21)*(z(2)*z(42)-z(3)*z(43)) + 2*rad*z(7)*z(427)*(z(2)*z(42)-z(3)*z(43)) + z(437)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(6)*z(15)*(4*e1*z(40)-8*e3*z(42)-z(2)*z(430)-z(4)*z(431)) + 0.5*z(7)*z(21)*(4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)) + 0.5*z(33)*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431)) - 4*e3*z(15)*(z(2)*z(42)+z(4)*z(40)) - 2*e1*z(21)*(z(5)*z(42)-z(6)*z(43)) - 2*rad*z(4)*z(427)*(z(5)*z(42)-z(6)*z(43)) - z(438)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(3)*z(15)*(4*e2*z(40)+4*e4*z(42)-z(5)*z(430)-z(7)*z(431)) - 0.5*z(4)*z(21)*(4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)) - 0.5*z(30)*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431));
z(440) = -2*z(23)*e2 - 2*e1*z(21) - 2*rad*z(8)*z(427);
z(441) = 4*e3*z(22) - 4*e3*z(16) - 2*rad*z(6)*z(428);
z(442) = 2*e4*z(16) + 2*rad*z(3)*z(428) - 2*e4*z(22);
z(443) = 2*z(23)*e2*(z(2)*z(42)-z(3)*z(43)) + 2*e4*z(15)*(z(3)*z(43)-z(4)*z(40)) + 4*e3*z(15)*(z(6)*z(43)-z(7)*z(40)) + z(441)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(23)*z(7)*(4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)) + 0.5*z(2)*z(15)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431)) + 0.5*z(34)*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431)) - 2*z(23)*e1*(z(5)*z(42)-z(6)*z(43)) - z(442)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(23)*z(4)*(4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)) - 0.5*z(5)*z(15)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431)) - 0.5*z(31)*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431));
z(444) = z(44)*z(440) + z(49)*z(443) + 2*e1*z(15)*z(45) - z(47)*z(435);
z(445) = (z(444)*(z(45)*z(48)-z(46)*z(49))-z(50)*(z(45)*z(436)+z(48)*z(435)-z(46)*z(440)-z(49)*z(439)))/z(50)^2;
z(446) = (z(444)*(z(44)*z(48)-z(46)*z(47))+z(50)*(z(47)*z(439)-z(44)*z(436)-z(48)*z(443)-2*e1*z(15)*z(46)))/z(50)^2;
z(447) = u1*(z(1)*z(445)*u3+2*z(1)*z(446)*u2-2*z(445)*z(102)-z(446)*z(101));
z(448) = u1*(z(1)*z(446)*u2+2*z(1)*z(445)*u3-2*z(446)*z(101)-z(445)*z(102));
z(449) = z(445)*u1*(z(15)*u2+z(109));
z(450) = 0.5*e2*z(446)*u1 + 0.5*e3*z(445)*u1 - 0.5*u2;
z(451) = (1-e1*z(446)-e4*z(445))*u1;
z(452) = -4*e1*z(450) - 2*e2*z(451);
z(453) = (z(12)+z(10)^2/z(11)^0.5)*z(452)/z(12)^2;
z(454) = u1*(z(15)*z(445)*u3+rad*z(453)-z(446)*z(109));
z(455) = rad*z(10)*z(452)/(z(11)^0.5*z(12)^2);
z(456) = 0.5*e2*z(445)*u1 - 0.5*u3 - 0.5*e3*z(446)*u1;
z(457) = (e1*z(445)-e4*z(446))*u1;
z(458) = e1*z(451) + 2*e2*z(450) + 2*e3*z(456) + 2*e4p - e4*z(457);
z(459) = z(5)*z(455) + z(16)*z(458) + 2*e4*z(61);
z(460) = 2*e3*z(457) - 4*e1*z(450) - 4*e3p;
z(461) = z(6)*z(455) + z(16)*z(460) - 4*e3*z(61);
z(462) = e3*z(451) + 2*e2p - 2*e1*z(456) - 2*e4*z(450) - e2*z(457);
z(463) = z(7)*z(455) + z(16)*z(462) + 2*e2*z(61);
z(464) = u1*z(459) + u2*z(461) + u3*z(463) + z(446)*u1*z(77) - z(445)*u1*z(72);
z(465) = rad*(u2*z(453)-z(445)*u1*z(57));
z(466) = -z(446)*u1*z(110) - z(465);
z(467) = 2*e3*z(457) - 2*e2*z(451) - 4*e3p;
z(468) = z(2)*z(455) + z(16)*z(467) - 4*e3*z(61);
z(469) = e1*z(451) + e4*z(457) + 2*e2*z(450) - 2*e3*z(456) - 2*e4p;
z(470) = z(3)*z(455) + z(16)*z(469) - 2*e4*z(61);
z(471) = e4*z(451) + 2*e2*z(456) + 2*e3*z(450) + 2*e1p - e1*z(457);
z(472) = z(4)*z(455) + z(16)*z(471) + 2*e1*z(61);
z(473) = z(445)*u1*z(74) - u1*z(468) - u2*z(470) - u3*z(472) - z(446)*u1*z(79);
z(474) = u1*(z(23)*z(446)*u1+z(15)*z(445)*u2+z(445)*z(109));
z(475) = u1*(z(23)*z(445)*u1+z(15)*z(445)*u3+rad*z(453)-z(446)*z(109));
z(476) = z(23)*z(446)*u1*u3 - z(23)*z(445)*u1*u2 - z(445)*u1*z(121) - z(446)*u1*z(120) - z(465);
z(477) = 0.25*z(447) + 0.25*z(448) - m*(z(4)*z(24)*z(449)+z(6)*z(27)*z(454)+2*e2*z(27)*z(118)+z(5)*z(27)*z(466)+2*e4*z(24)*z(119)+2*e4*z(27)*z(117)+4*e3*z(24)*z(117)+8*e3*z(15)*z(116)-z(15)*z(454)-2*e1*z(24)*z(118)-z(3)*z(24)*z(454)-z(7)*z(27)*z(449)-4*e3*z(7)*z(16)*z(118)-2*e4*z(4)*z(16)*z(118)-2*z(24)*z(464)-4*e3*z(27)*z(119)-4*e4*z(16)*z(115)-2*z(6)*z(15)*z(473)-z(2)*z(24)*z(466)-4*e3*z(5)*z(16)*z(117)-4*e3*z(6)*z(16)*z(119)-2*e4*z(2)*z(16)*z(117)-2*e4*z(3)*z(16)*z(119)) - m*(z(4)*z(24)*z(474)+z(6)*z(27)*z(475)+2*e2*z(27)*z(123)+2*z(27)*z(473)+z(5)*z(27)*z(476)+2*e4*z(24)*z(124)+2*e4*z(27)*z(122)+4*e3*z(24)*z(122)+4*e4*z(15)*z(115)-z(15)*z(475)-2*e1*z(24)*z(123)-z(3)*z(24)*z(475)-z(7)*z(27)*z(474)-4*e3*z(7)*z(16)*z(123)-2*e4*z(4)*z(16)*z(123)-8*e3*z(16)*z(116)-4*e3*z(27)*z(124)-2*z(3)*z(15)*z(464)-z(2)*z(24)*z(476)-4*e3*z(5)*z(16)*z(122)-4*e3*z(6)*z(16)*z(124)-2*e4*z(2)*z(16)*z(122)-2*e4*z(3)*z(16)*z(124));
z(478) = z(446)*u1*(2*z(1)*u1-z(100));
z(479) = z(446)*u1*(z(1)*u1-z(100));
z(480) = -0.25*z(478) - 0.25*z(479) - m*(z(4)*z(25)*z(449)+z(6)*z(28)*z(454)+2*e2*z(28)*z(118)+4*e3*z(4)*z(16)*z(118)+z(15)*z(466)+2*z(28)*z(473)+z(5)*z(28)*z(466)+2*e4*z(25)*z(119)+2*e4*z(28)*z(117)+2*z(2)*z(15)*z(464)+4*e3*z(25)*z(117)+8*e3*z(16)*z(115)+4*e3*z(2)*z(16)*z(117)+4*e3*z(3)*z(16)*z(119)-2*e1*z(25)*z(118)-z(3)*z(25)*z(454)-z(7)*z(28)*z(449)-2*e4*z(7)*z(16)*z(118)-2*z(25)*z(464)-8*e3*z(15)*z(115)-4*e3*z(28)*z(119)-4*e4*z(16)*z(116)-z(2)*z(25)*z(466)-2*e4*z(5)*z(16)*z(117)-2*e4*z(6)*z(16)*z(119)) - m*(z(23)*z(474)+z(4)*z(25)*z(474)+z(6)*z(28)*z(475)+2*e2*z(28)*z(123)+4*e3*z(4)*z(16)*z(123)+z(15)*z(476)+z(5)*z(28)*z(476)+2*e4*z(25)*z(124)+2*e4*z(28)*z(122)+2*z(5)*z(15)*z(473)+4*e3*z(25)*z(122)+4*e4*z(15)*z(116)+4*e3*z(2)*z(16)*z(122)+4*e3*z(3)*z(16)*z(124)-2*e1*z(25)*z(123)-z(3)*z(25)*z(475)-z(7)*z(28)*z(474)-2*e4*z(7)*z(16)*z(123)-4*e3*z(28)*z(124)-2*z(23)*e1*z(115)-2*z(23)*e2*z(116)-z(23)*z(4)*z(464)-z(23)*z(7)*z(473)-z(2)*z(25)*z(476)-2*e4*z(5)*z(16)*z(122)-2*e4*z(6)*z(16)*z(124));
z(481) = z(445)*u1*(z(1)*u1-z(100));
z(482) = z(445)*u1*(2*z(1)*u1-z(100));
z(483) = -0.25*z(481) - 0.25*z(482) - m*(z(4)*z(26)*z(449)+z(6)*z(29)*z(454)+2*e2*z(29)*z(118)+2*e1*z(7)*z(16)*z(118)+2*z(29)*z(473)+z(5)*z(29)*z(466)+2*e4*z(26)*z(119)+2*e4*z(29)*z(117)+4*e1*z(16)*z(116)+4*e3*z(26)*z(117)+2*e1*z(5)*z(16)*z(117)+2*e1*z(6)*z(16)*z(119)-2*e1*z(26)*z(118)-z(3)*z(26)*z(454)-z(7)*z(29)*z(449)-2*e2*z(4)*z(16)*z(118)-2*z(26)*z(464)-4*e2*z(16)*z(115)-4*e3*z(29)*z(119)-z(2)*z(26)*z(466)-2*e2*z(2)*z(16)*z(117)-2*e2*z(3)*z(16)*z(119)) - m*(z(23)*z(475)+z(4)*z(26)*z(474)+z(6)*z(29)*z(475)+2*e2*z(29)*z(123)+2*e1*z(7)*z(16)*z(123)+z(23)*z(3)*z(464)+z(23)*z(6)*z(473)+z(5)*z(29)*z(476)+2*e4*z(26)*z(124)+2*e4*z(29)*z(122)+4*e3*z(26)*z(122)+2*e1*z(5)*z(16)*z(122)+2*e1*z(6)*z(16)*z(124)-2*e1*z(26)*z(123)-z(3)*z(26)*z(475)-z(7)*z(29)*z(474)-2*e2*z(4)*z(16)*z(123)-4*z(23)*e3*z(116)-4*e3*z(29)*z(124)-2*z(23)*e4*z(115)-z(2)*z(26)*z(476)-2*e2*z(2)*z(16)*z(122)-2*e2*z(3)*z(16)*z(124));
z(484) = z(477) + z(51)*z(480) + z(446)*z(134) - z(52)*z(483) - z(445)*z(131);
z(485) = m*(z(23)*e2*z(27)+2*e4*z(15)*z(25)+2*e4*z(16)*z(27)+4*e3*z(16)*z(24)+4*e3*z(16)*z(28)+2*e4*z(2)*z(15)*z(16)+4*e3*z(3)*z(15)*z(16)+4*e3*z(5)*z(15)*z(16)-4*e3*z(15)*z(24)-4*e3*z(15)*z(28)-2*e4*z(15)*z(27)-2*e4*z(16)*z(25)-z(23)*e1*z(24)-2*z(23)*e3*z(7)*z(16)-2*e4*z(6)*z(15)*z(16)-z(23)*e4*z(4)*z(16));
z(486) = m*(z(23)*e1*z(25)+2*e4*z(15)*z(28)+4*e3*z(15)*z(25)+z(23)*e4*z(7)*z(16)+4*e3*z(2)*z(15)*z(16)-4*e3*z(16)*z(25)-2*e4*z(16)*z(28)-z(23)*e2*z(28)-2*z(23)*e3*z(4)*z(16)-2*e4*z(5)*z(15)*z(16));
z(487) = m*(z(23)*e1*z(26)+z(23)*e4*z(25)+2*e1*z(16)*z(28)+2*e2*z(16)*z(25)+2*e4*z(15)*z(29)+4*e3*z(15)*z(26)+z(23)*e2*z(4)*z(16)+2*z(23)*e3*z(3)*z(16)+2*e1*z(5)*z(15)*z(16)-4*e3*z(16)*z(26)-2*z(23)*e3*z(28)-2*e4*z(16)*z(29)-z(23)*e2*z(29)-2*e2*z(2)*z(15)*z(16)-z(23)*e1*z(7)*z(16)-z(23)*e4*z(6)*z(16));
z(488) = 4*z(51)*z(486) - 2*z(485) - 2*z(52)*z(487) - z(129)*z(445) - z(130)*z(446);
z(489) = e3*z(451) + 2*e1*z(456) + 2*e4*z(450) + 2*e2p - e2*z(457);
z(490) = (z(9)*z(17)*z(489)+2*e2*z(17)*z(58)+2*e2*z(9)^2*z(58))/z(17)^1.5;
z(491) = (z(18)*(z(9)*z(490)+z(18)*z(489)+2*e2*z(59))+2*e2*z(9)*(z(18)*z(58)+2*z(9)*z(59))/z(17)^0.5)/z(18)^3;
z(492) = z(21)*z(452) + rad*z(9)*z(453) + z(15)*z(489) + rad*z(10)*z(491) + 2*rad*e2*z(57) + 2*rad*z(427)*z(55);
z(493) = 2*e3*z(450) + 2*e1p - 2*e2*z(456) - e1*z(457) - e4*z(451);
z(494) = z(23)*z(452) - rad*z(8)*z(453) - z(15)*z(493) - 2*rad*e1*z(57);
z(495) = -z(23)*z(489) - z(21)*z(493) - 2*rad*e1*z(60) - 2*rad*z(427)*z(63) - rad*z(8)*z(491);
z(496) = u1*z(492) + u2*z(494) + u3*z(495) + z(446)*u1*z(90) - z(445)*u1*z(88);
z(497) = rad*(z(18)*z(490)+4*e2*z(9)*z(59)/z(17)^0.5)/z(18)^3;
z(498) = z(459) - 2*e4*z(62) - z(5)*z(497) - z(22)*z(458) - 2*rad*z(428)*z(68);
z(499) = 2*rad*z(15)*z(453) + 2*rad*z(21)*z(491) + 2*z(21)*z(36)*z(489) + 4*z(211)*z(428)*z(63) + 4*z(291)*z(427)*z(60) + 2*rad*z(9)*z(36)*z(491) + 4*rad*e2*z(36)*z(60) + 4*rad*z(15)*z(428)*z(55) + 4*rad*z(36)*z(427)*z(58) + 4*z(291)*z(10)*z(428)*z(57) + 2*z(36)*(z(455)-z(497)) + 2*z(9)*z(21)*(z(455)-z(497)) + 4*e2*z(21)*(z(61)-z(62)) + 4*rad*z(9)*z(427)*(z(61)-z(62)) - 2*z(15)*z(36)*z(452) - 2*rad*z(10)*z(36)*z(453) - 2*z(23)*z(36)*z(493) - 4*rad*z(21)*z(428)*z(58) - 4*z(291)*z(9)*z(428)*z(60) - 4*rad*z(428)*(z(61)-z(62)) - 4*z(23)*e1*(z(61)-z(62)) - 2*z(23)*z(8)*(z(455)-z(497)) - 2*z(10)*z(15)*(z(455)-z(497));
z(500) = z(23)*(2*z(499)/(z(38)^0.5*z(39)^2)-z(429)*(2+z(39)/z(38)^0.5)*z(64)/(z(38)*z(39)^3));
z(501) = (2*z(15)*z(429)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(453)+(z(15)*z(429)*z(64)-2*z(15)*z(38)*z(499)-2*rad*z(38)*z(429)*z(57))/z(38)^1.5))/z(39)^3;
z(502) = (2*z(21)*z(429)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(491)+(z(21)*z(429)*z(64)-2*z(21)*z(38)*z(499)-4*rad*z(38)*z(427)*z(64)-2*rad*z(38)*z(429)*z(60))/z(38)^1.5))/z(39)^3;
z(503) = z(2)*z(497) + z(22)*z(467) + 2*rad*z(428)*z(54) - z(468) - 4*e3*z(62);
z(504) = rad*z(6)*(z(2)*z(42)+z(4)*z(40))*z(453) + z(15)*(z(2)*z(42)+z(4)*z(40))*z(460) + z(21)*(z(2)*z(42)-z(3)*z(43))*z(462) + 2*rad*e4*(z(5)*z(42)+z(7)*z(40))*z(57) + rad*z(7)*(z(2)*z(42)-z(3)*z(43))*z(491) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(498) + 2*rad*e2*(z(2)*z(42)-z(3)*z(43))*z(60) + 2*rad*z(427)*(z(2)*z(42)-z(3)*z(43))*z(53) + 2*e3*z(15)*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) + e2*z(21)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*z(15)*(4*e1*z(40)-8*e3*z(42)-z(2)*z(430)-z(4)*z(431))*z(71) + 0.5*z(21)*(4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432))*z(53) + rad*z(7)*z(427)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*rad*z(6)*(4*e1*z(40)-8*e3*z(42)-z(2)*z(430)-z(4)*z(431))*z(57) + 0.5*rad*z(7)*(4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432))*z(60) + 0.5*z(438)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431))*z(85) + 0.25*z(6)*z(15)*(z(4)*z(501)+4*e1*z(70)+4*z(40)*z(471)+4*z(42)*z(467)+8*e3*z(65)-2*z(430)*z(54)-2*z(431)*z(69)-z(2)*z(500)) + 0.25*z(7)*z(21)*(4*e4*z(67)+4*z(42)*z(467)+8*e3*z(65)-4*z(43)*z(469)-2*z(430)*z(54)-2*z(432)*z(66)-z(2)*z(500)-z(3)*z(502)) + 0.25*z(33)*(z(4)*z(501)+4*e1*z(70)+4*e4*z(67)+4*z(40)*z(471)+4*z(42)*z(467)+8*e3*z(65)-4*z(43)*z(469)-2*z(430)*z(54)-2*z(431)*z(69)-2*z(432)*z(66)-z(2)*z(500)-z(3)*z(502)) - rad*z(3)*(z(5)*z(42)+z(7)*z(40))*z(453) - z(15)*(z(5)*z(42)+z(7)*z(40))*z(469) - 4*rad*e3*(z(2)*z(42)+z(4)*z(40))*z(57) - z(21)*(z(5)*z(42)-z(6)*z(43))*z(471) - 2*rad*e1*(z(5)*z(42)-z(6)*z(43))*z(60) - 2*rad*z(427)*(z(5)*z(42)-z(6)*z(43))*z(69) - rad*z(4)*(z(5)*z(42)-z(6)*z(43))*z(491) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(503) - e4*z(15)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e1*z(21)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(15)*(4*e2*z(40)+4*e4*z(42)-z(5)*z(430)-z(7)*z(431))*z(66) - 0.5*z(21)*(4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432))*z(69) - rad*z(4)*z(427)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*rad*z(3)*(4*e2*z(40)+4*e4*z(42)-z(5)*z(430)-z(7)*z(431))*z(57) - 0.5*rad*z(4)*(4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432))*z(60) - 0.5*z(437)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431))*z(83) - 0.25*z(3)*z(15)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+4*z(42)*z(458)-4*e4*z(65)-2*z(430)*z(68)-2*z(431)*z(53)-z(5)*z(500)) - 0.25*z(4)*z(21)*(4*z(42)*z(458)+8*e3*z(67)-4*e4*z(65)-4*z(43)*z(460)-2*z(430)*z(68)-2*z(432)*z(71)-z(5)*z(500)-z(6)*z(502)) - 0.25*z(30)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+4*z(42)*z(458)+8*e3*z(67)-4*e4*z(65)-4*z(43)*z(460)-2*z(430)*z(68)-2*z(431)*z(53)-2*z(432)*z(71)-z(5)*z(500)-z(6)*z(502));
z(505) = z(461) + 4*e3*z(62) - z(6)*z(497) - z(22)*z(460) - 2*rad*z(428)*z(71);
z(506) = z(3)*z(497) + z(22)*z(469) + 2*rad*z(428)*z(66) - z(470) - 2*e4*z(62);
z(507) = rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(453) + z(23)*(z(2)*z(42)-z(3)*z(43))*z(462) + z(15)*(z(3)*z(43)-z(4)*z(40))*z(458) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(505) + 2*rad*e4*(z(3)*z(43)-z(4)*z(40))*z(57) + 4*rad*e3*(z(6)*z(43)-z(7)*z(40))*z(57) + z(23)*e2*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*z(23)*(4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432))*z(53) + 0.5*z(15)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431))*z(54) + 0.5*rad*z(2)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431))*z(57) + 0.5*z(442)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431))*z(75) + 0.25*z(2)*z(15)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+8*e3*z(67)-4*z(43)*z(460)-2*z(431)*z(53)-2*z(432)*z(71)-z(6)*z(502)) + 0.25*z(23)*z(7)*(4*e4*z(67)+4*z(42)*z(467)+8*e3*z(65)-4*z(43)*z(469)-2*z(430)*z(54)-2*z(432)*z(66)-z(2)*z(500)-z(3)*z(502)) + 0.25*z(34)*(z(4)*z(501)+4*e1*z(70)+4*e4*z(67)+4*z(40)*z(471)+4*z(42)*z(467)+8*e3*z(65)-4*z(43)*z(469)-2*z(430)*z(54)-2*z(431)*z(69)-2*z(432)*z(66)-z(2)*z(500)-z(3)*z(502)) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(453) - z(23)*(z(5)*z(42)-z(6)*z(43))*z(471) - z(15)*(z(6)*z(43)-z(7)*z(40))*z(467) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(506) - 2*e3*z(15)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - e4*z(15)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - z(23)*e1*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(23)*(4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432))*z(69) - 0.5*z(15)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431))*z(68) - 0.5*rad*z(5)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431))*z(57) - 0.5*z(441)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431))*z(73) - 0.25*z(5)*z(15)*(z(4)*z(501)+4*e1*z(70)+4*e4*z(67)+4*z(40)*z(471)-4*z(43)*z(469)-2*z(431)*z(69)-2*z(432)*z(66)-z(3)*z(502)) - 0.25*z(23)*z(4)*(4*z(42)*z(458)+8*e3*z(67)-4*e4*z(65)-4*z(43)*z(460)-2*z(430)*z(68)-2*z(432)*z(71)-z(5)*z(500)-z(6)*z(502)) - 0.25*z(31)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+4*z(42)*z(458)+8*e3*z(67)-4*e4*z(65)-4*z(43)*z(460)-2*z(430)*z(68)-2*z(431)*z(53)-2*z(432)*z(71)-z(5)*z(500)-z(6)*z(502));
z(508) = z(463) - 2*e2*z(62) - z(7)*z(497) - z(22)*z(462) - 2*rad*z(428)*z(53);
z(509) = z(4)*z(497) + z(22)*z(471) + 2*e1*z(62) + 2*rad*z(428)*z(69) - z(472);
z(510) = z(23)*(z(5)*z(42)+z(7)*z(40))*z(469) + z(21)*(z(3)*z(43)-z(4)*z(40))*z(458) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(491) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(508) + 2*rad*e4*(z(3)*z(43)-z(4)*z(40))*z(60) + 2*rad*z(427)*(z(3)*z(43)-z(4)*z(40))*z(68) + 4*rad*e3*(z(6)*z(43)-z(7)*z(40))*z(60) + z(23)*e4*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) + rad*z(2)*z(427)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(23)*(4*e2*z(40)+4*e4*z(42)-z(5)*z(430)-z(7)*z(431))*z(66) + 0.5*z(21)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431))*z(54) + 0.5*rad*z(2)*(4*e2*z(40)+8*e3*z(43)-z(6)*z(432)-z(7)*z(431))*z(60) + 0.5*z(434)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e1*z(40)+4*e4*z(43)-8*e3*z(42)-z(2)*z(430)-z(3)*z(432)-z(4)*z(431))*z(80) + 0.25*z(23)*z(3)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+4*z(42)*z(458)-4*e4*z(65)-2*z(430)*z(68)-2*z(431)*z(53)-z(5)*z(500)) + 0.25*z(2)*z(21)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+8*e3*z(67)-4*z(43)*z(460)-2*z(431)*z(53)-2*z(432)*z(71)-z(6)*z(502)) + 0.25*z(35)*(z(4)*z(501)+4*e1*z(70)+4*e4*z(67)+4*z(40)*z(471)+4*z(42)*z(467)+8*e3*z(65)-4*z(43)*z(469)-2*z(430)*z(54)-2*z(431)*z(69)-2*z(432)*z(66)-z(2)*z(500)-z(3)*z(502)) - z(23)*(z(2)*z(42)+z(4)*z(40))*z(460) - z(21)*(z(6)*z(43)-z(7)*z(40))*z(467) - 2*rad*z(427)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(491) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(509) - 2*z(23)*e3*(z(2)*z(65)-2*z(40)*z(69)-2*z(42)*z(54)-z(4)*z(70)) - 2*e3*z(21)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) - e4*z(21)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - rad*z(5)*z(427)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(23)*(4*e1*z(40)-8*e3*z(42)-z(2)*z(430)-z(4)*z(431))*z(71) - 0.5*z(21)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431))*z(68) - 0.5*rad*z(5)*(4*e1*z(40)+4*e4*z(43)-z(3)*z(432)-z(4)*z(431))*z(60) - 0.5*z(433)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*(4*e2*z(40)+4*e4*z(42)+8*e3*z(43)-z(5)*z(430)-z(6)*z(432)-z(7)*z(431))*z(78) - 0.25*z(23)*z(6)*(z(4)*z(501)+4*e1*z(70)+4*z(40)*z(471)+4*z(42)*z(467)+8*e3*z(65)-2*z(430)*z(54)-2*z(431)*z(69)-z(2)*z(500)) - 0.25*z(5)*z(21)*(z(4)*z(501)+4*e1*z(70)+4*e4*z(67)+4*z(40)*z(471)-4*z(43)*z(469)-2*z(431)*z(69)-2*z(432)*z(66)-z(3)*z(502)) - 0.25*z(32)*(z(7)*z(501)+4*e2*z(70)+4*z(40)*z(462)+4*z(42)*z(458)+8*e3*z(67)-4*e4*z(65)-4*z(43)*z(460)-2*z(430)*z(68)-2*z(431)*z(53)-2*z(432)*z(71)-z(5)*z(500)-z(6)*z(502));
z(511) = u1*z(504) + u2*z(507) + u3*z(510) + z(446)*u1*z(81) - z(445)*u1*z(76);
z(512) = (z(444)*(z(45)*z(91)-z(49)*z(87))-z(50)*(z(45)*z(496)+z(435)*z(91)-z(49)*z(511)-z(440)*z(87)))/z(50)^2;
z(513) = m*(2*z(23)*e3*z(27)+2*e4*z(15)*z(26)+4*e3*z(16)*z(29)+z(23)*e4*z(3)*z(16)+2*z(23)*e3*z(6)*z(16)+2*e1*z(6)*z(15)*z(16)-4*e3*z(15)*z(29)-2*e1*z(16)*z(27)-2*e2*z(16)*z(24)-2*e4*z(16)*z(26)-z(23)*e4*z(24)-2*e2*z(3)*z(15)*z(16));
z(514) = 8*m*z(16)*(e1*z(29)+e2*z(26)) + 4*z(132)*(e4*z(26)+e1*z(6)*z(16)-2*e3*z(29)-e2*z(3)*z(16));
z(515) = z(130)*z(445) + z(133)*z(446) + 2*z(51)*z(487) - 2*z(513) - z(52)*z(514);
z(516) = (z(444)*(z(44)*z(91)-z(47)*z(87))+z(50)*(z(47)*z(511)-z(44)*z(496)-z(443)*z(91)-2*e1*z(15)*z(87)))/z(50)^2;
z(517) = z(484) + z(137)*z(516) + z(488)*z(92) - z(136)*z(512) - z(515)*z(93);
z(518) = m*(4*e1*z(96)+z(23)*z(9)*z(446)+z(23)*z(10)*z(445)-2*e2*z(15)-2*e2*z(97)-2*z(8)*z(15)*z(445));
z(519) = m*(e4*z(15)*z(24)+2*e3*z(16)*z(27)-2*e3*z(15)*z(27)-e4*z(16)*z(24)-2*e3*z(6)*z(15)*z(16)-e4*z(3)*z(15)*z(16));
z(520) = z(127)*z(446) + 2*z(52)*z(513) - 8*z(519) - 2*z(51)*z(485) - z(126)*z(445);
z(521) = z(520) + z(51)*z(488) + z(137)*z(446) - z(52)*z(515) - z(136)*z(445);
z(522) = (z(521)*z(141)+z(139)*(9.810000000000001*z(518)-z(517)))/z(139)^2;
a14 = z(522);
z(523) = e1*(z(18)+z(9)^2/z(17)^0.5)/z(18)^2;
z(524) = e1*z(9)/(z(17)^0.5*z(18)^2);
z(525) = 4*rad*z(21)*z(523) + 4*z(23)*e2*z(36) + 4*z(211)*z(8)*z(524) + 4*e1*z(21)*z(36) + 4*rad*z(9)*z(36)*z(523) + 4*rad*z(10)*z(15)*z(524) - 4*rad*z(36)*z(524) - 4*rad*z(9)*z(21)*z(524);
z(526) = z(23)*z(525)/(z(38)^0.5*z(39)^2);
z(527) = z(15)*z(525)/(z(38)^0.5*z(39)^2);
z(528) = (4*rad*z(39)*z(523)-z(21)*z(525)/z(38)^0.5)/z(39)^2;
z(529) = 2*e1*z(22) - 2*e1*z(16) - 2*rad*z(7)*z(524);
z(530) = 2*e2*z(22) + 2*rad*z(4)*z(524) - 2*e2*z(16);
z(531) = 2*e3*z(21)*(z(3)*z(43)-z(4)*z(40)) + 2*rad*z(5)*z(523)*(z(3)*z(43)-z(4)*z(40)) + z(529)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(23)*z(3)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(7)*z(527)) + 0.5*z(35)*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527)) - 2*z(23)*e3*(z(5)*z(42)+z(7)*z(40)) - 2*rad*z(2)*z(523)*(z(6)*z(43)-z(7)*z(40)) - z(530)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(2)*z(21)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40)) - 0.5*z(23)*z(6)*(4*e2*z(40)-z(2)*z(526)-z(4)*z(527)) - 0.5*z(5)*z(21)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527)) - 0.5*z(32)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527));
z(532) = 2*e1*z(15) + 2*rad*z(10)*z(523);
z(533) = 2*e3*z(16) - 2*e3*z(22) - 2*rad*z(5)*z(524);
z(534) = 2*e3*z(15)*(z(5)*z(42)+z(7)*z(40)) + 2*rad*z(7)*z(523)*(z(2)*z(42)-z(3)*z(43)) + z(533)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) + 0.5*z(6)*z(15)*(4*e2*z(40)-z(2)*z(526)-z(4)*z(527)) + 0.5*z(7)*z(21)*(4*e3*z(43)-z(2)*z(526)-z(3)*z(528)) + 0.5*z(33)*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527)) - 2*e1*z(21)*(z(2)*z(42)-z(3)*z(43)) - 2*e2*z(21)*(z(5)*z(42)-z(6)*z(43)) - 2*rad*z(4)*z(523)*(z(5)*z(42)-z(6)*z(43)) - 2*rad*z(2)*z(524)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(4)*z(21)*(4*e3*z(42)-z(5)*z(526)-z(6)*z(528)) - 0.5*z(3)*z(15)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(7)*z(527)) - 0.5*z(30)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527));
z(535) = 2*e2*z(21) - 2*z(23)*e1 - 2*rad*z(8)*z(523);
z(536) = 2*e3*z(16) + 2*rad*z(3)*z(524) - 2*e3*z(22);
z(537) = 2*e3*z(15)*(z(3)*z(43)-z(4)*z(40)) + 0.5*z(23)*z(7)*(4*e3*z(43)-z(2)*z(526)-z(3)*z(528)) + 0.5*z(34)*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527)) - 2*z(23)*e1*(z(2)*z(42)-z(3)*z(43)) - 2*z(23)*e2*(z(5)*z(42)-z(6)*z(43)) - z(536)*(z(3)*z(43)-z(2)*z(42)-z(4)*z(40)) - 0.5*z(2)*z(15)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40)) - 2*rad*z(6)*z(524)*(z(6)*z(43)-z(5)*z(42)-z(7)*z(40)) - 0.5*z(23)*z(4)*(4*e3*z(42)-z(5)*z(526)-z(6)*z(528)) - 0.5*z(5)*z(15)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527)) - 0.5*z(31)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527));
z(538) = z(44)*z(535) + z(49)*z(537) - z(47)*z(531) - 2*e2*z(15)*z(45);
z(539) = (z(538)*(z(45)*z(48)-z(46)*z(49))-z(50)*(z(45)*z(532)+z(48)*z(531)-z(46)*z(535)-z(49)*z(534)))/z(50)^2;
z(540) = (z(538)*(z(44)*z(48)-z(46)*z(47))-z(50)*(z(44)*z(532)+z(48)*z(537)-z(47)*z(534)-2*e2*z(15)*z(46)))/z(50)^2;
z(541) = u1*(z(1)*z(539)*u3+2*z(1)*z(540)*u2-2*z(539)*z(102)-z(540)*z(101));
z(542) = u1*(z(1)*z(540)*u2+2*z(1)*z(539)*u3-2*z(540)*z(101)-z(539)*z(102));
z(543) = z(539)*u1*(z(15)*u2+z(109));
z(544) = (1+e2*z(540)+e3*z(539))*u1;
z(545) = 0.5*u2 - 0.5*e1*z(540)*u1 - 0.5*e4*z(539)*u1;
z(546) = -4*e2*z(545) - 2*e1*z(544);
z(547) = (z(12)+z(10)^2/z(11)^0.5)*z(546)/z(12)^2;
z(548) = u1*(z(15)*z(539)*u3+rad*z(547)-z(540)*z(109));
z(549) = rad*z(10)*z(546)/(z(11)^0.5*z(12)^2);
z(550) = (e2*z(539)-e3*z(540))*u1;
z(551) = 0.5*u3 + 0.5*e4*z(540)*u1 - 0.5*e1*z(539)*u1;
z(552) = e2*z(544) + e3*z(550) + 2*e1*z(545) + 2*e4*z(551) + 2*e3p;
z(553) = z(5)*z(549) + z(16)*z(552) + 2*e3*z(61);
z(554) = -4*e3*z(551) - 2*e1*z(544);
z(555) = z(6)*z(549) + z(16)*z(554);
z(556) = 2*e2*z(551) + 2*e3*z(545) - e1*z(550) - e4*z(544) - 2*e1p;
z(557) = z(7)*z(549) + z(16)*z(556) - 2*e1*z(61);
z(558) = u2*z(555) + u1*z(553) + u3*z(557) + z(540)*u1*z(77) - z(539)*u1*z(72);
z(559) = rad*(u2*z(547)-z(539)*u1*z(57));
z(560) = -z(540)*u1*z(110) - z(559);
z(561) = -4*e2*z(545) - 4*e3*z(551);
z(562) = z(2)*z(549) + z(16)*z(561);
z(563) = e2*z(544) + 2*e1*z(545) - 2*e4*z(551) - e3*z(550) - 2*e3p;
z(564) = z(3)*z(549) + z(16)*z(563) - 2*e3*z(61);
z(565) = e2*z(550) + e3*z(544) + 2*e1*z(551) + 2*e4*z(545) + 2*e2p;
z(566) = z(4)*z(549) + z(16)*z(565) + 2*e2*z(61);
z(567) = z(539)*u1*z(74) - u1*z(562) - u2*z(564) - u3*z(566) - z(540)*u1*z(79);
z(568) = u1*(z(23)*z(540)*u1+z(15)*z(539)*u2+z(539)*z(109));
z(569) = u1*(z(23)*z(539)*u1+z(15)*z(539)*u3+rad*z(547)-z(540)*z(109));
z(570) = z(23)*z(540)*u1*u3 - z(23)*z(539)*u1*u2 - z(539)*u1*z(121) - z(540)*u1*z(120) - z(559);
z(571) = 0.25*z(541) + 0.25*z(542) - m*(z(4)*z(24)*z(543)+z(6)*z(27)*z(548)+z(5)*z(27)*z(560)+2*e3*z(24)*z(119)+2*e3*z(27)*z(117)-z(15)*z(548)-2*e1*z(27)*z(118)-2*e2*z(24)*z(118)-z(3)*z(24)*z(548)-z(7)*z(27)*z(543)-2*e3*z(4)*z(16)*z(118)-2*z(24)*z(558)-4*e3*z(16)*z(115)-2*z(6)*z(15)*z(567)-z(2)*z(24)*z(560)-2*e3*z(2)*z(16)*z(117)-2*e3*z(3)*z(16)*z(119)) - m*(z(4)*z(24)*z(568)+z(6)*z(27)*z(569)+2*z(27)*z(567)+z(5)*z(27)*z(570)+2*e3*z(24)*z(124)+2*e3*z(27)*z(122)+4*e3*z(15)*z(115)-z(15)*z(569)-2*e1*z(27)*z(123)-2*e2*z(24)*z(123)-z(3)*z(24)*z(569)-z(7)*z(27)*z(568)-2*e3*z(4)*z(16)*z(123)-2*z(3)*z(15)*z(558)-z(2)*z(24)*z(570)-2*e3*z(2)*z(16)*z(122)-2*e3*z(3)*z(16)*z(124));
z(572) = z(540)*u1*(2*z(1)*u1-z(100));
z(573) = z(540)*u1*(z(1)*u1-z(100));
z(574) = -0.25*z(572) - 0.25*z(573) - m*(z(4)*z(25)*z(543)+z(6)*z(28)*z(548)+z(15)*z(560)+2*z(28)*z(567)+z(5)*z(28)*z(560)+2*e3*z(25)*z(119)+2*e3*z(28)*z(117)+2*z(2)*z(15)*z(558)-2*e1*z(28)*z(118)-2*e2*z(25)*z(118)-z(3)*z(25)*z(548)-z(7)*z(28)*z(543)-2*e3*z(7)*z(16)*z(118)-2*z(25)*z(558)-4*e3*z(16)*z(116)-z(2)*z(25)*z(560)-2*e3*z(5)*z(16)*z(117)-2*e3*z(6)*z(16)*z(119)) - m*(z(23)*z(568)+z(4)*z(25)*z(568)+z(6)*z(28)*z(569)+z(15)*z(570)+z(5)*z(28)*z(570)+2*z(23)*e1*z(116)+2*e3*z(25)*z(124)+2*e3*z(28)*z(122)+2*z(5)*z(15)*z(567)+4*e3*z(15)*z(116)-2*e1*z(28)*z(123)-2*e2*z(25)*z(123)-z(3)*z(25)*z(569)-z(7)*z(28)*z(568)-2*e3*z(7)*z(16)*z(123)-2*z(23)*e2*z(115)-z(23)*z(4)*z(558)-z(23)*z(7)*z(567)-z(2)*z(25)*z(570)-2*e3*z(5)*z(16)*z(122)-2*e3*z(6)*z(16)*z(124));
z(575) = z(539)*u1*(z(1)*u1-z(100));
z(576) = z(539)*u1*(2*z(1)*u1-z(100));
z(577) = m*(z(3)*z(26)*z(548)+z(7)*z(29)*z(543)+2*e1*z(29)*z(118)+2*e2*z(26)*z(118)+2*z(26)*z(558)+z(2)*z(26)*z(560)-z(4)*z(26)*z(543)-z(6)*z(29)*z(548)-2*e1*z(4)*z(16)*z(118)-2*e2*z(7)*z(16)*z(118)-2*z(29)*z(567)-4*e1*z(16)*z(115)-4*e2*z(16)*z(116)-2*e3*z(26)*z(119)-2*e3*z(29)*z(117)-z(5)*z(29)*z(560)-2*e1*z(2)*z(16)*z(117)-2*e1*z(3)*z(16)*z(119)-2*e2*z(5)*z(16)*z(117)-2*e2*z(6)*z(16)*z(119)) + m*(z(3)*z(26)*z(569)+z(7)*z(29)*z(568)+2*e1*z(29)*z(123)+2*e2*z(26)*z(123)+z(2)*z(26)*z(570)+2*z(23)*e3*z(115)-z(23)*z(569)-z(4)*z(26)*z(568)-z(6)*z(29)*z(569)-2*e1*z(4)*z(16)*z(123)-2*e2*z(7)*z(16)*z(123)-2*e3*z(26)*z(124)-2*e3*z(29)*z(122)-z(23)*z(3)*z(558)-z(23)*z(6)*z(567)-z(5)*z(29)*z(570)-2*e1*z(2)*z(16)*z(122)-2*e1*z(3)*z(16)*z(124)-2*e2*z(5)*z(16)*z(122)-2*e2*z(6)*z(16)*z(124)) - 0.25*z(575) - 0.25*z(576);
z(578) = z(571) + z(51)*z(574) + z(540)*z(134) - z(52)*z(577) - z(539)*z(131);
z(579) = m*(2*e3*z(15)*z(25)+2*e3*z(16)*z(27)+2*e3*z(2)*z(15)*z(16)-2*e3*z(15)*z(27)-2*e3*z(16)*z(25)-z(23)*e1*z(27)-z(23)*e2*z(24)-2*e3*z(6)*z(15)*z(16)-z(23)*e3*z(4)*z(16));
z(580) = m*(z(23)*e1*z(28)+z(23)*e2*z(25)+2*e3*z(15)*z(28)+z(23)*e3*z(7)*z(16)-2*e3*z(16)*z(28)-2*e3*z(5)*z(15)*z(16));
z(581) = m*(2*e1*z(16)*z(25)+2*e3*z(16)*z(29)+z(23)*e1*z(4)*z(16)+z(23)*e2*z(7)*z(16)+z(23)*e3*z(6)*z(16)-2*e2*z(16)*z(28)-2*e3*z(15)*z(29)-z(23)*e1*z(29)-z(23)*e2*z(26)-z(23)*e3*z(25)-2*e1*z(2)*z(15)*z(16)-2*e2*z(5)*z(15)*z(16));
z(582) = 2*z(52)*z(581) + 4*z(51)*z(580) - 2*z(579) - z(129)*z(539) - z(130)*z(540);
z(583) = e1*z(550) + e4*z(544) + 2*e2*z(551) + 2*e3*z(545) + 2*e1p;
z(584) = (z(9)*z(17)*z(583)+2*e1*z(17)*z(58)+2*e1*z(9)^2*z(58))/z(17)^1.5;
z(585) = (z(18)*(z(9)*z(584)+z(18)*z(583)+2*e1*z(59))+2*e1*z(9)*(z(18)*z(58)+2*z(9)*z(59))/z(17)^0.5)/z(18)^3;
z(586) = z(21)*z(546) + rad*z(9)*z(547) + z(15)*z(583) + rad*z(10)*z(585) + 2*rad*e1*z(57) + 2*rad*z(523)*z(55);
z(587) = e3*z(544) + 2*e1*z(551) - 2*e4*z(545) - e2*z(550) - 2*e2p;
z(588) = z(23)*z(546) + 2*rad*e2*z(57) - rad*z(8)*z(547) - z(15)*z(587);
z(589) = 2*rad*e2*z(60) - z(23)*z(583) - z(21)*z(587) - 2*rad*z(523)*z(63) - rad*z(8)*z(585);
z(590) = u1*z(586) + u2*z(588) + u3*z(589) + z(540)*u1*z(90) - z(539)*u1*z(88);
z(591) = rad*(z(18)*z(584)+4*e1*z(9)*z(59)/z(17)^0.5)/z(18)^3;
z(592) = z(553) - 2*e3*z(62) - z(5)*z(591) - z(22)*z(552) - 2*rad*z(524)*z(68);
z(593) = 2*rad*z(15)*z(547) + 2*rad*z(21)*z(585) + 2*z(21)*z(36)*z(583) + 4*z(211)*z(524)*z(63) + 4*z(291)*z(523)*z(60) + 2*rad*z(9)*z(36)*z(585) + 4*rad*e1*z(36)*z(60) + 4*rad*z(15)*z(524)*z(55) + 4*rad*z(36)*z(523)*z(58) + 4*z(291)*z(10)*z(524)*z(57) + 2*z(36)*(z(549)-z(591)) + 2*z(9)*z(21)*(z(549)-z(591)) + 4*z(23)*e2*(z(61)-z(62)) + 4*e1*z(21)*(z(61)-z(62)) + 4*rad*z(9)*z(523)*(z(61)-z(62)) - 2*z(15)*z(36)*z(546) - 2*rad*z(10)*z(36)*z(547) - 2*z(23)*z(36)*z(587) - 4*rad*z(21)*z(524)*z(58) - 4*z(291)*z(9)*z(524)*z(60) - 4*rad*z(524)*(z(61)-z(62)) - 2*z(23)*z(8)*(z(549)-z(591)) - 2*z(10)*z(15)*(z(549)-z(591));
z(594) = z(23)*(2*z(593)/(z(38)^0.5*z(39)^2)-z(525)*(2+z(39)/z(38)^0.5)*z(64)/(z(38)*z(39)^3));
z(595) = (2*z(15)*z(525)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(547)+(z(15)*z(525)*z(64)-2*z(15)*z(38)*z(593)-2*rad*z(38)*z(525)*z(57))/z(38)^1.5))/z(39)^3;
z(596) = (2*z(21)*z(525)*z(64)/z(38)+z(39)*(4*rad*z(39)*z(585)+(z(21)*z(525)*z(64)-2*z(21)*z(38)*z(593)-4*rad*z(38)*z(523)*z(64)-2*rad*z(38)*z(525)*z(60))/z(38)^1.5))/z(39)^3;
z(597) = z(22)*z(561) + z(2)*z(591) + 2*rad*z(524)*z(54) - z(562);
z(598) = z(15)*(z(2)*z(42)+z(4)*z(40))*z(554) + rad*z(6)*(z(2)*z(42)+z(4)*z(40))*z(547) + z(21)*(z(2)*z(42)-z(3)*z(43))*z(556) + 2*rad*e3*(z(5)*z(42)+z(7)*z(40))*z(57) + rad*z(7)*(z(2)*z(42)-z(3)*z(43))*z(585) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(592) + 2*rad*z(523)*(z(2)*z(42)-z(3)*z(43))*z(53) + 0.5*z(15)*(4*e2*z(40)-z(2)*z(526)-z(4)*z(527))*z(71) + 0.5*z(21)*(4*e3*z(43)-z(2)*z(526)-z(3)*z(528))*z(53) + 0.5*rad*z(6)*(4*e2*z(40)-z(2)*z(526)-z(4)*z(527))*z(57) + 0.5*rad*z(7)*(4*e3*z(43)-z(2)*z(526)-z(3)*z(528))*z(60) + rad*z(7)*z(523)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.5*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527))*z(85) + rad*z(2)*z(524)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.25*z(6)*z(15)*(4*z(42)*z(561)+z(4)*z(595)+4*e2*z(70)+4*z(40)*z(565)-2*z(526)*z(54)-2*z(527)*z(69)-z(2)*z(594)) + 0.25*z(7)*z(21)*(4*z(42)*z(561)+4*e3*z(67)-4*z(43)*z(563)-2*z(526)*z(54)-2*z(528)*z(66)-z(2)*z(594)-z(3)*z(596)) + 0.25*z(33)*(4*z(42)*z(561)+z(4)*z(595)+4*e2*z(70)+4*e3*z(67)+4*z(40)*z(565)-4*z(43)*z(563)-2*z(526)*z(54)-2*z(527)*z(69)-2*z(528)*z(66)-z(2)*z(594)-z(3)*z(596)) - rad*z(3)*(z(5)*z(42)+z(7)*z(40))*z(547) - z(15)*(z(5)*z(42)+z(7)*z(40))*z(563) - z(21)*(z(5)*z(42)-z(6)*z(43))*z(565) - 2*rad*e1*(z(2)*z(42)-z(3)*z(43))*z(60) - 2*rad*e2*(z(5)*z(42)-z(6)*z(43))*z(60) - 2*rad*z(523)*(z(5)*z(42)-z(6)*z(43))*z(69) - rad*z(4)*(z(5)*z(42)-z(6)*z(43))*z(585) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(597) - 0.5*z(21)*(4*e3*z(42)-z(5)*z(526)-z(6)*z(528))*z(69) - 0.5*rad*z(4)*(4*e3*z(42)-z(5)*z(526)-z(6)*z(528))*z(60) - e3*z(15)*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) - e1*z(21)*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - e2*z(21)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(15)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(7)*z(527))*z(66) - rad*z(4)*z(523)*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*rad*z(3)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(7)*z(527))*z(57) - 0.5*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527))*z(83) - 0.5*z(533)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.25*z(4)*z(21)*(4*z(42)*z(552)-4*z(43)*z(554)-4*e3*z(65)-2*z(526)*z(68)-2*z(528)*z(71)-z(5)*z(594)-z(6)*z(596)) - 0.25*z(3)*z(15)*(z(7)*z(595)+4*z(40)*z(556)+4*z(42)*z(552)-4*e1*z(70)-4*e3*z(65)-2*z(526)*z(68)-2*z(527)*z(53)-z(5)*z(594)) - 0.25*z(30)*(z(7)*z(595)+4*z(40)*z(556)+4*z(42)*z(552)-4*z(43)*z(554)-4*e1*z(70)-4*e3*z(65)-2*z(526)*z(68)-2*z(527)*z(53)-2*z(528)*z(71)-z(5)*z(594)-z(6)*z(596));
z(599) = z(555) - z(22)*z(554) - z(6)*z(591) - 2*rad*z(524)*z(71);
z(600) = z(3)*z(591) + z(22)*z(563) + 2*rad*z(524)*z(66) - z(564) - 2*e3*z(62);
z(601) = rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(547) + z(23)*(z(2)*z(42)-z(3)*z(43))*z(556) + z(15)*(z(3)*z(43)-z(4)*z(40))*z(552) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(599) + 2*rad*e3*(z(3)*z(43)-z(4)*z(40))*z(57) + 0.5*z(23)*(4*e3*z(43)-z(2)*z(526)-z(3)*z(528))*z(53) + 0.5*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527))*z(75) + 0.5*z(536)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + rad*z(6)*z(524)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) + 0.25*z(2)*z(15)*(z(7)*z(595)+4*z(40)*z(556)-4*z(43)*z(554)-4*e1*z(70)-2*z(527)*z(53)-2*z(528)*z(71)-z(6)*z(596)) + 0.25*z(23)*z(7)*(4*z(42)*z(561)+4*e3*z(67)-4*z(43)*z(563)-2*z(526)*z(54)-2*z(528)*z(66)-z(2)*z(594)-z(3)*z(596)) + 0.25*z(34)*(4*z(42)*z(561)+z(4)*z(595)+4*e2*z(70)+4*e3*z(67)+4*z(40)*z(565)-4*z(43)*z(563)-2*z(526)*z(54)-2*z(527)*z(69)-2*z(528)*z(66)-z(2)*z(594)-z(3)*z(596)) - z(15)*(z(6)*z(43)-z(7)*z(40))*z(561) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(547) - z(23)*(z(5)*z(42)-z(6)*z(43))*z(565) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(600) - 0.5*z(15)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40))*z(54) - 0.5*rad*z(2)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40))*z(57) - 0.5*z(23)*(4*e3*z(42)-z(5)*z(526)-z(6)*z(528))*z(69) - e3*z(15)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - z(23)*e1*(2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) - z(23)*e2*(2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.5*z(15)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527))*z(68) - 0.5*rad*z(5)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527))*z(57) - 0.5*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527))*z(73) - 0.25*z(23)*z(4)*(4*z(42)*z(552)-4*z(43)*z(554)-4*e3*z(65)-2*z(526)*z(68)-2*z(528)*z(71)-z(5)*z(594)-z(6)*z(596)) - 0.25*z(5)*z(15)*(z(4)*z(595)+4*e2*z(70)+4*e3*z(67)+4*z(40)*z(565)-4*z(43)*z(563)-2*z(527)*z(69)-2*z(528)*z(66)-z(3)*z(596)) - 0.25*z(31)*(z(7)*z(595)+4*z(40)*z(556)+4*z(42)*z(552)-4*z(43)*z(554)-4*e1*z(70)-4*e3*z(65)-2*z(526)*z(68)-2*z(527)*z(53)-2*z(528)*z(71)-z(5)*z(594)-z(6)*z(596));
z(602) = z(557) + 2*e1*z(62) - z(7)*z(591) - z(22)*z(556) - 2*rad*z(524)*z(53);
z(603) = z(4)*z(591) + z(22)*z(565) + 2*e2*z(62) + 2*rad*z(524)*z(69) - z(566);
z(604) = z(23)*(z(5)*z(42)+z(7)*z(40))*z(563) + z(21)*(z(3)*z(43)-z(4)*z(40))*z(552) + rad*z(5)*(z(3)*z(43)-z(4)*z(40))*z(585) + (z(6)*z(43)-z(5)*z(42)-z(7)*z(40))*z(602) + 2*rad*e3*(z(3)*z(43)-z(4)*z(40))*z(60) + 2*rad*z(523)*(z(3)*z(43)-z(4)*z(40))*z(68) + z(23)*e3*(z(5)*z(65)-2*z(40)*z(53)-2*z(42)*z(68)-z(7)*z(70)) + rad*z(2)*z(523)*(z(7)*z(70)+2*z(40)*z(53)-2*z(43)*z(71)-z(6)*z(67)) + 0.5*z(23)*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(7)*z(527))*z(66) + 0.5*(4*e2*z(40)+4*e3*z(43)-z(2)*z(526)-z(3)*z(528)-z(4)*z(527))*z(80) + 0.5*z(530)*(z(4)*z(70)+2*z(40)*z(69)+2*z(42)*z(54)-2*z(43)*z(66)-z(2)*z(65)-z(3)*z(67)) + 0.25*z(2)*z(21)*(z(7)*z(595)+4*z(40)*z(556)-4*z(43)*z(554)-4*e1*z(70)-2*z(527)*z(53)-2*z(528)*z(71)-z(6)*z(596)) + 0.25*z(23)*z(3)*(z(7)*z(595)+4*z(40)*z(556)+4*z(42)*z(552)-4*e1*z(70)-4*e3*z(65)-2*z(526)*z(68)-2*z(527)*z(53)-z(5)*z(594)) + 0.25*z(35)*(4*z(42)*z(561)+z(4)*z(595)+4*e2*z(70)+4*e3*z(67)+4*z(40)*z(565)-4*z(43)*z(563)-2*z(526)*z(54)-2*z(527)*z(69)-2*z(528)*z(66)-z(2)*z(594)-z(3)*z(596)) - z(23)*(z(2)*z(42)+z(4)*z(40))*z(554) - z(21)*(z(6)*z(43)-z(7)*z(40))*z(561) - 2*rad*z(523)*(z(6)*z(43)-z(7)*z(40))*z(54) - rad*z(2)*(z(6)*z(43)-z(7)*z(40))*z(585) - (z(3)*z(43)-z(2)*z(42)-z(4)*z(40))*z(603) - 0.5*z(21)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40))*z(54) - 0.5*rad*z(2)*(z(6)*z(528)+z(7)*z(527)+4*e1*z(40))*z(60) - 0.5*z(23)*(4*e2*z(40)-z(2)*z(526)-z(4)*z(527))*z(71) - e3*z(21)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - rad*z(5)*z(523)*(z(4)*z(70)+2*z(40)*z(69)-2*z(43)*z(66)-z(3)*z(67)) - 0.5*z(21)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527))*z(68) - 0.5*rad*z(5)*(4*e2*z(40)+4*e3*z(43)-z(3)*z(528)-z(4)*z(527))*z(60) - 0.5*(4*e3*z(42)-4*e1*z(40)-z(5)*z(526)-z(6)*z(528)-z(7)*z(527))*z(78) - 0.5*z(529)*(z(7)*z(70)+2*z(40)*z(53)+2*z(42)*z(68)-2*z(43)*z(71)-z(5)*z(65)-z(6)*z(67)) - 0.25*z(23)*z(6)*(4*z(42)*z(561)+z(4)*z(595)+4*e2*z(70)+4*z(40)*z(565)-2*z(526)*z(54)-2*z(527)*z(69)-z(2)*z(594)) - 0.25*z(5)*z(21)*(z(4)*z(595)+4*e2*z(70)+4*e3*z(67)+4*z(40)*z(565)-4*z(43)*z(563)-2*z(527)*z(69)-2*z(528)*z(66)-z(3)*z(596)) - 0.25*z(32)*(z(7)*z(595)+4*z(40)*z(556)+4*z(42)*z(552)-4*z(43)*z(554)-4*e1*z(70)-4*e3*z(65)-2*z(526)*z(68)-2*z(527)*z(53)-2*z(528)*z(71)-z(5)*z(594)-z(6)*z(596));
z(605) = u1*z(598) + u2*z(601) + u3*z(604) + z(540)*u1*z(81) - z(539)*u1*z(76);
z(606) = (z(538)*(z(45)*z(91)-z(49)*z(87))-z(50)*(z(45)*z(590)+z(531)*z(91)-z(49)*z(605)-z(535)*z(87)))/z(50)^2;
z(607) = m*(z(23)*e3*z(24)+2*e2*z(16)*z(27)+2*e3*z(16)*z(26)-2*e1*z(16)*z(24)-2*e3*z(15)*z(26)-2*e1*z(3)*z(15)*z(16)-2*e2*z(6)*z(15)*z(16)-z(23)*e3*z(3)*z(16));
z(608) = 4*z(132)*(e3*z(26)+e1*z(3)*z(16)+e2*z(6)*z(16)) - 8*m*z(16)*(e1*z(26)-e2*z(29));
z(609) = 2*z(607) + z(130)*z(539) + z(133)*z(540) - 2*z(51)*z(581) - z(52)*z(608);
z(610) = (z(538)*(z(44)*z(91)-z(47)*z(87))-z(50)*(z(44)*z(590)+z(537)*z(91)-z(47)*z(605)-2*e2*z(15)*z(87)))/z(50)^2;
z(611) = z(578) + z(137)*z(610) + z(582)*z(92) - z(136)*z(606) - z(609)*z(93);
z(612) = m*(z(23)*z(9)*z(540)+z(23)*z(10)*z(539)-4*e2*z(96)-2*e1*z(15)-2*e1*z(97)-2*z(8)*z(15)*z(539));
z(613) = m*e3*(z(15)*z(24)-z(16)*z(24)-z(3)*z(15)*z(16));
z(614) = z(127)*z(540) - 8*z(613) - 2*z(51)*z(579) - 2*z(52)*z(607) - z(126)*z(539);
z(615) = z(614) + z(51)*z(582) + z(137)*z(540) - z(52)*z(609) - z(136)*z(539);
z(616) = (z(615)*z(141)+z(139)*(9.810000000000001*z(612)-z(611)))/z(139)^2;
a15 = z(616);
a21 = z(146);
a22 = 0.5*z(240);
a23 = z(350);
a24 = z(450);
a25 = 0.5*z(544);
a31 = z(147);
a32 = z(241);
a33 = -0.5*z(351);
a34 = 0.5*z(451);
a35 = z(545);
a41 = z(154);
a42 = z(248);
a43 = -0.5*z(358);
a44 = -0.5*z(457);
a45 = z(551);
a51 = z(153);
a52 = 0.5*z(247);
a53 = z(357);
a54 = z(456);
a55 = 0.5*z(550);
Output(1)=T; Output(2)=ke; Output(3)=pe; Output(4)=ke+pe; Output(5)=(-(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*z(15)*u2*z(5)-(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*z(15)*u2*z(2)+(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*(z(15)*u1-z(23)*u3)*z(6)+(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*(z(15)*u1-z(23)*u3)*z(3)+(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*z(23)*u2*z(7)+(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*z(23)*u2*z(4)-(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*z(21)*u3*z(5)-(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*z(21)*u3*z(2)+(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*z(21)*u1*z(7)+(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*z(21)*u1*z(4)+(z(6)*z(21)-z(23)*z(5)-z(7)*z(15))*(z(30)*u1+z(31)*u2+z(32)*u3)+(z(23)*z(2)+z(4)*z(15)-z(3)*z(21))*(z(33)*u1+z(34)*u2+z(35)*u3)); Output(6)=(-z(15)*u2*z(8)+(z(15)*u1-z(23)*u3)*z(9)+z(23)*u2*z(10)-z(21)*u3*z(8)+z(21)*u1*z(10)); 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; Output(21)=a11; Output(22)=a12; Output(23)=a13; Output(24)=a14; Output(25)=a15; Output(26)=a21; Output(27)=a22; Output(28)=a23; Output(29)=a24; Output(30)=a25; Output(31)=a31; Output(32)=a32; Output(33)=a33; Output(34)=a34; Output(35)=a35; Output(36)=a41; Output(37)=a42; Output(38)=a43; Output(39)=a44; Output(40)=a45; Output(41)=a51; Output(42)=a52; Output(43)=a53; Output(44)=a54; Output(45)=a55; Output(46)=q1; Output(47)=q2; Output(48)=q3;
FileIdentifier = fopen('all');
% WriteOutput( 1, Output(1:48) );
WriteOutput( FileIdentifier(1), Output(1:48) );
%===========================================================================
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 TwinDisksEig.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( ''TwinDisksEig.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 = 20; % sys(1) is the number of continuous states
sizes.NumDiscStates = 0; % sys(2) is the number of discrete states
sizes.NumOutputs = 48; % 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;