-
Notifications
You must be signed in to change notification settings - Fork 7
/
Copy pathsolver_H4lv1v2ks1_det_equations_from_macaulay2_k_116x125.m
603 lines (503 loc) · 73.9 KB
/
solver_H4lv1v2ks1_det_equations_from_macaulay2_k_116x125.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
% Generated using GBSolver generator Copyright Martin Bujnak,
% Zuzana Kukelova, Tomas Pajdla CTU Prague 2008.
%
% Please refer to the following paper, when using this code :
% Kukelova Z., Bujnak M., Pajdla T., Automatic Generator of Minimal Problem Solvers,
% ECCV 2008, Marseille, France, October 12-18, 2008
function [s, l1, l2, k, v1, v2, v3, u1, u2, u3] = solver_H4lv1v2ks1_det_equations_from_macaulay2_k_116x125(x11, x1p1, x21, x2p1, x31, x3p1, x41, x4p1)
v1 = 10000;
v2 = 10000;
v3 = 10000;
l1 = 10000;
l2 = 10000;
k = 10000;
u1 = 10000;
u2 = 10000;
u3 = 10000;
s = 10000;
% precalculate polynomial equations coefficients
c(1) = x11(1)*x21(1)*x1p1(1)^2*x2p1(2) - x11(1)*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)*x21(2)*x1p1(1)*x2p1(1)^2 + x11(1)*x21(2)*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(1)*x2p1(2) + x11(1)*x21(1)^2*x1p1(2)*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(1)*x2p1(2) - x11(1)*x21(1)*x1p1(2)*x2p1(2)^2 + x11(1)*x21(1)*x1p1(2)^2*x2p1(2) - x11(1)*x21(2)*x1p1(1)*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(1)*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(1)*x2p1(2);
c(2) = x11(1)*x21(2)*x2p1(2) - x11(1)*x21(2)*x1p1(2) + x11(1)*x1p1(2)*x2p1(2) - x21(1)*x1p1(2)*x2p1(2) - x21(2)*x1p1(1)*x2p1(2) + x21(2)*x1p1(2)*x2p1(1);
c(3) = x11(1)*x21(1)*x2p1(1) - x11(1)*x21(1)*x1p1(1) + x11(1)*x1p1(1)*x2p1(1) - x21(1)*x1p1(1)*x2p1(1);
c(4) = x11(1)*x21(1)^2*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(1)^2*x1p1(2)^2*x2p1(1)^2 + x11(1)*x21(2)^2*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(2)^2*x1p1(1)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(2)^2 - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(2)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(2)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(2)^2*x1p1(2)^2*x2p1(1);
c(5) = x11(1)*x21(1)^2*x1p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1) + x11(1)*x21(1)^2*x1p1(2)^2 + x11(1)*x21(2)^2*x1p1(1)^2 - x11(1)^2*x21(1)*x1p1(2)^2 - x11(1)^2*x21(2)^2*x1p1(1) - x11(2)^2*x21(1)*x1p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1) + x11(1)*x21(2)^2*x1p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2 - x11(2)^2*x21(2)^2*x1p1(1) + x11(1)*x21(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x2p1(1)^2 + x11(1)^2*x21(1)^2*x2p1(1) + x11(1)*x21(1)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x2p1(2)^2 + x11(1)^2*x21(2)^2*x2p1(1) - x11(2)^2*x21(1)*x2p1(1)^2 + x11(2)^2*x21(1)^2*x2p1(1) + x11(1)*x21(2)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x2p1(2)^2 + x11(2)^2*x21(2)^2*x2p1(1) + x11(1)*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(1)^2*x2p1(2)^2 + x11(1)*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x1p1(2)^2*x2p1(1) - x21(1)*x1p1(1)^2*x2p1(1)^2 - x21(1)^2*x1p1(1)*x2p1(1)^2 + x21(1)^2*x1p1(1)^2*x2p1(1) - x21(1)*x1p1(1)^2*x2p1(2)^2 - x21(1)*x1p1(2)^2*x2p1(1)^2 - x21(1)^2*x1p1(1)*x2p1(2)^2 + x21(1)^2*x1p1(2)^2*x2p1(1) - x21(2)^2*x1p1(1)*x2p1(1)^2 + x21(2)^2*x1p1(1)^2*x2p1(1) - x21(1)*x1p1(2)^2*x2p1(2)^2 - x21(2)^2*x1p1(1)*x2p1(2)^2 + x21(2)^2*x1p1(2)^2*x2p1(1);
c(6) = 2*x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2) - x11(1)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(2) - x11(1)^2*x21(2)*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(2)*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(2) + x11(1)^2*x21(1)^2*x1p1(2)*x2p1(1) + 2*x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2) + 2*x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2) - x11(1)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(2) - x11(1)^2*x21(2)*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(2)*x1p1(2)^2*x2p1(1) - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(2) + x11(1)^2*x21(2)^2*x1p1(2)*x2p1(1) - x11(2)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)^2*x21(2)*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(2)*x1p1(1)^2*x2p1(1) - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(2) + x11(2)^2*x21(1)^2*x1p1(2)*x2p1(1) + 2*x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2) - x11(2)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)^2*x21(2)*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(2)*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(2) + x11(2)^2*x21(2)^2*x1p1(2)*x2p1(1);
c(7) = x11(1)*x21(2)*x1p1(1)^2*x2p1(2) - x11(1)*x21(2)*x1p1(2)*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(2)*x2p1(2) - x11(1)^2*x21(1)*x1p1(2)*x2p1(2) - x11(1)^2*x21(2)*x1p1(1)*x2p1(2) + x11(1)^2*x21(2)*x1p1(2)*x2p1(1) - x11(1)*x21(2)*x1p1(2)*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(2) - x11(2)^2*x21(1)*x1p1(2)*x2p1(2) - x11(2)^2*x21(2)*x1p1(1)*x2p1(2) + x11(2)^2*x21(2)*x1p1(2)*x2p1(1);
c(8) = 2*x11(1)*x21(1)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(1)*x2p1(1)^2 + 2*x11(1)*x21(1)^2*x1p1(2)^2*x2p1(1) + 2*x11(1)*x21(2)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(1)*x2p1(2)^2 - 2*x11(2)^2*x21(1)*x1p1(1)*x2p1(1)^2 + 2*x11(1)*x21(2)^2*x1p1(2)^2*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(1)*x2p1(2)^2;
c(9) = x11(1)*x21(1)*x1p1(1)^2*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(1)*x2p1(1) - x11(1)^2*x21(1)*x1p1(1)*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(2)^2 + x11(1)*x21(1)*x1p1(2)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(1)*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)*x2p1(1);
c(10) = x11(1)*x21(1)*x2p1(2) - x11(1)*x21(2)*x1p1(1) - x11(1)*x21(1)*x1p1(2) + x11(1)*x21(2)*x2p1(1) + x11(1)*x1p1(1)*x2p1(2) + x11(1)*x1p1(2)*x2p1(1) - 2*x21(1)*x1p1(1)*x2p1(2);
c(11) = 2*x11(1)*x21(1)^2*x2p1(2) - x11(1)^2*x21(2)*x1p1(1) - x11(2)^2*x21(1)*x1p1(2) - x11(2)^2*x21(2)*x1p1(1) - x11(1)^2*x21(1)*x1p1(2) - x11(1)^2*x21(1)*x2p1(2) + x11(1)^2*x21(2)*x2p1(1) + 2*x11(1)*x21(2)^2*x2p1(2) - x11(2)^2*x21(1)*x2p1(2) + x11(2)^2*x21(2)*x2p1(1) + 2*x11(1)*x1p1(1)^2*x2p1(2) - x11(1)^2*x1p1(1)*x2p1(2) + x11(1)^2*x1p1(2)*x2p1(1) + 2*x11(1)*x1p1(2)^2*x2p1(2) - x11(2)^2*x1p1(1)*x2p1(2) + x11(2)^2*x1p1(2)*x2p1(1) - x21(1)*x1p1(2)*x2p1(1)^2 - x21(1)*x1p1(1)^2*x2p1(2) - x21(2)*x1p1(1)*x2p1(1)^2 + x21(2)*x1p1(1)^2*x2p1(1) - x21(1)^2*x1p1(1)*x2p1(2) + x21(1)^2*x1p1(2)*x2p1(1) - x21(1)*x1p1(2)*x2p1(2)^2 - x21(1)*x1p1(2)^2*x2p1(2) - x21(2)*x1p1(1)*x2p1(2)^2 + x21(2)*x1p1(2)^2*x2p1(1) - x21(2)^2*x1p1(1)*x2p1(2) + x21(2)^2*x1p1(2)*x2p1(1);
c(12) = 2*x11(1)*x21(1)^2*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(1) - 2*x11(1)^2*x21(1)*x1p1(1) + 2*x11(1)*x21(2)^2*x2p1(1) + 2*x11(1)*x1p1(1)^2*x2p1(1) + 2*x11(1)*x1p1(2)^2*x2p1(1) - 2*x21(1)*x1p1(1)*x2p1(1)^2 - 2*x21(1)*x1p1(1)*x2p1(2)^2;
c(13) = 2*x11(1)*x2p1(2) - x21(1)*x1p1(2) - x21(2)*x1p1(1) - x21(1)*x2p1(2) + x21(2)*x2p1(1) - x1p1(1)*x2p1(2) + x1p1(2)*x2p1(1);
c(14) = 2*x11(1)*x2p1(1) - 2*x21(1)*x1p1(1);
c(15) = x11(1)*x21(1)^2 - x11(1)^2*x21(1) + x11(1)*x21(2)^2 - x11(2)^2*x21(1) + x11(1)*x1p1(1)^2 - x11(1)^2*x1p1(1) + x11(1)*x1p1(2)^2 - x11(2)^2*x1p1(1) + x11(1)*x2p1(1)^2 - x21(1)*x1p1(1)^2 + x11(1)^2*x2p1(1) - x21(1)^2*x1p1(1) + x11(1)*x2p1(2)^2 - x21(1)*x1p1(2)^2 + x11(2)^2*x2p1(1) - x21(2)^2*x1p1(1) - x21(1)*x2p1(1)^2 + x21(1)^2*x2p1(1) - x21(1)*x2p1(2)^2 + x21(2)^2*x2p1(1) - x1p1(1)*x2p1(1)^2 + x1p1(1)^2*x2p1(1) - x1p1(1)*x2p1(2)^2 + x1p1(2)^2*x2p1(1);
c(16) = x11(1) - x21(1) - x1p1(1) + x2p1(1);
c(17) = x11(1)*x21(1)*x1p1(1)^2*x2p1(2) - x11(1)*x21(1)*x1p1(2)*x2p1(1)^2 + 2*x11(1)*x21(1)^2*x1p1(2)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(1)^2 + x11(2)*x21(1)*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(1)*x1p1(1)*x2p1(2) - x11(1)^2*x21(1)*x1p1(2)*x2p1(1) - x11(1)*x21(1)*x1p1(2)*x2p1(2)^2 + x11(1)*x21(1)*x1p1(2)^2*x2p1(2) + 2*x11(1)*x21(2)^2*x1p1(2)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(2)^2 + x11(2)*x21(1)*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)*x2p1(2) - x11(2)^2*x21(1)*x1p1(2)*x2p1(1);
c(18) = x11(2)*x21(1)*x2p1(2) - x11(2)*x21(1)*x1p1(2) + x11(1)*x1p1(2)*x2p1(2) - x11(2)*x1p1(1)*x2p1(2) + x11(2)*x1p1(2)*x2p1(1) - x21(1)*x1p1(2)*x2p1(2);
c(19) = x11(1)*x21(1)*x2p1(1) - x11(1)*x21(1)*x1p1(1) + x11(1)*x1p1(1)*x2p1(1) - x21(1)*x1p1(1)*x2p1(1);
c(20) = x11(1)*x21(1)^2*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(1)^2*x1p1(2)^2*x2p1(1)^2 + x11(1)*x21(2)^2*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(2)^2*x1p1(1)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(2)^2 - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(2)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(2)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(2)^2*x1p1(2)^2*x2p1(1);
c(21) = x11(1)*x21(1)^2*x1p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1) + x11(1)*x21(1)^2*x1p1(2)^2 + x11(1)*x21(2)^2*x1p1(1)^2 - x11(1)^2*x21(1)*x1p1(2)^2 - x11(1)^2*x21(2)^2*x1p1(1) - x11(2)^2*x21(1)*x1p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1) + x11(1)*x21(2)^2*x1p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2 - x11(2)^2*x21(2)^2*x1p1(1) + x11(1)*x21(1)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x2p1(1)^2 + x11(1)^2*x21(1)^2*x2p1(1) + x11(1)*x21(1)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x2p1(1)^2 - x11(1)^2*x21(1)*x2p1(2)^2 + x11(1)^2*x21(2)^2*x2p1(1) - x11(2)^2*x21(1)*x2p1(1)^2 + x11(2)^2*x21(1)^2*x2p1(1) + x11(1)*x21(2)^2*x2p1(2)^2 - x11(2)^2*x21(1)*x2p1(2)^2 + x11(2)^2*x21(2)^2*x2p1(1) + x11(1)*x1p1(1)^2*x2p1(1)^2 - x11(1)^2*x1p1(1)*x2p1(1)^2 + x11(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(1)^2*x2p1(2)^2 + x11(1)*x1p1(2)^2*x2p1(1)^2 - x11(1)^2*x1p1(1)*x2p1(2)^2 + x11(1)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)^2*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(2)^2*x2p1(2)^2 - x11(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)^2*x1p1(2)^2*x2p1(1) - x21(1)*x1p1(1)^2*x2p1(1)^2 - x21(1)^2*x1p1(1)*x2p1(1)^2 + x21(1)^2*x1p1(1)^2*x2p1(1) - x21(1)*x1p1(1)^2*x2p1(2)^2 - x21(1)*x1p1(2)^2*x2p1(1)^2 - x21(1)^2*x1p1(1)*x2p1(2)^2 + x21(1)^2*x1p1(2)^2*x2p1(1) - x21(2)^2*x1p1(1)*x2p1(1)^2 + x21(2)^2*x1p1(1)^2*x2p1(1) - x21(1)*x1p1(2)^2*x2p1(2)^2 - x21(2)^2*x1p1(1)*x2p1(2)^2 + x21(2)^2*x1p1(2)^2*x2p1(1);
c(22) = x11(1)*x21(1)^2*x1p1(2)*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(1)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)^2*x21(1)^2*x1p1(1)*x2p1(2) + x11(1)^2*x21(1)^2*x1p1(2)*x2p1(1) + x11(1)*x21(1)^2*x1p1(2)*x2p1(2)^2 + x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1)^2 + x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(2)*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(2)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(1)^2*x21(2)^2*x1p1(1)*x2p1(2) + x11(1)^2*x21(2)^2*x1p1(2)*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(2)^2*x21(1)^2*x1p1(1)*x2p1(2) + x11(2)^2*x21(1)^2*x1p1(2)*x2p1(1) + x11(1)*x21(2)^2*x1p1(2)*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(2)^2*x1p1(2)^2*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(2)^2*x21(2)^2*x1p1(1)*x2p1(2) + x11(2)^2*x21(2)^2*x1p1(2)*x2p1(1);
c(23) = x11(1)*x21(1)^2*x1p1(2)*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(1)^2 + x11(2)*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(2) + x11(2)*x21(1)^2*x1p1(2)*x2p1(1) - x11(1)^2*x21(1)*x1p1(2)*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(2)^2 + x11(2)*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)^2*x1p1(1)*x2p1(2) + x11(2)*x21(2)^2*x1p1(2)*x2p1(1) - x11(2)^2*x21(1)*x1p1(2)*x2p1(2);
c(24) = 2*x11(1)*x21(1)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(1)*x2p1(1)^2 + 2*x11(1)*x21(1)^2*x1p1(2)^2*x2p1(1) + 2*x11(1)*x21(2)^2*x1p1(1)^2*x2p1(1) - 2*x11(1)^2*x21(1)*x1p1(1)*x2p1(2)^2 - 2*x11(2)^2*x21(1)*x1p1(1)*x2p1(1)^2 + 2*x11(1)*x21(2)^2*x1p1(2)^2*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(1)*x2p1(2)^2;
c(25) = x11(1)*x21(1)*x1p1(1)^2*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(1)*x2p1(1) - x11(1)^2*x21(1)*x1p1(1)*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(2)^2 + x11(1)*x21(1)*x1p1(2)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(1)*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)*x2p1(1);
c(26) = x11(1)*x21(1)*x2p1(2) - x11(2)*x21(1)*x1p1(1) - x11(1)*x21(1)*x1p1(2) + x11(2)*x21(1)*x2p1(1) + 2*x11(1)*x1p1(2)*x2p1(1) - x21(1)*x1p1(1)*x2p1(2) - x21(1)*x1p1(2)*x2p1(1);
c(27) = x11(1)*x21(1)^2*x1p1(2) - x11(2)*x21(1)^2*x1p1(1) - 2*x11(1)^2*x21(1)*x1p1(2) + x11(1)*x21(2)^2*x1p1(2) - x11(2)*x21(2)^2*x1p1(1) - 2*x11(2)^2*x21(1)*x1p1(2) + x11(1)*x21(1)^2*x2p1(2) + x11(2)*x21(1)^2*x2p1(1) + x11(1)*x21(2)^2*x2p1(2) + x11(2)*x21(2)^2*x2p1(1) + x11(1)*x1p1(2)*x2p1(1)^2 + x11(1)*x1p1(1)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(1)^2 + x11(2)*x1p1(1)^2*x2p1(1) - x11(1)^2*x1p1(1)*x2p1(2) + x11(1)^2*x1p1(2)*x2p1(1) + x11(1)*x1p1(2)*x2p1(2)^2 + x11(1)*x1p1(2)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(2)^2 + x11(2)*x1p1(2)^2*x2p1(1) - x11(2)^2*x1p1(1)*x2p1(2) + x11(2)^2*x1p1(2)*x2p1(1) - 2*x21(1)*x1p1(2)*x2p1(1)^2 - x21(1)^2*x1p1(1)*x2p1(2) + x21(1)^2*x1p1(2)*x2p1(1) - 2*x21(1)*x1p1(2)*x2p1(2)^2 - x21(2)^2*x1p1(1)*x2p1(2) + x21(2)^2*x1p1(2)*x2p1(1);
c(28) = 2*x11(1)*x21(1)^2*x2p1(1) - 2*x11(2)^2*x21(1)*x1p1(1) - 2*x11(1)^2*x21(1)*x1p1(1) + 2*x11(1)*x21(2)^2*x2p1(1) + 2*x11(1)*x1p1(1)^2*x2p1(1) + 2*x11(1)*x1p1(2)^2*x2p1(1) - 2*x21(1)*x1p1(1)*x2p1(1)^2 - 2*x21(1)*x1p1(1)*x2p1(2)^2;
c(29) = x11(1)*x1p1(2) - x11(2)*x1p1(1) + x11(1)*x2p1(2) + x11(2)*x2p1(1) - 2*x21(1)*x1p1(2) - x1p1(1)*x2p1(2) + x1p1(2)*x2p1(1);
c(30) = 2*x11(1)*x2p1(1) - 2*x21(1)*x1p1(1);
c(31) = x11(1)*x21(1)^2 - x11(1)^2*x21(1) + x11(1)*x21(2)^2 - x11(2)^2*x21(1) + x11(1)*x1p1(1)^2 - x11(1)^2*x1p1(1) + x11(1)*x1p1(2)^2 - x11(2)^2*x1p1(1) + x11(1)*x2p1(1)^2 - x21(1)*x1p1(1)^2 + x11(1)^2*x2p1(1) - x21(1)^2*x1p1(1) + x11(1)*x2p1(2)^2 - x21(1)*x1p1(2)^2 + x11(2)^2*x2p1(1) - x21(2)^2*x1p1(1) - x21(1)*x2p1(1)^2 + x21(1)^2*x2p1(1) - x21(1)*x2p1(2)^2 + x21(2)^2*x2p1(1) - x1p1(1)*x2p1(1)^2 + x1p1(1)^2*x2p1(1) - x1p1(1)*x2p1(2)^2 + x1p1(2)^2*x2p1(1);
c(32) = x11(1) - x21(1) - x1p1(1) + x2p1(1);
c(33) = x11(1)*x21(2)*x1p1(2)*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(2) + x11(2)*x21(2)*x1p1(2)*x2p1(1);
c(34) = x11(1)*x21(1)*x1p1(2)*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(2) + x11(1)*x21(2)*x1p1(1)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(1);
c(35) = x11(1)*x21(2)*x1p1(1)^2*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(2)*x2p1(1)^2 - x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)*x1p1(1)^2*x2p1(1)^2 - x11(2)*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(1)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(1)^2 + x11(1)*x21(1)^2*x1p1(2)*x2p1(2)^2 - x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1)^2 - x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2*x2p1(1)^2 - x11(2)*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(2)*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(2)^2*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)*x1p1(2)^2*x2p1(2)^2 + x11(1)*x21(2)^2*x1p1(2)*x2p1(2)^2 - x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2) - x11(2)*x21(1)*x1p1(2)^2*x2p1(2)^2 - x11(2)*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(2)^2*x1p1(2)^2*x2p1(1);
c(36) = 2*x11(1)*x21(2)*x1p1(2)*x2p1(1) - 2*x11(2)*x21(1)*x1p1(1)*x2p1(2);
c(37) = 2*x11(1)*x21(2)*x1p1(2)*x2p1(1)^2 - x11(2)*x21(1)*x1p1(2)*x2p1(1)^2 - x11(2)*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(1)^2 + x11(2)*x21(2)*x1p1(1)^2*x2p1(1) - x11(2)*x21(1)^2*x1p1(1)*x2p1(2) + x11(2)*x21(1)^2*x1p1(2)*x2p1(1) + 2*x11(1)*x21(2)*x1p1(2)*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(2)^2 + x11(2)*x21(2)*x1p1(2)^2*x2p1(1) - x11(2)*x21(2)^2*x1p1(1)*x2p1(2) + x11(2)*x21(2)^2*x1p1(2)*x2p1(1);
c(38) = x11(1)*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)*x21(1)*x1p1(1)^2*x2p1(2) + x11(1)*x21(2)*x1p1(1)*x2p1(1)^2 + x11(1)*x21(2)*x1p1(1)^2*x2p1(1) - x11(1)*x21(1)^2*x1p1(1)*x2p1(2) + x11(1)*x21(1)^2*x1p1(2)*x2p1(1) - 2*x11(2)*x21(1)*x1p1(1)*x2p1(1)^2 + x11(1)*x21(1)*x1p1(2)*x2p1(2)^2 - x11(1)*x21(1)*x1p1(2)^2*x2p1(2) + x11(1)*x21(2)*x1p1(1)*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(1) - x11(1)*x21(2)^2*x1p1(1)*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1) - 2*x11(2)*x21(1)*x1p1(1)*x2p1(2)^2;
c(39) = 2*x11(1)*x21(2)*x1p1(2) - x11(2)*x21(1)*x1p1(2) - x11(2)*x21(2)*x1p1(1) - x11(2)*x21(1)*x2p1(2) + x11(2)*x21(2)*x2p1(1) - x11(2)*x1p1(1)*x2p1(2) + x11(2)*x1p1(2)*x2p1(1);
c(40) = x11(1)*x21(1)*x1p1(2) + x11(1)*x21(2)*x1p1(1) - 2*x11(2)*x21(1)*x1p1(1) - x11(1)*x21(1)*x2p1(2) + x11(1)*x21(2)*x2p1(1) - x11(1)*x1p1(1)*x2p1(2) + x11(1)*x1p1(2)*x2p1(1);
c(41) = x11(1)*x21(2)*x1p1(1)^2 + x11(1)*x21(1)^2*x1p1(2) - x11(2)*x21(1)*x1p1(1)^2 - x11(2)*x21(1)^2*x1p1(1) + x11(1)*x21(2)*x1p1(2)^2 + x11(1)*x21(2)^2*x1p1(2) - x11(2)*x21(1)*x1p1(2)^2 - x11(2)*x21(2)^2*x1p1(1) + x11(1)*x21(2)*x2p1(1)^2 - x11(1)*x21(1)^2*x2p1(2) - x11(2)*x21(1)*x2p1(1)^2 + x11(2)*x21(1)^2*x2p1(1) + x11(1)*x21(2)*x2p1(2)^2 - x11(1)*x21(2)^2*x2p1(2) - x11(2)*x21(1)*x2p1(2)^2 + x11(2)*x21(2)^2*x2p1(1) + x11(1)*x1p1(2)*x2p1(1)^2 - x11(1)*x1p1(1)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(1)^2 + x11(2)*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(2)*x2p1(2)^2 - x11(1)*x1p1(2)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(2)^2 + x11(2)*x1p1(2)^2*x2p1(1);
c(42) = x11(1)*x21(2) - x11(2)*x21(1) + x11(1)*x1p1(2) - x11(2)*x1p1(1) - x11(1)*x2p1(2) + x11(2)*x2p1(1);
c(43) = x11(1)*x21(2)*x1p1(2)*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(2) + x11(2)*x21(2)*x1p1(2)*x2p1(1);
c(44) = x11(1)*x21(1)*x1p1(2)*x2p1(1) - x11(1)*x21(1)*x1p1(1)*x2p1(2) + x11(1)*x21(2)*x1p1(1)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(1);
c(45) = x11(1)*x21(2)*x1p1(1)^2*x2p1(1)^2 - x11(2)*x21(1)*x1p1(1)^2*x2p1(1)^2 + x11(1)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)^2*x21(1)*x1p1(1)^2*x2p1(2) - x11(1)^2*x21(2)*x1p1(1)*x2p1(1)^2 + x11(1)^2*x21(2)*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)*x1p1(1)^2*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(1)^2 - x11(2)*x21(1)*x1p1(1)^2*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2*x2p1(1)^2 + x11(1)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(1)^2*x21(1)*x1p1(2)^2*x2p1(2) - x11(1)^2*x21(2)*x1p1(1)*x2p1(2)^2 + x11(1)^2*x21(2)*x1p1(2)^2*x2p1(1) + x11(2)^2*x21(1)*x1p1(2)*x2p1(1)^2 - x11(2)^2*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)^2*x21(2)*x1p1(1)*x2p1(1)^2 + x11(2)^2*x21(2)*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)*x1p1(2)^2*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2*x2p1(2)^2 + x11(2)^2*x21(1)*x1p1(2)*x2p1(2)^2 - x11(2)^2*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)^2*x21(2)*x1p1(1)*x2p1(2)^2 + x11(2)^2*x21(2)*x1p1(2)^2*x2p1(1);
c(46) = 2*x11(1)*x21(2)*x1p1(2)*x2p1(1) - 2*x11(2)*x21(1)*x1p1(1)*x2p1(2);
c(47) = x11(1)*x21(2)*x1p1(2)*x2p1(1)^2 + x11(1)*x21(2)*x1p1(1)^2*x2p1(2) - 2*x11(2)*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(1)^2 + x11(2)*x21(2)*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(2)*x1p1(1)*x2p1(2) + x11(1)^2*x21(2)*x1p1(2)*x2p1(1) + x11(1)*x21(2)*x1p1(2)*x2p1(2)^2 + x11(1)*x21(2)*x1p1(2)^2*x2p1(2) - 2*x11(2)*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)*x1p1(1)*x2p1(2)^2 + x11(2)*x21(2)*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(2)*x1p1(1)*x2p1(2) + x11(2)^2*x21(2)*x1p1(2)*x2p1(1);
c(48) = x11(1)*x21(1)*x1p1(2)*x2p1(1)^2 - x11(1)*x21(1)*x1p1(1)^2*x2p1(2) + 2*x11(1)*x21(2)*x1p1(1)^2*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(1)^2 - x11(2)*x21(1)*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(1)*x1p1(1)*x2p1(2) + x11(1)^2*x21(1)*x1p1(2)*x2p1(1) + x11(1)*x21(1)*x1p1(2)*x2p1(2)^2 - x11(1)*x21(1)*x1p1(2)^2*x2p1(2) + 2*x11(1)*x21(2)*x1p1(2)^2*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(1)*x2p1(2) + x11(2)^2*x21(1)*x1p1(2)*x2p1(1);
c(49) = x11(1)*x21(2)*x1p1(2) - x11(2)*x21(2)*x1p1(1) + x11(1)*x21(2)*x2p1(2) - 2*x11(2)*x21(1)*x2p1(2) + x11(2)*x21(2)*x2p1(1) - x21(2)*x1p1(1)*x2p1(2) + x21(2)*x1p1(2)*x2p1(1);
c(50) = x11(1)*x21(1)*x1p1(2) - x11(2)*x21(1)*x1p1(1) - x11(1)*x21(1)*x2p1(2) + 2*x11(1)*x21(2)*x2p1(1) - x11(2)*x21(1)*x2p1(1) - x21(1)*x1p1(1)*x2p1(2) + x21(1)*x1p1(2)*x2p1(1);
c(51) = x11(1)*x21(2)*x1p1(1)^2 - x11(2)*x21(1)*x1p1(1)^2 + x11(1)^2*x21(1)*x1p1(2) - x11(1)^2*x21(2)*x1p1(1) + x11(1)*x21(2)*x1p1(2)^2 - x11(2)*x21(1)*x1p1(2)^2 + x11(2)^2*x21(1)*x1p1(2) - x11(2)^2*x21(2)*x1p1(1) + x11(1)*x21(2)*x2p1(1)^2 - x11(2)*x21(1)*x2p1(1)^2 - x11(1)^2*x21(1)*x2p1(2) + x11(1)^2*x21(2)*x2p1(1) + x11(1)*x21(2)*x2p1(2)^2 - x11(2)*x21(1)*x2p1(2)^2 - x11(2)^2*x21(1)*x2p1(2) + x11(2)^2*x21(2)*x2p1(1) + x21(1)*x1p1(2)*x2p1(1)^2 - x21(1)*x1p1(1)^2*x2p1(2) - x21(2)*x1p1(1)*x2p1(1)^2 + x21(2)*x1p1(1)^2*x2p1(1) + x21(1)*x1p1(2)*x2p1(2)^2 - x21(1)*x1p1(2)^2*x2p1(2) - x21(2)*x1p1(1)*x2p1(2)^2 + x21(2)*x1p1(2)^2*x2p1(1);
c(52) = x11(1)*x21(2) - x11(2)*x21(1) + x21(1)*x1p1(2) - x21(2)*x1p1(1) - x21(1)*x2p1(2) + x21(2)*x2p1(1);
c(53) = x11(1)*x21(1)^2*x1p1(2)*x2p1(1)^2 - x11(1)*x21(1)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(1)^2*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(1)*x1p1(2)*x2p1(1)^2 + x11(1)^2*x21(1)*x1p1(1)^2*x2p1(2) + x11(1)^2*x21(2)*x1p1(1)*x2p1(1)^2 - x11(1)^2*x21(2)*x1p1(1)^2*x2p1(1) + x11(1)*x21(1)^2*x1p1(2)*x2p1(2)^2 - x11(1)*x21(1)^2*x1p1(2)^2*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1)^2 - x11(1)*x21(2)^2*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(1)^2*x1p1(2)^2*x2p1(1) - x11(2)*x21(2)^2*x1p1(1)*x2p1(1)^2 + x11(2)*x21(2)^2*x1p1(1)^2*x2p1(1) - x11(1)^2*x21(1)*x1p1(2)*x2p1(2)^2 + x11(1)^2*x21(1)*x1p1(2)^2*x2p1(2) + x11(1)^2*x21(2)*x1p1(1)*x2p1(2)^2 - x11(1)^2*x21(2)*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(2)*x2p1(1)^2 + x11(2)^2*x21(1)*x1p1(1)^2*x2p1(2) + x11(2)^2*x21(2)*x1p1(1)*x2p1(1)^2 - x11(2)^2*x21(2)*x1p1(1)^2*x2p1(1) + x11(1)*x21(2)^2*x1p1(2)*x2p1(2)^2 - x11(1)*x21(2)^2*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)^2*x1p1(1)*x2p1(2)^2 + x11(2)*x21(2)^2*x1p1(2)^2*x2p1(1) - x11(2)^2*x21(1)*x1p1(2)*x2p1(2)^2 + x11(2)^2*x21(1)*x1p1(2)^2*x2p1(2) + x11(2)^2*x21(2)*x1p1(1)*x2p1(2)^2 - x11(2)^2*x21(2)*x1p1(2)^2*x2p1(1);
c(54) = x11(1)*x21(2)*x1p1(2)*x2p1(1)^2 - x11(1)*x21(2)*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(1)^2 + x11(2)*x21(1)*x1p1(1)^2*x2p1(2) - x11(2)*x21(1)^2*x1p1(1)*x2p1(2) + x11(2)*x21(1)^2*x1p1(2)*x2p1(1) + x11(1)^2*x21(2)*x1p1(1)*x2p1(2) - x11(1)^2*x21(2)*x1p1(2)*x2p1(1) + x11(1)*x21(2)*x1p1(2)*x2p1(2)^2 - x11(1)*x21(2)*x1p1(2)^2*x2p1(2) - x11(2)*x21(1)*x1p1(2)*x2p1(2)^2 + x11(2)*x21(1)*x1p1(2)^2*x2p1(2) - x11(2)*x21(2)^2*x1p1(1)*x2p1(2) + x11(2)*x21(2)^2*x1p1(2)*x2p1(1) + x11(2)^2*x21(2)*x1p1(1)*x2p1(2) - x11(2)^2*x21(2)*x1p1(2)*x2p1(1);
c(55) = x11(1)*x21(2)*x1p1(1)*x2p1(1)^2 - x11(1)*x21(2)*x1p1(1)^2*x2p1(1) - x11(1)*x21(1)^2*x1p1(1)*x2p1(2) + x11(1)*x21(1)^2*x1p1(2)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(1)^2 + x11(2)*x21(1)*x1p1(1)^2*x2p1(1) + x11(1)^2*x21(1)*x1p1(1)*x2p1(2) - x11(1)^2*x21(1)*x1p1(2)*x2p1(1) + x11(1)*x21(2)*x1p1(1)*x2p1(2)^2 - x11(1)*x21(2)*x1p1(2)^2*x2p1(1) - x11(1)*x21(2)^2*x1p1(1)*x2p1(2) + x11(1)*x21(2)^2*x1p1(2)*x2p1(1) - x11(2)*x21(1)*x1p1(1)*x2p1(2)^2 + x11(2)*x21(1)*x1p1(2)^2*x2p1(1) + x11(2)^2*x21(1)*x1p1(1)*x2p1(2) - x11(2)^2*x21(1)*x1p1(2)*x2p1(1);
c(56) = x11(1)*x21(2)*x1p1(2) - x11(2)*x21(1)*x1p1(2) - x11(1)*x21(2)*x2p1(2) + x11(2)*x21(1)*x2p1(2) - x11(2)*x1p1(1)*x2p1(2) + x11(2)*x1p1(2)*x2p1(1) + x21(2)*x1p1(1)*x2p1(2) - x21(2)*x1p1(2)*x2p1(1);
c(57) = x11(1)*x21(2)*x1p1(1) - x11(2)*x21(1)*x1p1(1) - x11(1)*x21(2)*x2p1(1) + x11(2)*x21(1)*x2p1(1) - x11(1)*x1p1(1)*x2p1(2) + x11(1)*x1p1(2)*x2p1(1) + x21(1)*x1p1(1)*x2p1(2) - x21(1)*x1p1(2)*x2p1(1);
c(58) = x11(1)*x21(1)^2*x1p1(2) - x11(2)*x21(1)^2*x1p1(1) - x11(1)^2*x21(1)*x1p1(2) + x11(1)^2*x21(2)*x1p1(1) + x11(1)*x21(2)^2*x1p1(2) - x11(2)*x21(2)^2*x1p1(1) - x11(2)^2*x21(1)*x1p1(2) + x11(2)^2*x21(2)*x1p1(1) - x11(1)*x21(1)^2*x2p1(2) + x11(2)*x21(1)^2*x2p1(1) + x11(1)^2*x21(1)*x2p1(2) - x11(1)^2*x21(2)*x2p1(1) - x11(1)*x21(2)^2*x2p1(2) + x11(2)*x21(2)^2*x2p1(1) + x11(2)^2*x21(1)*x2p1(2) - x11(2)^2*x21(2)*x2p1(1) + x11(1)*x1p1(2)*x2p1(1)^2 - x11(1)*x1p1(1)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(1)^2 + x11(2)*x1p1(1)^2*x2p1(1) + x11(1)*x1p1(2)*x2p1(2)^2 - x11(1)*x1p1(2)^2*x2p1(2) - x11(2)*x1p1(1)*x2p1(2)^2 + x11(2)*x1p1(2)^2*x2p1(1) - x21(1)*x1p1(2)*x2p1(1)^2 + x21(1)*x1p1(1)^2*x2p1(2) + x21(2)*x1p1(1)*x2p1(1)^2 - x21(2)*x1p1(1)^2*x2p1(1) - x21(1)*x1p1(2)*x2p1(2)^2 + x21(1)*x1p1(2)^2*x2p1(2) + x21(2)*x1p1(1)*x2p1(2)^2 - x21(2)*x1p1(2)^2*x2p1(1);
c(59) = x11(1)*x1p1(2) - x11(2)*x1p1(1) - x11(1)*x2p1(2) + x11(2)*x2p1(1) - x21(1)*x1p1(2) + x21(2)*x1p1(1) + x21(1)*x2p1(2) - x21(2)*x2p1(1);
c(60) = x31(1)^2*x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(2) + x31(1)^2*x41(2)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(61) = x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x4p1(1) - x31(1)*x41(1)*x4p1(2) + x41(2)*x3p1(1)*x4p1(1);
c(62) = x31(1)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x3p1(1)^2 - x31(1)*x41(2)*x3p1(2)^2 + x31(2)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(2) - x31(1)*x41(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(2) - x31(1)*x3p1(1)^2*x4p1(2) + x31(1)^2*x3p1(1)*x4p1(2) - x31(1)*x3p1(2)^2*x4p1(2) + x31(2)^2*x3p1(1)*x4p1(2) + x41(2)*x3p1(1)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(2) + x41(2)*x3p1(1)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(2);
c(63) = x31(1)^2*x41(1)*x3p1(1) - x31(1)*x41(1)*x3p1(1)^2 - x31(1)*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1) - x31(1)*x41(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(1) - x31(1)*x41(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(1) + x31(1)^2*x3p1(1)*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(1) + x31(2)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(1);
c(64) = x31(1)*x41(1)*x3p1(2)*x4p1(1)^2 + x31(1)*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(1)^2*x3p1(2)*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)*x4p1(2)^2 + x31(1)*x41(2)*x3p1(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)*x4p1(2) - x31(1)*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(65) = x31(1)^2*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2) + x31(1)^2*x41(2)*x3p1(1)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(2);
c(66) = x31(1)^2*x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(2)^2*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(2);
c(67) = x31(1)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1);
c(68) = x31(1)^2*x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1);
c(69) = x31(1)*x41(2)*x3p1(2) - x31(1)*x3p1(2)*x4p1(2) + x41(1)*x3p1(2)*x4p1(2) - x41(2)*x3p1(2)*x4p1(1);
c(70) = x31(1)*x41(1)*x3p1(1) - x31(1)*x3p1(1)*x4p1(1);
c(71) = x31(1)^2*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)^2*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)^2*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(2)^2*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)^2*x41(2)^2*x3p1(2)^2*x4p1(1);
c(72) = x31(1)^2*x41(1)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1)^2 + x31(2)^2*x41(1)*x3p1(2)^2 + x31(1)^2*x41(1)*x4p1(1)^2 - x31(1)^2*x41(1)^2*x4p1(1) + x31(1)^2*x41(1)*x4p1(2)^2 - x31(1)^2*x41(2)^2*x4p1(1) + x31(2)^2*x41(1)*x4p1(1)^2 - x31(2)^2*x41(1)^2*x4p1(1) + x31(2)^2*x41(1)*x4p1(2)^2 - x31(2)^2*x41(2)^2*x4p1(1) - x31(1)^2*x3p1(1)^2*x4p1(1) - x31(1)^2*x3p1(2)^2*x4p1(1) - x31(2)^2*x3p1(1)^2*x4p1(1) - x31(2)^2*x3p1(2)^2*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(1)^2 - x41(1)^2*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(2)^2 + x41(1)*x3p1(2)^2*x4p1(1)^2 - x41(1)^2*x3p1(2)^2*x4p1(1) - x41(2)^2*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(2)^2*x4p1(2)^2 - x41(2)^2*x3p1(2)^2*x4p1(1);
c(73) = x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x4p1(2);
c(74) = x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x4p1(1);
c(75) = x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(2)^2;
c(76) = x31(1)^2*x41(1)^2*x3p1(1) - x31(1)*x41(1)^2*x3p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)^2 + x31(1)^2*x41(2)^2*x3p1(1) + x31(2)^2*x41(1)^2*x3p1(1) - x31(1)*x41(2)^2*x3p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1) - x31(1)*x41(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1)^2 - x31(1)*x41(2)^2*x4p1(2)^2 - x31(1)*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x3p1(2)^2*x4p1(1)^2 + x31(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x3p1(2)^2*x4p1(2)^2 + x31(2)^2*x3p1(1)*x4p1(2)^2 + x41(1)^2*x3p1(1)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(1)^2 + x41(2)^2*x3p1(1)*x4p1(2)^2;
c(77) = x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2) + x31(1)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(1)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)^2*x4p1(1) - x31(1)^2*x41(1)^2*x3p1(2)*x4p1(1) + x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 + x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)^2*x4p1(1) - x31(1)^2*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)^2*x4p1(1) - x31(2)^2*x41(1)^2*x3p1(2)*x4p1(1) + x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)^2*x4p1(1) - x31(2)^2*x41(2)^2*x3p1(2)*x4p1(1);
c(78) = x31(1)*x41(2)*x3p1(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(2) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)*x4p1(1) + x31(1)*x41(2)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(2) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)*x4p1(1);
c(79) = x31(1)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(1)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2)^2 - x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1);
c(80) = x31(1)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)*x4p1(1);
c(81) = x41(2)*x3p1(1) - x31(1)*x4p1(2) - x31(1)*x41(2) + x3p1(1)*x4p1(2);
c(82) = x41(1)*x3p1(1) - x31(1)*x4p1(1) - x31(1)*x41(1) + x3p1(1)*x4p1(1);
c(83) = x31(1)^2*x3p1(1) - x31(1)*x41(2)^2 - x31(1)*x3p1(1)^2 - x31(1)*x41(1)^2 - x31(1)*x3p1(2)^2 + x31(2)^2*x3p1(1) - x31(1)*x4p1(1)^2 + x41(1)^2*x3p1(1) - x31(1)*x4p1(2)^2 + x41(2)^2*x3p1(1) + x3p1(1)*x4p1(1)^2 + x3p1(1)*x4p1(2)^2;
c(84) = x31(1)*x41(1)*x3p1(2) + x31(1)*x41(2)*x3p1(1) - x31(1)*x3p1(1)*x4p1(2) - x31(1)*x3p1(2)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2) - x41(2)*x3p1(1)*x4p1(1);
c(85) = x31(1)*x41(2)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(2) + x31(1)*x41(2)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(2) + x31(1)*x41(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(2) + x31(1)^2*x41(1)*x4p1(2) - x31(1)^2*x41(2)*x4p1(1) + x31(1)*x41(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(2) + x31(2)^2*x41(1)*x4p1(2) - x31(2)^2*x41(2)*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(2) - x31(1)^2*x3p1(2)*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(2) - x31(2)^2*x3p1(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(1)^2 + x41(1)*x3p1(1)^2*x4p1(2) - x41(2)*x3p1(1)^2*x4p1(1) - x41(1)^2*x3p1(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(2)^2 + x41(1)*x3p1(2)^2*x4p1(2) - x41(2)*x3p1(2)^2*x4p1(1) - x41(2)^2*x3p1(2)*x4p1(1);
c(86) = x31(1)*x41(1)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(1) + x31(1)*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1) + x31(1)*x41(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(1) + x31(1)*x41(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(1) - x31(1)^2*x3p1(1)*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(1) - x31(2)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(1)^2 - x41(1)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2)^2 - x41(2)^2*x3p1(1)*x4p1(1);
c(87) = x3p1(1) - x31(1);
c(88) = x31(1)*x41(2) - x31(1)*x4p1(2) + x41(1)*x3p1(2) + x41(1)*x4p1(2) - x41(2)*x4p1(1) - x3p1(2)*x4p1(1);
c(89) = x31(1)*x41(1) - x31(1)*x4p1(1) + x41(1)*x3p1(1) - x3p1(1)*x4p1(1);
c(90) = x31(1)^2*x41(1) + x31(2)^2*x41(1) + x41(1)*x3p1(1)^2 - x31(1)^2*x4p1(1) + x41(1)*x3p1(2)^2 - x31(2)^2*x4p1(1) + x41(1)*x4p1(1)^2 - x41(1)^2*x4p1(1) + x41(1)*x4p1(2)^2 - x41(2)^2*x4p1(1) - x3p1(1)^2*x4p1(1) - x3p1(2)^2*x4p1(1);
c(91) = x41(1) - x4p1(1);
c(92) = x31(2)*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)^2*x3p1(2)*x4p1(1) - x31(2)*x41(1)*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(2) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2) + x31(1)^2*x41(1)*x3p1(2)*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)^2*x3p1(2)*x4p1(1) - x31(2)*x41(1)*x3p1(2)^2*x4p1(1) + x31(2)*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2) + x31(2)^2*x41(1)*x3p1(2)*x4p1(1);
c(93) = x31(2)*x3p1(1)*x4p1(1) - x31(2)*x41(1)*x4p1(1) - x31(1)*x3p1(2)*x4p1(1) - x31(1)*x41(1)*x4p1(2) + x41(1)*x3p1(1)*x4p1(2) + x41(1)*x3p1(2)*x4p1(1);
c(94) = x31(2)*x41(1)^2*x3p1(1) - x31(2)*x41(1)*x3p1(1)^2 - x31(1)*x41(1)^2*x3p1(2) + x31(1)^2*x41(1)*x3p1(2) - x31(1)*x41(2)^2*x3p1(2) - x31(2)*x41(1)*x3p1(2)^2 + x31(2)*x41(2)^2*x3p1(1) + x31(2)^2*x41(1)*x3p1(2) - x31(1)*x41(1)^2*x4p1(2) - x31(2)*x41(1)*x4p1(1)^2 - x31(1)*x41(2)^2*x4p1(2) - x31(2)*x41(1)*x4p1(2)^2 - x31(1)*x3p1(2)*x4p1(1)^2 - x31(1)*x3p1(1)^2*x4p1(2) + x31(2)*x3p1(1)*x4p1(1)^2 + x31(1)^2*x3p1(1)*x4p1(2) - x31(1)*x3p1(2)*x4p1(2)^2 - x31(1)*x3p1(2)^2*x4p1(2) + x31(2)*x3p1(1)*x4p1(2)^2 + x31(2)^2*x3p1(1)*x4p1(2) + x41(1)*x3p1(2)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(2) + x41(1)*x3p1(2)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(2);
c(95) = x31(1)^2*x41(1)*x3p1(1) - x31(1)*x41(1)*x3p1(1)^2 - x31(1)*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1) - x31(1)*x41(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(1) - x31(1)*x41(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(1) + x31(1)^2*x3p1(1)*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(1) + x31(2)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(1);
c(96) = x31(1)*x41(1)*x3p1(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(1)*x4p1(1);
c(97) = x31(2)*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2) - x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(1)^2 + x31(1)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(1)^2*x3p1(2)*x4p1(2)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)^2*x3p1(2)*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2) - x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 - x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(1)^2 + x31(1)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(2) + x31(2)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(2)^2*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2) - x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(2);
c(98) = x31(2)*x41(1)^2*x3p1(1)*x4p1(2) - x31(2)*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)*x41(1)^2*x3p1(2)*x4p1(2) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2) - x31(1)*x41(2)^2*x3p1(2)*x4p1(2) - x31(2)*x41(1)*x3p1(2)^2*x4p1(2) + x31(2)*x41(2)^2*x3p1(1)*x4p1(2) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2);
c(99) = x31(1)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1);
c(100) = x31(1)^2*x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1);
c(101) = x31(2)*x41(1)*x3p1(2) - x31(2)*x3p1(2)*x4p1(1);
c(102) = x31(1)*x41(1)*x3p1(1) - x31(1)*x3p1(1)*x4p1(1);
c(103) = x31(1)^2*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)^2*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)^2*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(2)^2*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)^2*x41(2)^2*x3p1(2)^2*x4p1(1);
c(104) = x31(1)^2*x41(1)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1)^2 + x31(2)^2*x41(1)*x3p1(2)^2 + x31(1)^2*x41(1)*x4p1(1)^2 - x31(1)^2*x41(1)^2*x4p1(1) + x31(1)^2*x41(1)*x4p1(2)^2 - x31(1)^2*x41(2)^2*x4p1(1) + x31(2)^2*x41(1)*x4p1(1)^2 - x31(2)^2*x41(1)^2*x4p1(1) + x31(2)^2*x41(1)*x4p1(2)^2 - x31(2)^2*x41(2)^2*x4p1(1) - x31(1)^2*x3p1(1)^2*x4p1(1) - x31(1)^2*x3p1(2)^2*x4p1(1) - x31(2)^2*x3p1(1)^2*x4p1(1) - x31(2)^2*x3p1(2)^2*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(1)^2 - x41(1)^2*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(2)^2 + x41(1)*x3p1(2)^2*x4p1(1)^2 - x41(1)^2*x3p1(2)^2*x4p1(1) - x41(2)^2*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(2)^2*x4p1(2)^2 - x41(2)^2*x3p1(2)^2*x4p1(1);
c(105) = x31(2)*x3p1(1)*x4p1(2) - x31(1)*x3p1(2)*x4p1(2) - x31(2)*x41(1)*x4p1(2) + x41(1)*x3p1(2)*x4p1(2);
c(106) = x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x4p1(1);
c(107) = x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1)^2 + x31(1)^2*x41(2)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1)*x4p1(2)^2;
c(108) = x31(1)^2*x41(1)^2*x3p1(1) - x31(1)*x41(1)^2*x3p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)^2 + x31(1)^2*x41(2)^2*x3p1(1) + x31(2)^2*x41(1)^2*x3p1(1) - x31(1)*x41(2)^2*x3p1(2)^2 + x31(2)^2*x41(2)^2*x3p1(1) - x31(1)*x41(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1)^2 - x31(1)*x41(2)^2*x4p1(2)^2 - x31(1)*x3p1(1)^2*x4p1(1)^2 + x31(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x3p1(2)^2*x4p1(1)^2 + x31(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x3p1(2)^2*x4p1(2)^2 + x31(2)^2*x3p1(1)*x4p1(2)^2 + x41(1)^2*x3p1(1)*x4p1(1)^2 + x41(1)^2*x3p1(1)*x4p1(2)^2 + x41(2)^2*x3p1(1)*x4p1(1)^2 + x41(2)^2*x3p1(1)*x4p1(2)^2;
c(109) = x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(2)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2)^2 - x31(1)^2*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2)^2 - x31(2)^2*x41(2)^2*x3p1(2)*x4p1(1);
c(110) = x31(2)*x41(1)*x3p1(2)*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(2)*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(2)*x4p1(1);
c(111) = x31(1)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)^2*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(1)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2)^2 - x31(1)^2*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(1)^2 - x31(2)^2*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(2)^2*x3p1(1)*x4p1(1);
c(112) = x31(1)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)*x4p1(1);
c(113) = x31(2)*x3p1(1) - x31(1)*x3p1(2) - x31(2)*x41(1) - x31(1)*x4p1(2) + x41(1)*x3p1(2) + x3p1(1)*x4p1(2);
c(114) = x41(1)*x3p1(1) - x31(1)*x4p1(1) - x31(1)*x41(1) + x3p1(1)*x4p1(1);
c(115) = x31(1)^2*x3p1(1) - x31(1)*x41(2)^2 - x31(1)*x3p1(1)^2 - x31(1)*x41(1)^2 - x31(1)*x3p1(2)^2 + x31(2)^2*x3p1(1) - x31(1)*x4p1(1)^2 + x41(1)^2*x3p1(1) - x31(1)*x4p1(2)^2 + x41(2)^2*x3p1(1) + x3p1(1)*x4p1(1)^2 + x3p1(1)*x4p1(2)^2;
c(116) = x31(1)*x41(1)*x3p1(2) + x31(2)*x41(1)*x3p1(1) - x31(1)*x3p1(2)*x4p1(1) - x31(2)*x3p1(1)*x4p1(1);
c(117) = x31(2)*x41(1)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(2) + x31(2)*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(2) + x31(2)*x41(1)*x4p1(1)^2 - x31(2)*x41(1)^2*x4p1(1) + x31(2)*x41(1)*x4p1(2)^2 - x31(2)*x41(2)^2*x4p1(1) - x31(2)*x3p1(1)^2*x4p1(1) - x31(1)^2*x3p1(2)*x4p1(1) - x31(2)*x3p1(2)^2*x4p1(1) - x31(2)^2*x3p1(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(1)^2 - x41(1)^2*x3p1(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(2)^2 - x41(2)^2*x3p1(2)*x4p1(1);
c(118) = x31(1)*x41(1)*x3p1(1)^2 + x31(1)^2*x41(1)*x3p1(1) + x31(1)*x41(1)*x3p1(2)^2 + x31(2)^2*x41(1)*x3p1(1) + x31(1)*x41(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(1) + x31(1)*x41(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(1) - x31(1)^2*x3p1(1)*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(1) - x31(2)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(1)^2 - x41(1)^2*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2)^2 - x41(2)^2*x3p1(1)*x4p1(1);
c(119) = x3p1(1) - x31(1);
c(120) = x31(2)*x41(1) - x31(2)*x4p1(1) + x41(1)*x3p1(2) - x3p1(2)*x4p1(1);
c(121) = x31(1)*x41(1) - x31(1)*x4p1(1) + x41(1)*x3p1(1) - x3p1(1)*x4p1(1);
c(122) = x31(1)^2*x41(1) + x31(2)^2*x41(1) + x41(1)*x3p1(1)^2 - x31(1)^2*x4p1(1) + x41(1)*x3p1(2)^2 - x31(2)^2*x4p1(1) + x41(1)*x4p1(1)^2 - x41(1)^2*x4p1(1) + x41(1)*x4p1(2)^2 - x41(2)^2*x4p1(1) - x3p1(1)^2*x4p1(1) - x3p1(2)^2*x4p1(1);
c(123) = x41(1) - x4p1(1);
c(124) = x31(2)*x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(2) + x31(2)*x41(2)*x3p1(1)*x4p1(1);
c(125) = x31(2)*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(2) + x31(2)*x41(2)*x3p1(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(2);
c(126) = x31(2)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(1);
c(127) = x31(2)*x41(1)*x3p1(2)*x4p1(2) - x31(2)*x41(2)*x3p1(2)*x4p1(1);
c(128) = x31(1)*x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(1)*x4p1(1);
c(129) = x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)^2*x4p1(1) - x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 + x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2) + x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2) + x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(2)^2*x4p1(1) - x31(2)*x41(2)^2*x3p1(1)^2*x4p1(1) - x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 + x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2) + x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(2)^2*x4p1(1);
c(130) = x31(2)*x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(2);
c(131) = x31(2)*x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(1);
c(132) = x31(2)*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(2)^2;
c(133) = x31(2)*x41(2)*x3p1(1) - x31(1)*x41(2)*x3p1(2) - x31(1)*x3p1(2)*x4p1(2) + x31(2)*x3p1(1)*x4p1(2);
c(134) = x31(2)*x41(1)*x3p1(1) - x31(1)*x41(1)*x3p1(2) - x31(1)*x3p1(2)*x4p1(1) + x31(2)*x3p1(1)*x4p1(1);
c(135) = x31(2)*x41(1)^2*x3p1(1) - x31(1)*x41(1)^2*x3p1(2) - x31(1)*x41(2)^2*x3p1(2) + x31(2)*x41(2)^2*x3p1(1) - x31(1)*x3p1(2)*x4p1(1)^2 + x31(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x3p1(2)*x4p1(2)^2 + x31(2)*x3p1(1)*x4p1(2)^2;
c(136) = x31(1)*x41(1)*x3p1(2)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(2) - x31(2)*x41(2)*x3p1(1)*x4p1(1);
c(137) = x31(1)*x41(1)^2*x3p1(2)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)*x41(2)*x3p1(1)^2*x4p1(1) - x31(2)*x41(1)^2*x3p1(2)*x4p1(1) - x31(1)*x41(2)*x3p1(2)*x4p1(2)^2 + x31(1)*x41(2)^2*x3p1(2)*x4p1(2) + x31(2)*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)*x41(2)*x3p1(2)^2*x4p1(1) - x31(2)*x41(2)^2*x3p1(2)*x4p1(1);
c(138) = x31(1)*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)*x41(1)^2*x3p1(1)*x4p1(2) + x31(2)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1) + x31(1)*x41(2)^2*x3p1(1)*x4p1(2) + x31(2)*x41(1)*x3p1(1)*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(1)*x4p1(1);
c(139) = x31(2)*x3p1(1) - x31(1)*x3p1(2);
c(140) = x31(2)*x41(1)*x3p1(2) - x31(1)*x41(2)*x3p1(2) + x31(2)*x41(1)*x4p1(2) - x31(2)*x41(2)*x4p1(1) + x31(1)*x3p1(2)*x4p1(2) - x31(2)*x3p1(2)*x4p1(1);
c(141) = x31(2)*x41(1)*x3p1(1) - x31(1)*x41(2)*x3p1(1) + x31(1)*x41(1)*x4p1(2) - x31(1)*x41(2)*x4p1(1) + x31(1)*x3p1(1)*x4p1(2) - x31(2)*x3p1(1)*x4p1(1);
c(142) = x31(2)*x41(1)*x3p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2 - x31(1)*x41(2)*x3p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2 - x31(1)*x41(2)*x4p1(1)^2 + x31(1)*x41(1)^2*x4p1(2) + x31(2)*x41(1)*x4p1(1)^2 - x31(2)*x41(1)^2*x4p1(1) - x31(1)*x41(2)*x4p1(2)^2 + x31(1)*x41(2)^2*x4p1(2) + x31(2)*x41(1)*x4p1(2)^2 - x31(2)*x41(2)^2*x4p1(1) + x31(1)*x3p1(1)^2*x4p1(2) - x31(2)*x3p1(1)^2*x4p1(1) + x31(1)*x3p1(2)^2*x4p1(2) - x31(2)*x3p1(2)^2*x4p1(1);
c(143) = x31(2)*x41(1) - x31(1)*x41(2) + x31(1)*x4p1(2) - x31(2)*x4p1(1);
c(144) = x31(2)*x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(2) + x31(2)*x41(2)*x3p1(1)*x4p1(1);
c(145) = x31(2)*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)^2*x41(1)*x3p1(2)*x4p1(2) + x31(1)^2*x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(2) + x31(2)*x41(1)*x3p1(2)^2*x4p1(2) + x31(2)*x41(2)*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(1)*x3p1(2)*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(2);
c(146) = x31(2)*x41(1)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(1) - x31(1)^2*x41(1)*x3p1(2)*x4p1(1) + x31(1)^2*x41(2)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(1) - x31(2)^2*x41(1)*x3p1(2)*x4p1(1) + x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(147) = x31(2)*x41(1)*x3p1(2)*x4p1(2) - x31(2)*x41(2)*x3p1(2)*x4p1(1);
c(148) = x31(1)*x41(1)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(1)*x4p1(1);
c(149) = x31(1)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)^2*x4p1(1);
c(150) = x31(2)*x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(2);
c(151) = x31(2)*x41(1)*x3p1(1)*x4p1(1) - x31(1)*x41(1)*x3p1(2)*x4p1(1);
c(152) = x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 - x31(1)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(1)^2*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(1)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(1)^2*x41(2)*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)^2*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)^2*x41(2)*x3p1(1)*x4p1(2)^2;
c(153) = x31(2)*x41(2)*x3p1(1) - x31(1)*x41(2)*x3p1(2) - x31(1)*x41(2)*x4p1(2) + x31(2)*x41(1)*x4p1(2) - x41(1)*x3p1(2)*x4p1(2) + x41(2)*x3p1(1)*x4p1(2);
c(154) = x31(2)*x41(1)*x3p1(1) - x31(1)*x41(1)*x3p1(2) - x31(1)*x41(2)*x4p1(1) + x31(2)*x41(1)*x4p1(1) - x41(1)*x3p1(2)*x4p1(1) + x41(2)*x3p1(1)*x4p1(1);
c(155) = x31(2)*x41(1)*x3p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2 - x31(1)^2*x41(1)*x3p1(2) + x31(1)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x3p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2 - x31(2)^2*x41(1)*x3p1(2) + x31(2)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x4p1(1)^2 + x31(2)*x41(1)*x4p1(1)^2 - x31(1)*x41(2)*x4p1(2)^2 + x31(2)*x41(1)*x4p1(2)^2 - x41(1)*x3p1(2)*x4p1(1)^2 + x41(2)*x3p1(1)*x4p1(1)^2 - x41(1)*x3p1(2)*x4p1(2)^2 + x41(2)*x3p1(1)*x4p1(2)^2;
c(156) = x31(1)*x41(1)*x3p1(2)*x4p1(2) - x31(1)*x41(2)*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(1)*x4p1(2) - x31(2)*x41(2)*x3p1(1)*x4p1(1);
c(157) = x31(2)*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)*x4p1(1);
c(158) = x31(1)*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)*x4p1(1) + x31(1)*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(159) = x31(2)*x41(1) - x31(1)*x41(2) - x41(1)*x3p1(2) + x41(2)*x3p1(1);
c(160) = x31(2)*x41(1)*x4p1(2) - x31(2)*x41(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(2) - x41(2)*x3p1(2)*x4p1(1);
c(161) = x31(1)*x41(1)*x4p1(2) - x31(1)*x41(2)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2) - x41(2)*x3p1(1)*x4p1(1);
c(162) = x31(1)^2*x41(1)*x4p1(2) - x31(1)^2*x41(2)*x4p1(1) + x31(2)^2*x41(1)*x4p1(2) - x31(2)^2*x41(2)*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(2) - x41(2)*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(2)^2*x4p1(2) - x41(2)*x3p1(2)^2*x4p1(1);
c(163) = x41(1)*x4p1(2) - x41(2)*x4p1(1);
c(164) = x31(1)*x41(1)^2*x3p1(2)*x4p1(2) - x31(1)*x41(2)*x3p1(1)^2*x4p1(2) + x31(2)*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)*x41(1)^2*x3p1(1)*x4p1(2) - x31(1)^2*x41(1)*x3p1(2)*x4p1(2) + x31(1)^2*x41(2)*x3p1(1)*x4p1(2) - x31(1)*x41(2)*x3p1(2)^2*x4p1(2) + x31(1)*x41(2)^2*x3p1(2)*x4p1(2) + x31(2)*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)*x41(2)^2*x3p1(1)*x4p1(2) - x31(2)^2*x41(1)*x3p1(2)*x4p1(2) + x31(2)^2*x41(2)*x3p1(1)*x4p1(2);
c(165) = x31(1)*x41(1)^2*x3p1(2)*x4p1(1) - x31(1)*x41(2)*x3p1(1)^2*x4p1(1) + x31(2)*x41(1)*x3p1(1)^2*x4p1(1) - x31(2)*x41(1)^2*x3p1(1)*x4p1(1) - x31(1)^2*x41(1)*x3p1(2)*x4p1(1) + x31(1)^2*x41(2)*x3p1(1)*x4p1(1) - x31(1)*x41(2)*x3p1(2)^2*x4p1(1) + x31(1)*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)*x41(1)*x3p1(2)^2*x4p1(1) - x31(2)*x41(2)^2*x3p1(1)*x4p1(1) - x31(2)^2*x41(1)*x3p1(2)*x4p1(1) + x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(166) = x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)^2*x4p1(2) - x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 + x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)^2*x4p1(2) - x31(1)*x41(2)^2*x3p1(1)^2*x4p1(2) - x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 - x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(2)^2*x4p1(1) + x31(2)*x41(2)^2*x3p1(1)^2*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)^2*x4p1(1) + x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)^2*x4p1(2) - x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(2)^2*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)^2*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)^2*x4p1(1);
c(167) = x31(1)*x41(1)^2*x3p1(2)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)*x4p1(1)^2 - x31(1)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(1)^2*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(1)^2*x4p1(2)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(1)^2 + x31(1)*x41(1)^2*x3p1(2)*x4p1(2)^2 + x31(1)*x41(2)^2*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(1)^2 - x31(2)*x41(1)^2*x3p1(1)*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(1)*x4p1(1)^2 - x31(1)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(1)^2*x41(2)*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)^2*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(2)*x3p1(2)^2*x4p1(2)^2 + x31(1)*x41(2)^2*x3p1(2)*x4p1(2)^2 + x31(2)*x41(1)*x3p1(2)^2*x4p1(2)^2 - x31(2)*x41(2)^2*x3p1(1)*x4p1(2)^2 - x31(2)^2*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)^2*x41(2)*x3p1(1)*x4p1(2)^2;
c(168) = x31(2)*x41(1)*x4p1(2) - x31(1)*x41(2)*x4p1(2) + x31(1)*x3p1(2)*x4p1(2) - x31(2)*x3p1(1)*x4p1(2) - x41(1)*x3p1(2)*x4p1(2) + x41(2)*x3p1(1)*x4p1(2);
c(169) = x31(2)*x41(1)*x4p1(1) - x31(1)*x41(2)*x4p1(1) + x31(1)*x3p1(2)*x4p1(1) - x31(2)*x3p1(1)*x4p1(1) - x41(1)*x3p1(2)*x4p1(1) + x41(2)*x3p1(1)*x4p1(1);
c(170) = x31(1)*x41(1)^2*x3p1(2) - x31(1)*x41(2)*x3p1(1)^2 + x31(2)*x41(1)*x3p1(1)^2 - x31(2)*x41(1)^2*x3p1(1) - x31(1)^2*x41(1)*x3p1(2) + x31(1)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x3p1(2)^2 + x31(1)*x41(2)^2*x3p1(2) + x31(2)*x41(1)*x3p1(2)^2 - x31(2)*x41(2)^2*x3p1(1) - x31(2)^2*x41(1)*x3p1(2) + x31(2)^2*x41(2)*x3p1(1) - x31(1)*x41(2)*x4p1(1)^2 + x31(2)*x41(1)*x4p1(1)^2 - x31(1)*x41(2)*x4p1(2)^2 + x31(2)*x41(1)*x4p1(2)^2 + x31(1)*x3p1(2)*x4p1(1)^2 - x31(2)*x3p1(1)*x4p1(1)^2 + x31(1)*x3p1(2)*x4p1(2)^2 - x31(2)*x3p1(1)*x4p1(2)^2 - x41(1)*x3p1(2)*x4p1(1)^2 + x41(2)*x3p1(1)*x4p1(1)^2 - x41(1)*x3p1(2)*x4p1(2)^2 + x41(2)*x3p1(1)*x4p1(2)^2;
c(171) = x31(1)*x41(2)*x3p1(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(2)*x4p1(2) - x31(2)*x41(1)*x3p1(2)*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(2)*x4p1(1) + x31(1)^2*x41(1)*x3p1(2)*x4p1(2) - x31(1)^2*x41(2)*x3p1(2)*x4p1(1) + x31(1)*x41(2)*x3p1(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(2)*x4p1(2) - x31(2)*x41(1)*x3p1(2)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(2)*x4p1(1) + x31(2)^2*x41(1)*x3p1(2)*x4p1(2) - x31(2)^2*x41(2)*x3p1(2)*x4p1(1);
c(172) = x31(1)*x41(2)*x3p1(1)*x4p1(1)^2 - x31(1)*x41(1)^2*x3p1(1)*x4p1(2) - x31(2)*x41(1)*x3p1(1)*x4p1(1)^2 + x31(2)*x41(1)^2*x3p1(1)*x4p1(1) + x31(1)^2*x41(1)*x3p1(1)*x4p1(2) - x31(1)^2*x41(2)*x3p1(1)*x4p1(1) + x31(1)*x41(2)*x3p1(1)*x4p1(2)^2 - x31(1)*x41(2)^2*x3p1(1)*x4p1(2) - x31(2)*x41(1)*x3p1(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x3p1(1)*x4p1(1) + x31(2)^2*x41(1)*x3p1(1)*x4p1(2) - x31(2)^2*x41(2)*x3p1(1)*x4p1(1);
c(173) = x31(2)*x41(1) - x31(1)*x41(2) + x31(1)*x3p1(2) - x31(2)*x3p1(1) - x41(1)*x3p1(2) + x41(2)*x3p1(1);
c(174) = x31(1)*x41(2)*x3p1(2) - x31(2)*x41(1)*x3p1(2) - x31(1)*x3p1(2)*x4p1(2) + x31(2)*x3p1(2)*x4p1(1) + x41(1)*x3p1(2)*x4p1(2) - x41(2)*x3p1(2)*x4p1(1);
c(175) = x31(1)*x41(2)*x3p1(1) - x31(2)*x41(1)*x3p1(1) - x31(1)*x3p1(1)*x4p1(2) + x31(2)*x3p1(1)*x4p1(1) + x41(1)*x3p1(1)*x4p1(2) - x41(2)*x3p1(1)*x4p1(1);
c(176) = x31(1)*x41(2)*x3p1(1)^2 - x31(2)*x41(1)*x3p1(1)^2 + x31(1)*x41(2)*x3p1(2)^2 - x31(2)*x41(1)*x3p1(2)^2 + x31(1)*x41(2)*x4p1(1)^2 - x31(1)*x41(1)^2*x4p1(2) - x31(2)*x41(1)*x4p1(1)^2 + x31(2)*x41(1)^2*x4p1(1) + x31(1)^2*x41(1)*x4p1(2) - x31(1)^2*x41(2)*x4p1(1) + x31(1)*x41(2)*x4p1(2)^2 - x31(1)*x41(2)^2*x4p1(2) - x31(2)*x41(1)*x4p1(2)^2 + x31(2)*x41(2)^2*x4p1(1) + x31(2)^2*x41(1)*x4p1(2) - x31(2)^2*x41(2)*x4p1(1) - x31(1)*x3p1(1)^2*x4p1(2) + x31(2)*x3p1(1)^2*x4p1(1) - x31(1)*x3p1(2)^2*x4p1(2) + x31(2)*x3p1(2)^2*x4p1(1) + x41(1)*x3p1(1)^2*x4p1(2) - x41(2)*x3p1(1)^2*x4p1(1) + x41(1)*x3p1(2)^2*x4p1(2) - x41(2)*x3p1(2)^2*x4p1(1);
c(177) = x31(1)*x41(2) - x31(2)*x41(1) - x31(1)*x4p1(2) + x31(2)*x4p1(1) + x41(1)*x4p1(2) - x41(2)*x4p1(1);
M = zeros(116, 125);
M(12289) = c(1);
M(13449) = c(2);
M(13101) = c(3);
M(12869) = c(4);
M(13913) = c(5);
M(12753) = c(6);
M(12405) = c(7);
M(12637) = c(8);
M(12057) = c(9);
M(13333) = c(10);
M(13797) = c(11);
M(13681) = c(12);
M(14261) = c(13);
M(14145) = c(14);
M(14377) = c(15);
M(14493) = c(16);
M(12290) = c(17);
M(13450) = c(18);
M(13102) = c(19);
M(12870) = c(20);
M(13914) = c(21);
M(12754) = c(22);
M(12406) = c(23);
M(12638) = c(24);
M(12058) = c(25);
M(13334) = c(26);
M(13798) = c(27);
M(13682) = c(28);
M(14262) = c(29);
M(14146) = c(30);
M(14378) = c(31);
M(14494) = c(32);
ci = [];
M(ci) = c(33);
ci = [];
M(ci) = c(34);
ci = [];
M(ci) = c(35);
ci = [];
M(ci) = c(36);
ci = [];
M(ci) = c(37);
ci = [];
M(ci) = c(38);
ci = [];
M(ci) = c(39);
ci = [];
M(ci) = c(40);
ci = [];
M(ci) = c(41);
ci = [];
M(ci) = c(42);
ci = [1165, 2803, 3500, 4454, 4571, 5268, 6197, 6430, 7127, 8201, 8318, 8435, 9248, 9481, 10178, 11471, 11588, 11705, 13451];
M(ci) = c(43);
ci = [933, 2455, 3152, 4222, 4339, 4804, 5733, 6082, 6779, 7737, 7854, 8087, 8668, 9133, 9830, 10891, 11008, 11357, 13103];
M(ci) = c(44);
ci = [1513, 3615, 3964, 5614, 5847, 6544, 7241, 7474, 7591, 9593, 9710, 9943, 10292, 10525, 10642, 12515, 12632, 12749, 13915];
M(ci) = c(45);
ci = [1049, 2571, 3384, 4338, 4455, 5036, 5965, 6314, 7011, 7969, 8086, 8319, 9016, 9365, 10062, 11239, 11356, 11589, 13335];
M(ci) = c(46);
ci = [1397, 3267, 3848, 5034, 5267, 5848, 6893, 7126, 7475, 8897, 9014, 9247, 9944, 10177, 10526, 12167, 12284, 12401, 13799];
M(ci) = c(47);
ci = [1281, 3035, 3732, 4802, 5035, 5616, 6661, 7010, 7359, 8549, 8666, 9015, 9712, 10061, 10410, 11935, 12052, 12285, 13683];
M(ci) = c(48);
ci = [3253, 5819, 6864, 7934, 8167, 8864, 9909, 10142, 10491, 11101, 11218, 11451, 12148, 12381, 12730, 13211, 13328, 13445, 14263];
M(ci) = c(49);
ci = [3021, 5587, 6632, 7702, 7935, 8516, 9677, 10026, 10375, 10753, 10870, 11219, 11916, 12265, 12614, 12979, 13096, 13329, 14147];
M(ci) = c(50);
ci = [3601, 6515, 7212, 8514, 8863, 9560, 10257, 10490, 10607, 11797, 11914, 12147, 12496, 12729, 12846, 13559, 13676, 13793, 14379];
M(ci) = c(51);
ci = [6501, 9531, 10228, 10718, 11067, 11764, 12461, 12694, 12811, 12958, 13191, 13540, 13773, 13890, 14023, 14140, 14257, 14495];
M(ci) = c(52);
ci = [1398, 1515, 2573, 2806, 3271, 3620, 3853, 3970, 5620, 5853, 5970, 6203, 6552, 6901, 7018, 7135, 7252, 7369, 7486, 7603, 9599, 9716, 9833, 9950, 10067, 10184, 10301, 10418, 10535, 10652, 12518, 12635, 12868, 13916];
M(ci) = c(53);
ci = [1166, 1399, 2109, 2226, 2807, 3272, 3505, 3854, 5040, 5273, 5390, 5507, 5856, 6205, 6322, 6439, 6904, 7021, 7138, 7487, 8903, 9020, 9137, 9254, 9371, 9488, 9953, 10070, 10187, 10536, 12170, 12287, 12752, 13800];
M(ci) = c(54);
ci = [1050, 1283, 1993, 2110, 2575, 3040, 3389, 3738, 4808, 5041, 5158, 5391, 5624, 5973, 6090, 6323, 6672, 6789, 7022, 7371, 8555, 8672, 8789, 9022, 9139, 9372, 9721, 9838, 10071, 10420, 11938, 12055, 12636, 13684];
M(ci) = c(55);
ci = [2790, 3255, 4429, 4546, 5243, 5824, 6173, 6870, 7940, 8173, 8290, 8407, 8872, 9221, 9338, 9455, 9920, 10037, 10154, 10503, 11107, 11224, 11341, 11458, 11575, 11692, 12157, 12274, 12391, 12740, 13214, 13331, 13796, 14264];
M(ci) = c(56);
ci = [2558, 3023, 4313, 4430, 5011, 5592, 5941, 6638, 7708, 7941, 8058, 8291, 8524, 8989, 9106, 9339, 9688, 9805, 10038, 10387, 10759, 10876, 10993, 11226, 11343, 11576, 11925, 12042, 12275, 12624, 12982, 13099, 13680, 14148];
M(ci) = c(57);
ci = [3254, 3603, 5009, 5242, 5823, 6520, 6869, 7218, 8520, 8869, 8986, 9219, 9568, 9917, 10034, 10151, 10268, 10385, 10502, 10619, 11803, 11920, 12037, 12154, 12271, 12388, 12505, 12622, 12739, 12856, 13562, 13679, 13912, 14380];
M(ci) = c(58);
ci = [5806, 6503, 7909, 8142, 8839, 9536, 9885, 10234, 10724, 11073, 11190, 11423, 11772, 12121, 12238, 12355, 12472, 12589, 12706, 12823, 12964, 13081, 13198, 13315, 13432, 13549, 13666, 13783, 13900, 14026, 14143, 14376, 14496];
M(ci) = c(59);
M(9045) = c(60);
M(11249) = c(61);
M(12177) = c(62);
M(11945) = c(63);
M(12293) = c(64);
M(9973) = c(65);
M(9277) = c(66);
M(9741) = c(67);
M(8697) = c(68);
M(13453) = c(69);
M(13105) = c(70);
M(12873) = c(71);
M(13917) = c(72);
M(11481) = c(73);
M(10901) = c(74);
M(10321) = c(75);
M(12525) = c(76);
M(12757) = c(77);
M(12409) = c(78);
M(12641) = c(79);
M(12061) = c(80);
M(13221) = c(81);
M(12989) = c(82);
M(13569) = c(83);
M(13337) = c(84);
M(13801) = c(85);
M(13685) = c(86);
M(14033) = c(87);
M(14265) = c(88);
M(14149) = c(89);
M(14381) = c(90);
M(14497) = c(91);
M(9046) = c(92);
M(11250) = c(93);
M(12178) = c(94);
M(11946) = c(95);
M(12294) = c(96);
M(9974) = c(97);
M(9278) = c(98);
M(9742) = c(99);
M(8698) = c(100);
M(13454) = c(101);
M(13106) = c(102);
M(12874) = c(103);
M(13918) = c(104);
M(11482) = c(105);
M(10902) = c(106);
M(10322) = c(107);
M(12526) = c(108);
M(12758) = c(109);
M(12410) = c(110);
M(12642) = c(111);
M(12062) = c(112);
M(13222) = c(113);
M(12990) = c(114);
M(13570) = c(115);
M(13338) = c(116);
M(13802) = c(117);
M(13686) = c(118);
M(14034) = c(119);
M(14266) = c(120);
M(14150) = c(121);
M(14382) = c(122);
M(14498) = c(123);
ci = [];
M(ci) = c(124);
ci = [];
M(ci) = c(125);
ci = [];
M(ci) = c(126);
ci = [];
M(ci) = c(127);
ci = [];
M(ci) = c(128);
ci = [];
M(ci) = c(129);
ci = [];
M(ci) = c(130);
ci = [];
M(ci) = c(131);
ci = [];
M(ci) = c(132);
ci = [];
M(ci) = c(133);
ci = [];
M(ci) = c(134);
ci = [];
M(ci) = c(135);
ci = [];
M(ci) = c(136);
ci = [];
M(ci) = c(137);
ci = [];
M(ci) = c(138);
ci = [];
M(ci) = c(139);
ci = [];
M(ci) = c(140);
ci = [];
M(ci) = c(141);
ci = [];
M(ci) = c(142);
ci = [];
M(ci) = c(143);
ci = [117, 472, 705, 1054, 1884, 2117, 2582, 2931, 3396, 4357, 4474, 5055, 5172, 5405, 5986, 7985, 8102, 8335, 11251];
M(ci) = c(144);
ci = [1402, 2813, 3278, 3511, 3860, 5053, 5286, 5867, 5984, 6217, 6914, 8913, 9030, 9263, 12179];
M(ci) = c(145);
ci = [233, 704, 821, 1286, 1883, 2348, 2581, 3046, 3395, 3744, 4704, 4821, 5054, 5635, 5752, 5985, 6682, 8565, 8682, 9031, 11947];
M(ci) = c(146);
ci = [1161, 2212, 2793, 3490, 4551, 5248, 5481, 6178, 6411, 7108, 8184, 8301, 8418, 9231, 9348, 9465, 10162, 11465, 11582, 11699, 13455];
M(ci) = c(147);
ci = [929, 1980, 2445, 3142, 4319, 4784, 5133, 5714, 6063, 6760, 7720, 7837, 8070, 8651, 8768, 9117, 9814, 10885, 11002, 11351, 13107];
M(ci) = c(148);
ci = [1509, 3256, 3605, 3954, 5827, 6524, 6873, 7222, 7455, 7572, 9576, 9693, 9926, 10275, 10392, 10509, 10626, 12509, 12626, 12743, 13919];
M(ci) = c(149);
ci = [1170, 2233, 2814, 3512, 4473, 4590, 5287, 5404, 5521, 6218, 8217, 8334, 8451, 11483];
M(ci) = c(150);
ci = [1, 356, 589, 938, 1651, 1768, 2001, 2466, 2699, 3164, 4124, 4241, 4358, 4823, 4940, 5173, 5754, 7753, 7870, 8103, 10903];
M(ci) = c(151);
ci = [1518, 3277, 3626, 3859, 3976, 5633, 5866, 6563, 6680, 6913, 7262, 9609, 9726, 9959, 12527];
M(ci) = c(152);
ci = [3258, 5249, 5830, 6179, 6876, 7953, 8186, 8883, 9000, 9233, 9930, 11117, 11234, 11467, 13223];
M(ci) = c(153);
ci = [813, 1864, 2329, 3026, 4668, 5017, 5598, 5947, 6644, 7721, 7954, 8535, 8652, 9001, 9698, 10769, 10886, 11235, 12991];
M(ci) = c(154);
ci = [3606, 5829, 6526, 6875, 7224, 8533, 8882, 9579, 9696, 9929, 10278, 11813, 11930, 12163, 13571];
M(ci) = c(155);
ci = [1045, 2096, 2561, 3374, 4435, 5016, 5365, 5946, 6295, 6992, 7952, 8069, 8302, 8999, 9116, 9349, 10046, 11233, 11350, 11583, 13339];
M(ci) = c(156);
ci = [1393, 2792, 3257, 3838, 5247, 5828, 6177, 6874, 7107, 7456, 8880, 8997, 9230, 9927, 10044, 10161, 10510, 12161, 12278, 12395, 13803];
M(ci) = c(157);
ci = [1277, 2560, 3025, 3722, 5015, 5596, 5945, 6642, 6991, 7340, 8532, 8649, 8998, 9695, 9812, 10045, 10394, 11929, 12046, 12279, 13687];
M(ci) = c(158);
ci = [6506, 8845, 9542, 9891, 10240, 10737, 11086, 11783, 11900, 12133, 12482, 12974, 13207, 14035];
M(ci) = c(159);
ci = [3249, 5228, 5809, 6854, 8147, 8844, 9193, 9890, 10123, 10472, 11084, 11201, 11434, 12131, 12248, 12365, 12714, 13205, 13322, 13439, 14267];
M(ci) = c(160);
ci = [3017, 4996, 5577, 6622, 7915, 8496, 8961, 9658, 10007, 10356, 10736, 10853, 11202, 11899, 12016, 12249, 12598, 12973, 13090, 13323, 14151];
M(ci) = c(161);
ci = [3597, 5808, 6505, 7202, 8843, 9540, 9889, 10238, 10471, 10588, 11780, 11897, 12130, 12479, 12596, 12713, 12830, 13553, 13670, 13787, 14383];
M(ci) = c(162);
ci = [6497, 8824, 9521, 10218, 11047, 11744, 12093, 12442, 12675, 12792, 12941, 13174, 13523, 13640, 13757, 13874, 14017, 14134, 14251, 14499];
M(ci) = c(163);
ci = [118, 475, 709, 1060, 1177, 1410, 1889, 2123, 2240, 2590, 2823, 2940, 3290, 3407, 3524, 3873, 5060, 5177, 5294, 5411, 5528, 5877, 5994, 6227, 6924, 8916, 9033, 9962, 12180];
M(ci) = c(164);
ci = [2, 119, 236, 359, 476, 593, 710, 827, 944, 1061, 1294, 1773, 1890, 2007, 2124, 2357, 2474, 2591, 2708, 2941, 3058, 3175, 3408, 3757, 4711, 4828, 4945, 5062, 5179, 5412, 5645, 5762, 5995, 6692, 8568, 8685, 9730, 11948];
M(ci) = c(165);
ci = [1278, 1395, 1512, 2563, 2796, 3029, 3262, 3611, 3728, 3845, 3962, 5601, 5834, 5951, 6184, 6533, 6650, 6883, 7000, 7117, 7234, 7351, 7468, 7585, 9583, 9700, 9817, 9934, 10051, 10168, 10285, 10402, 10519, 10636, 12512, 12629, 12862, 13920];
M(ci) = c(166);
ci = [234, 707, 825, 1292, 1409, 1526, 2353, 2587, 2820, 3054, 3287, 3404, 3521, 3638, 3755, 3872, 3989, 5640, 5757, 5874, 5991, 6224, 6573, 6690, 6923, 7272, 9612, 9729, 10310, 12528];
M(ci) = c(167);
ci = [698, 1869, 2568, 2801, 3266, 4443, 4560, 5026, 5259, 5376, 5493, 5842, 5959, 6192, 6889, 7960, 8077, 8194, 8311, 8428, 8893, 9010, 9243, 9940, 11120, 11237, 12166, 13224];
M(ci) = c(168);
ci = [582, 699, 816, 1635, 1753, 1870, 2335, 2452, 2569, 3034, 4093, 4327, 4444, 4677, 4794, 5027, 5144, 5377, 5610, 5727, 5960, 6657, 7728, 7845, 7962, 8079, 8312, 8545, 8662, 9011, 9708, 10772, 10889, 11934, 12992];
M(ci) = c(169);
ci = [814, 1867, 2333, 3032, 3265, 3614, 4673, 5023, 5256, 5606, 5839, 5956, 6189, 6538, 6655, 6888, 7237, 8540, 8657, 8890, 9007, 9240, 9589, 9706, 9939, 10288, 11816, 11933, 12514, 13572];
M(ci) = c(170);
ci = [1046, 1163, 1396, 2099, 2216, 2565, 2798, 3263, 3380, 3497, 3846, 5021, 5254, 5371, 5488, 5837, 5954, 6187, 6304, 6421, 6886, 7003, 7120, 7469, 8887, 9004, 9121, 9238, 9355, 9472, 9937, 10054, 10171, 10520, 12164, 12281, 12746, 13804];
M(ci) = c(171);
ci = [930, 1047, 1280, 1983, 2100, 2449, 2566, 3031, 3148, 3381, 3730, 4789, 5022, 5139, 5372, 5605, 5722, 5955, 6072, 6305, 6654, 6771, 7004, 7353, 8539, 8656, 8773, 9006, 9123, 9356, 9705, 9822, 10055, 10404, 11932, 12049, 12630, 13688];
M(ci) = c(172);
ci = [2322, 4653, 5584, 5817, 6514, 7923, 8156, 8506, 8855, 8972, 9205, 9554, 9671, 9904, 10253, 10744, 10861, 11094, 11211, 11444, 11793, 11910, 12143, 12492, 12977, 13558, 14036];
M(ci) = c(173);
ci = [2554, 2787, 3252, 4419, 4536, 5001, 5234, 5815, 5932, 6165, 6862, 7921, 8154, 8271, 8388, 8853, 8970, 9203, 9320, 9437, 9902, 10019, 10136, 10485, 11091, 11208, 11325, 11442, 11559, 11676, 12141, 12258, 12375, 12724, 13208, 13325, 13790, 14268];
M(ci) = c(174);
ci = [2438, 2555, 3020, 4303, 4420, 4769, 5002, 5583, 5700, 5933, 6630, 7689, 7922, 8039, 8272, 8505, 8622, 8971, 9088, 9321, 9670, 9787, 10020, 10369, 10743, 10860, 10977, 11210, 11327, 11560, 11909, 12026, 12259, 12608, 12976, 13093, 13674, 14152];
M(ci) = c(175);
ci = [3018, 3251, 3600, 4999, 5232, 5581, 5814, 6511, 6628, 6861, 7210, 8501, 8850, 8967, 9200, 9549, 9666, 9899, 10016, 10133, 10250, 10367, 10484, 10601, 11787, 11904, 12021, 12138, 12255, 12372, 12489, 12606, 12723, 12840, 13556, 13673, 13906, 14384];
M(ci) = c(176);
ci = [5570, 5803, 6500, 7899, 8132, 8481, 8830, 9527, 9644, 9877, 10226, 10705, 11054, 11171, 11404, 11753, 11870, 12103, 12220, 12337, 12454, 12571, 12688, 12805, 12948, 13065, 13182, 13299, 13416, 13533, 13650, 13767, 13884, 14020, 14137, 14370, 14500];
M(ci) = c(177);
M = gj(M);
A = zeros(9);
amcols = [125 124 123 122 121 120 119 117 116];
A(1, 4) = 1;
A(2, :) = -M(116, amcols);
A(3, :) = -M(115, amcols);
A(4, :) = -M(113, amcols);
A(5, :) = -M(112, amcols);
A(6, :) = -M(109, amcols);
A(7, :) = -M(106, amcols);
A(8, :) = -M(103, amcols);
A(9, :) = -M(100, amcols);
[V D] = eig(A);
sol = V([5, 4, 3, 2],:)./(ones(4, 1)*V(1,:));
% if (find(isnan(sol(:))) > 0)
%
% s = zeros(1, 0);
% l1 = zeros(1, 0);
% l2 = zeros(1, 0);
% k = zeros(1, 0);
% else
I = find(not(imag( sol(1,:) )));
if isempty(I) == 0
s = sol(1,I);
l1 = sol(2,I);
l2 = sol(3,I);
k = sol(4,I);
for i = 1:size(l1,2)
M1 = [ - x11(2)*(k(i)*(x1p1(1)^2 + x1p1(2)^2) + 1) - l1(i)*x11(2)*x1p1(1) - l2(i)*x11(2)*x1p1(2), x11(1)*(k(i)*(x1p1(1)^2 + x1p1(2)^2) + 1) + l1(i)*x11(1)*x1p1(1) + l2(i)*x11(1)*x1p1(2), 0, x11(1)*x1p1(2) - x11(2)*x1p1(1); ...
- x21(2)*(k(i)*(x2p1(1)^2 + x2p1(2)^2) + 1) - l1(i)*x21(2)*x2p1(1) - l2(i)*x21(2)*x2p1(2), x21(1)*(k(i)*(x2p1(1)^2 + x2p1(2)^2) + 1) + l1(i)*x21(1)*x2p1(1) + l2(i)*x21(1)*x2p1(2), 0, x21(1)*x2p1(2) - x21(2)*x2p1(1) ; ...
(k(i)*(x11(1)^2 + x11(2)^2) + 1)*(k(i)*(x1p1(1)^2 + x1p1(2)^2) + 1) + l1(i)*x1p1(1)*(k(i)*(x11(1)^2 + x11(2)^2) + 1) + l2(i)*x1p1(2)*(k(i)*(x11(1)^2 + x11(2)^2) + 1), 0, - x11(1)*(k(i)*(x1p1(1)^2 + x1p1(2)^2) + 1) - l1(i)*x11(1)*x1p1(1) - l2(i)*x11(1)*x1p1(2), x1p1(1)*(k(i)*(x11(1)^2 + x11(2)^2) + 1) - x11(1)*(k(i)*(x1p1(1)^2 + x1p1(2)^2) + 1) ; ...
(k(i)*(x21(1)^2 + x21(2)^2) + 1)*(k(i)*(x2p1(1)^2 + x2p1(2)^2) + 1) + l1(i)*x2p1(1)*(k(i)*(x21(1)^2 + x21(2)^2) + 1) + l2(i)*x2p1(2)*(k(i)*(x21(1)^2 + x21(2)^2) + 1), 0, - x21(1)*(k(i)*(x2p1(1)^2 + x2p1(2)^2) + 1) - l1(i)*x21(1)*x2p1(1) - l2(i)*x21(1)*x2p1(2), x2p1(1)*(k(i)*(x21(1)^2 + x21(2)^2) + 1) - x21(1)*(k(i)*(x2p1(1)^2 + x2p1(2)^2) + 1) ; ...
l1(i), l2(i), 1, 0 ];
[~,~,V] = svd(M1);
v1(i) = V(1,end)/V(4,end);
v2(i) = V(2,end)/V(4,end);
v3(i) = V(3,end)/V(4,end);
M2 = [ - x31(2)*(k(i)*(x3p1(1)^2 + x3p1(2)^2) + 1) - l1(i)*x31(2)*x3p1(1) - l2(i)*x31(2)*x3p1(2), x31(1)*(k(i)*(x3p1(1)^2 + x3p1(2)^2) + 1) + l1(i)*x31(1)*x3p1(1) + l2(i)*x31(1)*x3p1(2), 0, x31(1)*x3p1(2) - x31(2)*x3p1(1); ...
- s(i)*x41(2)*(k(i)*(x4p1(1)^2 + x4p1(2)^2) + 1) - l1(i)*s(i)*x41(2)*x4p1(1) - l2(i)*s(i)*x41(2)*x4p1(2), s(i)*x41(1)*(k(i)*(x4p1(1)^2 + x4p1(2)^2) + 1) + l1(i)*s(i)*x41(1)*x4p1(1) + l2(i)*s(i)*x41(1)*x4p1(2), 0, x41(1)*x4p1(2) - x41(2)*x4p1(1); ...
(k(i)*(x31(1)^2 + x31(2)^2) + 1)*(k(i)*(x3p1(1)^2 + x3p1(2)^2) + 1) + l1(i)*x3p1(1)*(k(i)*(x31(1)^2 + x31(2)^2) + 1) + l2(i)*x3p1(2)*(k(i)*(x31(1)^2 + x31(2)^2) + 1), 0, - x31(1)*(k(i)*(x3p1(1)^2 + x3p1(2)^2) + 1) - l1(i)*x31(1)*x3p1(1) - l2(i)*x31(1)*x3p1(2), x3p1(1)*(k(i)*(x31(1)^2 + x31(2)^2) + 1) - x31(1)*(k(i)*(x3p1(1)^2 + x3p1(2)^2) + 1); ...
s(i)*(k(i)*(x41(1)^2 + x41(2)^2) + 1)*(k(i)*(x4p1(1)^2 + x4p1(2)^2) + 1) + l1(i)*s(i)*x4p1(1)*(k(i)*(x41(1)^2 + x41(2)^2) + 1) + l2(i)*s(i)*x4p1(2)*(k(i)*(x41(1)^2 + x41(2)^2) + 1), 0, - s(i)*x41(1)*(k(i)*(x4p1(1)^2 + x4p1(2)^2) + 1) - l1(i)*s(i)*x41(1)*x4p1(1) - l2(i)*s(i)*x41(1)*x4p1(2), x4p1(1)*(k(i)*(x41(1)^2 + x41(2)^2) + 1) - x41(1)*(k(i)*(x4p1(1)^2 + x4p1(2)^2) + 1); ...
l1(i), l2(i), 1, 0];
[~,~,U] = svd(M2);
u1(i) = U(1,end)/U(4,end);
u2(i) = U(2,end)/U(4,end);
u3(i) = U(3,end)/U(4,end);
end
end
end