forked from qucontrol/krotov
-
Notifications
You must be signed in to change notification settings - Fork 1
/
gaussian.py
2984 lines (2936 loc) · 258 KB
/
gaussian.py
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
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
import numpy as np
from IPython.display import display, Markdown
def separate(V,row):
return V[row][0],V[row][1],V[row][2],V[row][3]
def coeffscuadraticos11(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(v13*(v24**2))+((0.5*(v14*v34))+((v14*(v22*v34))+((-0.5*(v13*v44))+(v12*(v23*v44)))));
aux1=((-0.25*v13)+((-0.5*(v13*v22))+((0.5*(v12*v23))+aux0)))-(v13*(v22*v44));
aux2=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux3=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux2));
aux4=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux3)));
aux5=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux4)));
aux6=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux5)));
aux7=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux6));
aux8=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux7));
aux9=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux8)));
aux10=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux9)));
aux11=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux10));
aux12=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux11)));
aux13=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux12)));
aux14=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux13))))-((v12**2)*(v33*v44));
aux15=((aux14-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux16=((0.+-1.j)*((aux1-(v12*(v24*v34)))-(v14*(v23*v24))))/((aux15-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux17=(-0.5*(v34**2))+((0.25*v44)+((0.5*(v22*v44))+((0.5*(v33*v44))+(v22*(v33*v44)))));
aux18=(-0.5*(v24**2))+((0.25*v33)+((0.5*(v22*v33))+((2.*(v23*(v24*v34)))+aux17)));
aux19=((0.125+((0.25*v22)+((-0.5*(v23**2))+aux18)))-((v23**2)*v44))-(v22*(v34**2));
aux20=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux21=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux20));
aux22=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux21)));
aux23=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux22)));
aux24=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux23)));
aux25=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux24));
aux26=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux25));
aux27=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux26)));
aux28=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux27)));
aux29=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux28));
aux30=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux29)));
aux31=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux30)));
aux32=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux31))))-((v12**2)*(v33*v44));
aux33=((aux32-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux34=(aux19-((v24**2)*v33))/((aux33-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux35=(b13*(b24**2))+((0.5*(b14*b34))+((b14*(b22*b34))+((-0.5*(b13*b44))+(b12*(b23*b44)))));
aux36=((-0.25*b13)+((-0.5*(b13*b22))+((0.5*(b12*b23))+aux35)))-(b13*(b22*b44));
aux37=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux38=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux37));
aux39=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux38)));
aux40=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux39)));
aux41=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux40)));
aux42=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux41));
aux43=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux42));
aux44=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux43)));
aux45=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux44)));
aux46=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux45));
aux47=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux46)));
aux48=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux47)));
aux49=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux48))))-((b12**2)*(b33*b44));
aux50=((aux49-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux51=((0.+1.j)*((aux36-(b12*(b24*b34)))-(b14*(b23*b24))))/((aux50-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux52=(-0.5*(b34**2))+((0.25*b44)+((0.5*(b22*b44))+((0.5*(b33*b44))+(b22*(b33*b44)))));
aux53=(-0.5*(b24**2))+((0.25*b33)+((0.5*(b22*b33))+((2.*(b23*(b24*b34)))+aux52)));
aux54=((0.125+((0.25*b22)+((-0.5*(b23**2))+aux53)))-((b23**2)*b44))-(b22*(b34**2));
aux55=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux56=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux55));
aux57=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux56)));
aux58=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux57)));
aux59=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux58)));
aux60=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux59));
aux61=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux60));
aux62=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux61)));
aux63=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux62)));
aux64=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux63));
aux65=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux64)));
aux66=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux65)));
aux67=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux66))))-((b12**2)*(b33*b44));
aux68=((aux67-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux69=(aux54-((b24**2)*b33))/((aux68-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
output=(0.25*(aux16-aux34))+(-0.25*(np.conjugate((aux51+aux69))));
return output
def coeffscuadraticos12(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(0.5*(b12*b24))+((-0.5*(b14*b33))+((b12*(b24*b33))+((0.5*(b13*b34))+(b13*(b22*b34)))));
aux1=(((-0.25*b14)+((-0.5*(b14*b22))+((b14*(b23**2))+aux0)))-(b12*(b23*b34)))-(b14*(b22*b33));
aux2=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux3=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux2));
aux4=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux3)));
aux5=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux4)));
aux6=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux5)));
aux7=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux6));
aux8=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux7));
aux9=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux8)));
aux10=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux9)));
aux11=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux10));
aux12=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux11)));
aux13=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux12)));
aux14=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux13))))-((b12**2)*(b33*b44));
aux15=((aux14-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux16=(aux1-(b13*(b23*b24)))/((aux15-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux17=((b14**2)*b23)+((0.5*(b24*b34))+((b11*(b24*b34))+((b12*(b13*b44))+(-0.5*(b23*b44)))));
aux18=((0.5*(b12*b13))+((-0.25*b23)+((-0.5*(b11*b23))+aux17)))-(b11*(b23*b44));
aux19=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux20=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux19));
aux21=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux20)));
aux22=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux21)));
aux23=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux22)));
aux24=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux23));
aux25=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux24));
aux26=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux25)));
aux27=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux26)));
aux28=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux27));
aux29=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux28)));
aux30=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux29)));
aux31=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux30))))-((b12**2)*(b33*b44));
aux32=((aux31-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux33=((aux18-(b12*(b14*b34)))-(b13*(b14*b24)))/((aux32-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux34=(-0.5*(b12*b33))+((b14*(b24*b33))+((b12*(b34**2))+((-0.5*(b12*b44))+(b13*(b23*b44)))));
aux35=((-0.25*b12)+((0.5*(b13*b23))+((0.5*(b14*b24))+aux34)))-(b12*(b33*b44));
aux36=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux37=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux36));
aux38=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux37)));
aux39=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux38)));
aux40=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux39)));
aux41=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux40));
aux42=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux41));
aux43=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux42)));
aux44=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux43)));
aux45=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux44));
aux46=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux45)));
aux47=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux46)));
aux48=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux47))))-((b12**2)*(b33*b44));
aux49=((aux48-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux50=((aux35-(b13*(b24*b34)))-(b14*(b23*b34)))/((aux49-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux51=(0.5*(v12*v24))+((-0.5*(v14*v33))+((v12*(v24*v33))+((0.5*(v13*v34))+(v13*(v22*v34)))));
aux52=((-0.25*v14)+((-0.5*(v14*v22))+((v14*(v23**2))+aux51)))-(v12*(v23*v34));
aux53=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux54=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux53));
aux55=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux54)));
aux56=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux55)));
aux57=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux56)));
aux58=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux57));
aux59=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux58));
aux60=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux59)));
aux61=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux60)));
aux62=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux61));
aux63=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux62)));
aux64=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux63)));
aux65=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux64))))-((v12**2)*(v33*v44));
aux66=((aux65-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux67=((aux52-(v14*(v22*v33)))-(v13*(v23*v24)))/((aux66-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux68=((v14**2)*v23)+((0.5*(v24*v34))+((v11*(v24*v34))+((v12*(v13*v44))+(-0.5*(v23*v44)))));
aux69=((0.5*(v12*v13))+((-0.25*v23)+((-0.5*(v11*v23))+aux68)))-(v11*(v23*v44));
aux70=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux71=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux70));
aux72=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux71)));
aux73=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux72)));
aux74=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux73)));
aux75=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux74));
aux76=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux75));
aux77=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux76)));
aux78=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux77)));
aux79=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux78));
aux80=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux79)));
aux81=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux80)));
aux82=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux81))))-((v12**2)*(v33*v44));
aux83=((aux82-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux84=((aux69-(v12*(v14*v34)))-(v13*(v14*v24)))/((aux83-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux85=(-0.5*(v12*v33))+((v14*(v24*v33))+((v12*(v34**2))+((-0.5*(v12*v44))+(v13*(v23*v44)))));
aux86=((-0.25*v12)+((0.5*(v13*v23))+((0.5*(v14*v24))+aux85)))-(v12*(v33*v44));
aux87=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux88=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux87));
aux89=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux88)));
aux90=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux89)));
aux91=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux90)));
aux92=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux91));
aux93=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux92));
aux94=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux93)));
aux95=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux94)));
aux96=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux95));
aux97=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux96)));
aux98=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux97)));
aux99=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux98))))-((v12**2)*(v33*v44));
aux100=((aux99-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux101=((aux86-(v13*(v24*v34)))-(v14*(v23*v34)))/((aux100-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
output=(0.5*(((0.+0.5j)*(aux16+aux33))-aux50))+(0.5*(((0.+-0.5j)*(aux67+aux84))-aux101));
return output
def coeffscuadraticos13(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(-0.5*(b34**2))+((0.25*b44)+((0.5*(b22*b44))+((0.5*(b33*b44))+(b22*(b33*b44)))));
aux1=(-0.5*(b24**2))+((0.25*b33)+((0.5*(b22*b33))+((2.*(b23*(b24*b34)))+aux0)));
aux2=((0.125+((0.25*b22)+((-0.5*(b23**2))+aux1)))-((b23**2)*b44))-(b22*(b34**2));
aux3=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux4=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux3));
aux5=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux4)));
aux6=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux5)));
aux7=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux6)));
aux8=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux7));
aux9=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux8));
aux10=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux9)));
aux11=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux10)));
aux12=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux11));
aux13=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux12)));
aux14=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux13)));
aux15=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux14))))-((b12**2)*(b33*b44));
aux16=((aux15-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux17=(aux2-((b24**2)*b33))/((aux16-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux18=(-0.5*(b24**2))+((0.25*b44)+((0.5*(b11*b44))+((0.5*(b22*b44))+(b11*(b22*b44)))));
aux19=(-0.5*(b14**2))+((0.25*b22)+((0.5*(b11*b22))+((2.*(b12*(b14*b24)))+aux18)));
aux20=((0.125+((0.25*b11)+((-0.5*(b12**2))+aux19)))-((b12**2)*b44))-(b11*(b24**2));
aux21=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux22=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux21));
aux23=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux22)));
aux24=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux23)));
aux25=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux24)));
aux26=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux25));
aux27=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux26));
aux28=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux27)));
aux29=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux28)));
aux30=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux29));
aux31=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux30)));
aux32=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux31)));
aux33=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux32))))-((b12**2)*(b33*b44));
aux34=((aux33-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux35=(aux20-((b14**2)*b22))/((aux34-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux36=(b13*(b24**2))+((0.5*(b14*b34))+((b14*(b22*b34))+((-0.5*(b13*b44))+(b12*(b23*b44)))));
aux37=((-0.25*b13)+((-0.5*(b13*b22))+((0.5*(b12*b23))+aux36)))-(b13*(b22*b44));
aux38=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux39=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux38));
aux40=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux39)));
aux41=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux40)));
aux42=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux41)));
aux43=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux42));
aux44=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux43));
aux45=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux44)));
aux46=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux45)));
aux47=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux46));
aux48=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux47)));
aux49=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux48)));
aux50=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux49))))-((b12**2)*(b33*b44));
aux51=((aux50-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux52=((aux37-(b12*(b24*b34)))-(b14*(b23*b24)))/((aux51-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux53=(-0.5*(v34**2))+((0.25*v44)+((0.5*(v22*v44))+((0.5*(v33*v44))+(v22*(v33*v44)))));
aux54=(-0.5*(v24**2))+((0.25*v33)+((0.5*(v22*v33))+((2.*(v23*(v24*v34)))+aux53)));
aux55=((0.125+((0.25*v22)+((-0.5*(v23**2))+aux54)))-((v23**2)*v44))-(v22*(v34**2));
aux56=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux57=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux56));
aux58=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux57)));
aux59=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux58)));
aux60=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux59)));
aux61=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux60));
aux62=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux61));
aux63=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux62)));
aux64=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux63)));
aux65=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux64));
aux66=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux65)));
aux67=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux66)));
aux68=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux67))))-((v12**2)*(v33*v44));
aux69=((aux68-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux70=(aux55-((v24**2)*v33))/((aux69-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux71=(-0.5*(v24**2))+((0.25*v44)+((0.5*(v11*v44))+((0.5*(v22*v44))+(v11*(v22*v44)))));
aux72=(-0.5*(v14**2))+((0.25*v22)+((0.5*(v11*v22))+((2.*(v12*(v14*v24)))+aux71)));
aux73=((0.125+((0.25*v11)+((-0.5*(v12**2))+aux72)))-((v12**2)*v44))-(v11*(v24**2));
aux74=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux75=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux74));
aux76=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux75)));
aux77=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux76)));
aux78=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux77)));
aux79=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux78));
aux80=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux79));
aux81=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux80)));
aux82=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux81)));
aux83=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux82));
aux84=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux83)));
aux85=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux84)));
aux86=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux85))))-((v12**2)*(v33*v44));
aux87=((aux86-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux88=(aux73-((v14**2)*v22))/((aux87-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux89=(v13*(v24**2))+((0.5*(v14*v34))+((v14*(v22*v34))+((-0.5*(v13*v44))+(v12*(v23*v44)))));
aux90=((-0.25*v13)+((-0.5*(v13*v22))+((0.5*(v12*v23))+aux89)))-(v13*(v22*v44));
aux91=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux92=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux91));
aux93=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux92)));
aux94=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux93)));
aux95=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux94)));
aux96=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux95));
aux97=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux96));
aux98=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux97)));
aux99=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux98)));
aux100=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux99));
aux101=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux100)));
aux102=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux101)));
aux103=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux102))))-((v12**2)*(v33*v44));
aux104=((aux103-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux105=((aux90-(v12*(v24*v34)))-(v14*(v23*v24)))/((aux104-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
output=(0.5*(((0.+-0.5j)*(aux17-aux35))-aux52))+(0.5*(((0.+0.5j)*(aux70-aux88))-aux105));
return output
def coeffscuadraticos14(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(-0.5*(b12*b33))+((b14*(b24*b33))+((b12*(b34**2))+((-0.5*(b12*b44))+(b13*(b23*b44)))));
aux1=(((-0.25*b12)+((0.5*(b13*b23))+((0.5*(b14*b24))+aux0)))-(b12*(b33*b44)))-(b13*(b24*b34));
aux2=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux3=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux2));
aux4=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux3)));
aux5=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux4)));
aux6=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux5)));
aux7=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux6));
aux8=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux7));
aux9=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux8)));
aux10=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux9)));
aux11=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux10));
aux12=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux11)));
aux13=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux12)));
aux14=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux13))))-((b12**2)*(b33*b44));
aux15=((aux14-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux16=(aux1-(b14*(b23*b34)))/((aux15-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux17=(b11*(b23*b24))+((-0.25*b34)+((-0.5*(b11*b34))+(((b12**2)*b34)+(-0.5*(b22*b34)))));
aux18=((0.5*(b13*b14))+((b13*(b14*b22))+((0.5*(b23*b24))+aux17)))-(b11*(b22*b34));
aux19=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux20=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux19));
aux21=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux20)));
aux22=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux21)));
aux23=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux22)));
aux24=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux23));
aux25=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux24));
aux26=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux25)));
aux27=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux26)));
aux28=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux27));
aux29=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux28)));
aux30=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux29)));
aux31=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux30))))-((b12**2)*(b33*b44));
aux32=((aux31-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux33=((aux18-(b12*(b13*b24)))-(b12*(b14*b23)))/((aux32-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux34=(0.5*(b12*b24))+((-0.5*(b14*b33))+((b12*(b24*b33))+((0.5*(b13*b34))+(b13*(b22*b34)))));
aux35=((-0.25*b14)+((-0.5*(b14*b22))+((b14*(b23**2))+aux34)))-(b12*(b23*b34));
aux36=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux37=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux36));
aux38=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux37)));
aux39=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux38)));
aux40=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux39)));
aux41=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux40));
aux42=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux41));
aux43=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux42)));
aux44=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux43)));
aux45=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux44));
aux46=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux45)));
aux47=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux46)));
aux48=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux47))))-((b12**2)*(b33*b44));
aux49=((aux48-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux50=((aux35-(b14*(b22*b33)))-(b13*(b23*b24)))/((aux49-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux51=(-0.5*(v12*v33))+((v14*(v24*v33))+((v12*(v34**2))+((-0.5*(v12*v44))+(v13*(v23*v44)))));
aux52=((-0.25*v12)+((0.5*(v13*v23))+((0.5*(v14*v24))+aux51)))-(v12*(v33*v44));
aux53=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux54=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux53));
aux55=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux54)));
aux56=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux55)));
aux57=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux56)));
aux58=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux57));
aux59=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux58));
aux60=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux59)));
aux61=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux60)));
aux62=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux61));
aux63=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux62)));
aux64=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux63)));
aux65=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux64))))-((v12**2)*(v33*v44));
aux66=((aux65-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux67=((aux52-(v13*(v24*v34)))-(v14*(v23*v34)))/((aux66-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux68=(v11*(v23*v24))+((-0.25*v34)+((-0.5*(v11*v34))+(((v12**2)*v34)+(-0.5*(v22*v34)))));
aux69=((0.5*(v13*v14))+((v13*(v14*v22))+((0.5*(v23*v24))+aux68)))-(v11*(v22*v34));
aux70=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux71=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux70));
aux72=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux71)));
aux73=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux72)));
aux74=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux73)));
aux75=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux74));
aux76=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux75));
aux77=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux76)));
aux78=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux77)));
aux79=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux78));
aux80=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux79)));
aux81=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux80)));
aux82=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux81))))-((v12**2)*(v33*v44));
aux83=((aux82-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux84=((aux69-(v12*(v13*v24)))-(v12*(v14*v23)))/((aux83-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux85=(0.5*(v12*v24))+((-0.5*(v14*v33))+((v12*(v24*v33))+((0.5*(v13*v34))+(v13*(v22*v34)))));
aux86=((-0.25*v14)+((-0.5*(v14*v22))+((v14*(v23**2))+aux85)))-(v12*(v23*v34));
aux87=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux88=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux87));
aux89=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux88)));
aux90=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux89)));
aux91=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux90)));
aux92=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux91));
aux93=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux92));
aux94=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux93)));
aux95=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux94)));
aux96=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux95));
aux97=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux96)));
aux98=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux97)));
aux99=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux98))))-((v12**2)*(v33*v44));
aux100=((aux99-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux101=((aux86-(v14*(v22*v33)))-(v13*(v23*v24)))/((aux100-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
output=(0.5*(((0.+-0.5j)*(aux16-aux33))-aux50))+(0.5*(((0.+0.5j)*(aux67-aux84))-aux101));
return output
def coeffscuadraticos22(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=((v13**2)*v24)+((v12*(v14*v33))+((-0.5*(v24*v33))+((0.5*(v23*v34))+(v11*(v23*v34)))));
aux1=((0.5*(v12*v14))+((-0.25*v24)+((-0.5*(v11*v24))+aux0)))-(v12*(v13*v34));
aux2=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux3=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux2));
aux4=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux3)));
aux5=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux4)));
aux6=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux5)));
aux7=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux6));
aux8=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux7));
aux9=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux8)));
aux10=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux9)));
aux11=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux10));
aux12=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux11)));
aux13=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux12)));
aux14=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux13))))-((v12**2)*(v33*v44));
aux15=((aux14-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux16=((0.+-1.j)*((aux1-(v11*(v24*v33)))-(v13*(v14*v23))))/((aux15-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux17=(-0.5*(v34**2))+((0.25*v44)+((0.5*(v11*v44))+((0.5*(v33*v44))+(v11*(v33*v44)))));
aux18=(-0.5*(v14**2))+((0.25*v33)+((0.5*(v11*v33))+((2.*(v13*(v14*v34)))+aux17)));
aux19=((0.125+((0.25*v11)+((-0.5*(v13**2))+aux18)))-((v13**2)*v44))-(v11*(v34**2));
aux20=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux21=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux20));
aux22=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux21)));
aux23=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux22)));
aux24=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux23)));
aux25=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux24));
aux26=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux25));
aux27=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux26)));
aux28=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux27)));
aux29=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux28));
aux30=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux29)));
aux31=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux30)));
aux32=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux31))))-((v12**2)*(v33*v44));
aux33=((aux32-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux34=(aux19-((v14**2)*v33))/((aux33-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux35=((b13**2)*b24)+((b12*(b14*b33))+((-0.5*(b24*b33))+((0.5*(b23*b34))+(b11*(b23*b34)))));
aux36=((0.5*(b12*b14))+((-0.25*b24)+((-0.5*(b11*b24))+aux35)))-(b12*(b13*b34));
aux37=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux38=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux37));
aux39=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux38)));
aux40=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux39)));
aux41=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux40)));
aux42=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux41));
aux43=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux42));
aux44=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux43)));
aux45=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux44)));
aux46=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux45));
aux47=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux46)));
aux48=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux47)));
aux49=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux48))))-((b12**2)*(b33*b44));
aux50=((aux49-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux51=((0.+1.j)*((aux36-(b11*(b24*b33)))-(b13*(b14*b23))))/((aux50-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux52=(-0.5*(b34**2))+((0.25*b44)+((0.5*(b11*b44))+((0.5*(b33*b44))+(b11*(b33*b44)))));
aux53=(-0.5*(b14**2))+((0.25*b33)+((0.5*(b11*b33))+((2.*(b13*(b14*b34)))+aux52)));
aux54=((0.125+((0.25*b11)+((-0.5*(b13**2))+aux53)))-((b13**2)*b44))-(b11*(b34**2));
aux55=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux56=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux55));
aux57=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux56)));
aux58=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux57)));
aux59=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux58)));
aux60=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux59));
aux61=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux60));
aux62=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux61)));
aux63=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux62)));
aux64=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux63));
aux65=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux64)));
aux66=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux65)));
aux67=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux66))))-((b12**2)*(b33*b44));
aux68=((aux67-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux69=(aux54-((b14**2)*b33))/((aux68-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
output=(0.25*(aux16-aux34))+(-0.25*(np.conjugate((aux51+aux69))));
return output
def coeffscuadraticos23(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(-0.5*(b12*b33))+((b14*(b24*b33))+((b12*(b34**2))+((-0.5*(b12*b44))+(b13*(b23*b44)))));
aux1=(((-0.25*b12)+((0.5*(b13*b23))+((0.5*(b14*b24))+aux0)))-(b12*(b33*b44)))-(b13*(b24*b34));
aux2=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux3=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux2));
aux4=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux3)));
aux5=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux4)));
aux6=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux5)));
aux7=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux6));
aux8=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux7));
aux9=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux8)));
aux10=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux9)));
aux11=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux10));
aux12=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux11)));
aux13=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux12)));
aux14=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux13))))-((b12**2)*(b33*b44));
aux15=((aux14-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux16=(aux1-(b14*(b23*b34)))/((aux15-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux17=(b11*(b23*b24))+((-0.25*b34)+((-0.5*(b11*b34))+(((b12**2)*b34)+(-0.5*(b22*b34)))));
aux18=((0.5*(b13*b14))+((b13*(b14*b22))+((0.5*(b23*b24))+aux17)))-(b11*(b22*b34));
aux19=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux20=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux19));
aux21=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux20)));
aux22=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux21)));
aux23=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux22)));
aux24=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux23));
aux25=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux24));
aux26=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux25)));
aux27=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux26)));
aux28=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux27));
aux29=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux28)));
aux30=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux29)));
aux31=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux30))))-((b12**2)*(b33*b44));
aux32=((aux31-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux33=((aux18-(b12*(b13*b24)))-(b12*(b14*b23)))/((aux32-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux34=((b14**2)*b23)+((0.5*(b24*b34))+((b11*(b24*b34))+((b12*(b13*b44))+(-0.5*(b23*b44)))));
aux35=((0.5*(b12*b13))+((-0.25*b23)+((-0.5*(b11*b23))+aux34)))-(b11*(b23*b44));
aux36=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux37=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux36));
aux38=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux37)));
aux39=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux38)));
aux40=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux39)));
aux41=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux40));
aux42=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux41));
aux43=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux42)));
aux44=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux43)));
aux45=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux44));
aux46=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux45)));
aux47=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux46)));
aux48=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux47))))-((b12**2)*(b33*b44));
aux49=((aux48-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux50=((aux35-(b12*(b14*b34)))-(b13*(b14*b24)))/((aux49-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux51=(-0.5*(v12*v33))+((v14*(v24*v33))+((v12*(v34**2))+((-0.5*(v12*v44))+(v13*(v23*v44)))));
aux52=((-0.25*v12)+((0.5*(v13*v23))+((0.5*(v14*v24))+aux51)))-(v12*(v33*v44));
aux53=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux54=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux53));
aux55=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux54)));
aux56=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux55)));
aux57=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux56)));
aux58=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux57));
aux59=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux58));
aux60=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux59)));
aux61=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux60)));
aux62=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux61));
aux63=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux62)));
aux64=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux63)));
aux65=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux64))))-((v12**2)*(v33*v44));
aux66=((aux65-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux67=((aux52-(v13*(v24*v34)))-(v14*(v23*v34)))/((aux66-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux68=(v11*(v23*v24))+((-0.25*v34)+((-0.5*(v11*v34))+(((v12**2)*v34)+(-0.5*(v22*v34)))));
aux69=((0.5*(v13*v14))+((v13*(v14*v22))+((0.5*(v23*v24))+aux68)))-(v11*(v22*v34));
aux70=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux71=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux70));
aux72=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux71)));
aux73=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux72)));
aux74=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux73)));
aux75=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux74));
aux76=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux75));
aux77=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux76)));
aux78=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux77)));
aux79=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux78));
aux80=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux79)));
aux81=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux80)));
aux82=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux81))))-((v12**2)*(v33*v44));
aux83=((aux82-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux84=((aux69-(v12*(v13*v24)))-(v12*(v14*v23)))/((aux83-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux85=((v14**2)*v23)+((0.5*(v24*v34))+((v11*(v24*v34))+((v12*(v13*v44))+(-0.5*(v23*v44)))));
aux86=((0.5*(v12*v13))+((-0.25*v23)+((-0.5*(v11*v23))+aux85)))-(v11*(v23*v44));
aux87=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux88=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux87));
aux89=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux88)));
aux90=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux89)));
aux91=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux90)));
aux92=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux91));
aux93=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux92));
aux94=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux93)));
aux95=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux94)));
aux96=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux95));
aux97=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux96)));
aux98=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux97)));
aux99=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux98))))-((v12**2)*(v33*v44));
aux100=((aux99-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux101=((aux86-(v12*(v14*v34)))-(v13*(v14*v24)))/((aux100-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
output=(0.5*(((0.+-0.5j)*(aux16-aux33))-aux50))+(0.5*(((0.+0.5j)*(aux67-aux84))-aux101));
return output
def coeffscuadraticos24(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(-0.5*(b34**2))+((0.25*b44)+((0.5*(b11*b44))+((0.5*(b33*b44))+(b11*(b33*b44)))));
aux1=(-0.5*(b14**2))+((0.25*b33)+((0.5*(b11*b33))+((2.*(b13*(b14*b34)))+aux0)));
aux2=((0.125+((0.25*b11)+((-0.5*(b13**2))+aux1)))-((b13**2)*b44))-(b11*(b34**2));
aux3=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux4=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux3));
aux5=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux4)));
aux6=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux5)));
aux7=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux6)));
aux8=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux7));
aux9=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux8));
aux10=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux9)));
aux11=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux10)));
aux12=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux11));
aux13=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux12)));
aux14=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux13)));
aux15=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux14))))-((b12**2)*(b33*b44));
aux16=((aux15-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux17=(aux2-((b14**2)*b33))/((aux16-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux18=(-0.5*(b23**2))+((0.25*b33)+((0.5*(b11*b33))+((0.5*(b22*b33))+(b11*(b22*b33)))));
aux19=(-0.5*(b13**2))+((0.25*b22)+((0.5*(b11*b22))+((2.*(b12*(b13*b23)))+aux18)));
aux20=((0.125+((0.25*b11)+((-0.5*(b12**2))+aux19)))-((b12**2)*b33))-(b11*(b23**2));
aux21=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux22=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux21));
aux23=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux22)));
aux24=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux23)));
aux25=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux24)));
aux26=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux25));
aux27=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux26));
aux28=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux27)));
aux29=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux28)));
aux30=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux29));
aux31=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux30)));
aux32=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux31)));
aux33=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux32))))-((b12**2)*(b33*b44));
aux34=((aux33-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux35=(aux20-((b13**2)*b22))/((aux34-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux36=((b13**2)*b24)+((b12*(b14*b33))+((-0.5*(b24*b33))+((0.5*(b23*b34))+(b11*(b23*b34)))));
aux37=((0.5*(b12*b14))+((-0.25*b24)+((-0.5*(b11*b24))+aux36)))-(b12*(b13*b34));
aux38=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux39=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux38));
aux40=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux39)));
aux41=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux40)));
aux42=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux41)));
aux43=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux42));
aux44=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux43));
aux45=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux44)));
aux46=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux45)));
aux47=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux46));
aux48=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux47)));
aux49=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux48)));
aux50=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux49))))-((b12**2)*(b33*b44));
aux51=((aux50-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux52=((aux37-(b11*(b24*b33)))-(b13*(b14*b23)))/((aux51-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux53=(-0.5*(v34**2))+((0.25*v44)+((0.5*(v11*v44))+((0.5*(v33*v44))+(v11*(v33*v44)))));
aux54=(-0.5*(v14**2))+((0.25*v33)+((0.5*(v11*v33))+((2.*(v13*(v14*v34)))+aux53)));
aux55=((0.125+((0.25*v11)+((-0.5*(v13**2))+aux54)))-((v13**2)*v44))-(v11*(v34**2));
aux56=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux57=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux56));
aux58=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux57)));
aux59=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux58)));
aux60=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux59)));
aux61=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux60));
aux62=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux61));
aux63=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux62)));
aux64=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux63)));
aux65=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux64));
aux66=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux65)));
aux67=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux66)));
aux68=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux67))))-((v12**2)*(v33*v44));
aux69=((aux68-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux70=(aux55-((v14**2)*v33))/((aux69-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux71=(-0.5*(v23**2))+((0.25*v33)+((0.5*(v11*v33))+((0.5*(v22*v33))+(v11*(v22*v33)))));
aux72=(-0.5*(v13**2))+((0.25*v22)+((0.5*(v11*v22))+((2.*(v12*(v13*v23)))+aux71)));
aux73=((0.125+((0.25*v11)+((-0.5*(v12**2))+aux72)))-((v12**2)*v33))-(v11*(v23**2));
aux74=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux75=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux74));
aux76=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux75)));
aux77=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux76)));
aux78=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux77)));
aux79=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux78));
aux80=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux79));
aux81=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux80)));
aux82=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux81)));
aux83=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux82));
aux84=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux83)));
aux85=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux84)));
aux86=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux85))))-((v12**2)*(v33*v44));
aux87=((aux86-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux88=(aux73-((v13**2)*v22))/((aux87-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux89=((v13**2)*v24)+((v12*(v14*v33))+((-0.5*(v24*v33))+((0.5*(v23*v34))+(v11*(v23*v34)))));
aux90=((0.5*(v12*v14))+((-0.25*v24)+((-0.5*(v11*v24))+aux89)))-(v12*(v13*v34));
aux91=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux92=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux91));
aux93=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux92)));
aux94=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux93)));
aux95=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux94)));
aux96=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux95));
aux97=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux96));
aux98=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux97)));
aux99=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux98)));
aux100=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux99));
aux101=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux100)));
aux102=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux101)));
aux103=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux102))))-((v12**2)*(v33*v44));
aux104=((aux103-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux105=((aux90-(v11*(v24*v33)))-(v13*(v14*v23)))/((aux104-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
output=(0.5*(((0.+-0.5j)*(aux17-aux35))-aux52))+(0.5*(((0.+0.5j)*(aux70-aux88))-aux105));
return output
def coeffscuadraticos33(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=(v13*(v24**2))+((0.5*(v14*v34))+((v14*(v22*v34))+((-0.5*(v13*v44))+(v12*(v23*v44)))));
aux1=((-0.25*v13)+((-0.5*(v13*v22))+((0.5*(v12*v23))+aux0)))-(v13*(v22*v44));
aux2=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux3=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux2));
aux4=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux3)));
aux5=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux4)));
aux6=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux5)));
aux7=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux6));
aux8=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux7));
aux9=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux8)));
aux10=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux9)));
aux11=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux10));
aux12=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux11)));
aux13=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux12)));
aux14=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux13))))-((v12**2)*(v33*v44));
aux15=((aux14-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux16=((0.+1.j)*((aux1-(v12*(v24*v34)))-(v14*(v23*v24))))/((aux15-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux17=(-0.5*(v24**2))+((0.25*v44)+((0.5*(v11*v44))+((0.5*(v22*v44))+(v11*(v22*v44)))));
aux18=(-0.5*(v14**2))+((0.25*v22)+((0.5*(v11*v22))+((2.*(v12*(v14*v24)))+aux17)));
aux19=((0.125+((0.25*v11)+((-0.5*(v12**2))+aux18)))-((v12**2)*v44))-(v11*(v24**2));
aux20=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux21=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux20));
aux22=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux21)));
aux23=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux22)));
aux24=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux23)));
aux25=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux24));
aux26=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux25));
aux27=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux26)));
aux28=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux27)));
aux29=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux28));
aux30=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux29)));
aux31=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux30)));
aux32=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux31))))-((v12**2)*(v33*v44));
aux33=((aux32-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux34=(aux19-((v14**2)*v22))/((aux33-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux35=(-0.5*(b24**2))+((0.25*b44)+((0.5*(b11*b44))+((0.5*(b22*b44))+(b11*(b22*b44)))));
aux36=(-0.5*(b14**2))+((0.25*b22)+((0.5*(b11*b22))+((2.*(b12*(b14*b24)))+aux35)));
aux37=((0.125+((0.25*b11)+((-0.5*(b12**2))+aux36)))-((b12**2)*b44))-(b11*(b24**2));
aux38=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux39=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux38));
aux40=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux39)));
aux41=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux40)));
aux42=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux41)));
aux43=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux42));
aux44=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux43));
aux45=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux44)));
aux46=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux45)));
aux47=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux46));
aux48=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux47)));
aux49=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux48)));
aux50=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux49))))-((b12**2)*(b33*b44));
aux51=((aux50-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux52=(aux37-((b14**2)*b22))/((aux51-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux53=(b13*(b24**2))+((0.5*(b14*b34))+((b14*(b22*b34))+((-0.5*(b13*b44))+(b12*(b23*b44)))));
aux54=((-0.25*b13)+((-0.5*(b13*b22))+((0.5*(b12*b23))+aux53)))-(b13*(b22*b44));
aux55=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux56=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux55));
aux57=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux56)));
aux58=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux57)));
aux59=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux58)));
aux60=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux59));
aux61=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux60));
aux62=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux61)));
aux63=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux62)));
aux64=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux63));
aux65=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux64)));
aux66=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux65)));
aux67=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux66))))-((b12**2)*(b33*b44));
aux68=((aux67-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux69=((0.+-1.j)*((aux54-(b12*(b24*b34)))-(b14*(b23*b24))))/((aux68-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
output=(0.25*(aux16-aux34))+(-0.25*(np.conjugate((aux52+aux69))));
return output
def coeffscuadraticos44(V,B):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
aux0=((v13**2)*v24)+((v12*(v14*v33))+((-0.5*(v24*v33))+((0.5*(v23*v34))+(v11*(v23*v34)))));
aux1=((0.5*(v12*v14))+((-0.25*v24)+((-0.5*(v11*v24))+aux0)))-(v12*(v13*v34));
aux2=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux3=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux2));
aux4=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux3)));
aux5=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux4)));
aux6=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux5)));
aux7=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux6));
aux8=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux7));
aux9=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux8)));
aux10=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux9)));
aux11=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux10));
aux12=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux11)));
aux13=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux12)));
aux14=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux13))))-((v12**2)*(v33*v44));
aux15=((aux14-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux16=((0.+1.j)*((aux1-(v11*(v24*v33)))-(v13*(v14*v23))))/((aux15-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux17=(-0.5*(v23**2))+((0.25*v33)+((0.5*(v11*v33))+((0.5*(v22*v33))+(v11*(v22*v33)))));
aux18=(-0.5*(v13**2))+((0.25*v22)+((0.5*(v11*v22))+((2.*(v12*(v13*v23)))+aux17)));
aux19=((0.125+((0.25*v11)+((-0.5*(v12**2))+aux18)))-((v12**2)*v33))-(v11*(v23**2));
aux20=(0.25*(v33*v44))+((0.5*(v11*(v33*v44)))+((0.5*(v22*(v33*v44)))+(v11*(v22*(v33*v44)))));
aux21=(0.5*(v11*(v22*v44)))+((2.*(v12*(v13*(v23*v44))))+((-0.5*((v23**2)*v44))+aux20));
aux22=(0.25*(v11*v44))+((-0.5*((v12**2)*v44))+((-0.5*((v13**2)*v44))+((0.25*(v22*v44))+aux21)));
aux23=(-0.5*(v11*(v34**2)))+(((v12**2)*(v34**2))+((-0.5*(v22*(v34**2)))+((0.125*v44)+aux22)));
aux24=(-2.*(v12*(v13*(v24*v34))))+((v23*(v24*v34))+((2.*(v11*(v23*(v24*v34))))+((-0.25*(v34**2))+aux23)));
aux25=(v13*(v14*v34))+((2.*(v13*(v14*(v22*v34))))+((-2.*(v12*(v14*(v23*v34))))+aux24));
aux26=(0.5*(v11*(v22*v33)))+((2.*(v12*(v14*(v24*v33))))+((-0.5*((v24**2)*v33))+aux25));
aux27=(0.25*(v11*v33))+((-0.5*((v12**2)*v33))+((-0.5*((v14**2)*v33))+((0.25*(v22*v33))+aux26)));
aux28=(-0.25*(v24**2))+((-0.5*(v11*(v24**2)))+(((v13**2)*(v24**2))+((0.125*v33)+aux27)));
aux29=((v14**2)*(v23**2))+((v12*(v14*v24))+((-2.*(v13*(v14*(v23*v24))))+aux28));
aux30=(-0.5*((v14**2)*v22))+((v12*(v13*v23))+((-0.25*(v23**2))+((-0.5*(v11*(v23**2)))+aux29)));
aux31=(-0.25*(v14**2))+((0.125*v22)+((0.25*(v11*v22))+((-0.5*((v13**2)*v22))+aux30)));
aux32=(0.0625+((0.125*v11)+((-0.25*(v12**2))+((-0.25*(v13**2))+aux31))))-((v12**2)*(v33*v44));
aux33=((aux32-(v11*((v23**2)*v44)))-((v13**2)*(v22*v44)))-(v11*(v22*(v34**2)));
aux34=(aux19-((v13**2)*v22))/((aux33-(v11*((v24**2)*v33)))-((v14**2)*(v22*v33)));
aux35=(-0.5*(b23**2))+((0.25*b33)+((0.5*(b11*b33))+((0.5*(b22*b33))+(b11*(b22*b33)))));
aux36=(-0.5*(b13**2))+((0.25*b22)+((0.5*(b11*b22))+((2.*(b12*(b13*b23)))+aux35)));
aux37=((0.125+((0.25*b11)+((-0.5*(b12**2))+aux36)))-((b12**2)*b33))-(b11*(b23**2));
aux38=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux39=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux38));
aux40=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux39)));
aux41=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux40)));
aux42=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux41)));
aux43=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux42));
aux44=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux43));
aux45=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux44)));
aux46=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux45)));
aux47=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux46));
aux48=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux47)));
aux49=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux48)));
aux50=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux49))))-((b12**2)*(b33*b44));
aux51=((aux50-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux52=(aux37-((b13**2)*b22))/((aux51-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
aux53=((b13**2)*b24)+((b12*(b14*b33))+((-0.5*(b24*b33))+((0.5*(b23*b34))+(b11*(b23*b34)))));
aux54=((0.5*(b12*b14))+((-0.25*b24)+((-0.5*(b11*b24))+aux53)))-(b12*(b13*b34));
aux55=(0.25*(b33*b44))+((0.5*(b11*(b33*b44)))+((0.5*(b22*(b33*b44)))+(b11*(b22*(b33*b44)))));
aux56=(0.5*(b11*(b22*b44)))+((2.*(b12*(b13*(b23*b44))))+((-0.5*((b23**2)*b44))+aux55));
aux57=(0.25*(b11*b44))+((-0.5*((b12**2)*b44))+((-0.5*((b13**2)*b44))+((0.25*(b22*b44))+aux56)));
aux58=(-0.5*(b11*(b34**2)))+(((b12**2)*(b34**2))+((-0.5*(b22*(b34**2)))+((0.125*b44)+aux57)));
aux59=(-2.*(b12*(b13*(b24*b34))))+((b23*(b24*b34))+((2.*(b11*(b23*(b24*b34))))+((-0.25*(b34**2))+aux58)));
aux60=(b13*(b14*b34))+((2.*(b13*(b14*(b22*b34))))+((-2.*(b12*(b14*(b23*b34))))+aux59));
aux61=(0.5*(b11*(b22*b33)))+((2.*(b12*(b14*(b24*b33))))+((-0.5*((b24**2)*b33))+aux60));
aux62=(0.25*(b11*b33))+((-0.5*((b12**2)*b33))+((-0.5*((b14**2)*b33))+((0.25*(b22*b33))+aux61)));
aux63=(-0.25*(b24**2))+((-0.5*(b11*(b24**2)))+(((b13**2)*(b24**2))+((0.125*b33)+aux62)));
aux64=((b14**2)*(b23**2))+((b12*(b14*b24))+((-2.*(b13*(b14*(b23*b24))))+aux63));
aux65=(-0.5*((b14**2)*b22))+((b12*(b13*b23))+((-0.25*(b23**2))+((-0.5*(b11*(b23**2)))+aux64)));
aux66=(-0.25*(b14**2))+((0.125*b22)+((0.25*(b11*b22))+((-0.5*((b13**2)*b22))+aux65)));
aux67=(0.0625+((0.125*b11)+((-0.25*(b12**2))+((-0.25*(b13**2))+aux66))))-((b12**2)*(b33*b44));
aux68=((aux67-(b11*((b23**2)*b44)))-((b13**2)*(b22*b44)))-(b11*(b22*(b34**2)));
aux69=((0.+-1.j)*((aux54-(b11*(b24*b33)))-(b13*(b14*b23))))/((aux68-(b11*((b24**2)*b33)))-((b14**2)*(b22*b33)));
output=(0.25*(aux16-aux34))+(-0.25*(np.conjugate((aux52+aux69))));
return output
def lineal1(V,B,R,S):
v11,v12,v13,v14=separate(V,0)
v21,v22,v23,v24=separate(V,1)
v31,v32,v33,v34=separate(V,2)
v41,v42,v43,v44=separate(V,3)
b11,b12,b13,b14=separate(B,0)
b21,b22,b23,b24=separate(B,1)
b31,b32,b33,b34=separate(B,2)
b41,b42,b43,b44=separate(B,3)
r1,r2,r3,r4=R[0],R[1],R[2],R[3]
s1,s2,s3,s4=S[0],S[1],S[2],S[3]
aux0=(2.*b33)+(((0.+-8.j)*(b24*b33))+((b22*(2.+(4.*b33)))+((0.+8.j)*(b23*b34))));
aux1=((0.+-2.j)*(b14*b34))+(((0.+-2.j)*(b23*b34))+((-4.*(b34**2))+((2.*b44)+(4.*(b33*b44)))));
aux2=b13+(b24+((2.*(b24*b33))+((2.*(b13*b44))+((0.+-1.j)*(1.+((2.*b33)+aux1))))));
aux3=(2.*(b23*b44))+((4.*(b11*(b23*b44)))+(((0.+-8.j)*(b13*(b23*b44)))+(-2.*(b12*aux2))));
aux4=(b14*(1.+((-4.*(b23**2))+(((0.+-4.j)*b24)+((4.*(b13*b24))+aux0)))))+aux3;
aux5=(-2.*(b24*b34))+((-4.*(b11*(b24*b34)))+(((0.+8.j)*(b13*(b24*b34)))+aux4));
aux6=(-4.*((b14**2)*b23))+((4.*(b13*(b23*b24)))+((-2.*(b13*b34))+((-4.*(b13*(b22*b34)))+aux5)));
aux7=(b14*(b24+((2.*(b24*b33))+(-2.*(b23*b34)))))+(b13*(b23+((-2.*(b24*b34))+(2.*(b23*b44)))));
aux8=(-4.*((b23**2)*(1.+(2.*b44))))+(b22*(2.+((-8.*(b34**2))+((4.*b44)+(b33*(4.+(8.*b44)))))));
aux9=(16.*(b23*(b24*b34)))+((-4.*(b34**2))+((2.*b44)+((4.*(b33*b44))+aux8)));
aux10=(16.*(b12*aux7))+(2.*(b11*(1.+((-4.*(b24**2))+((2.*b33)+((-8.*((b24**2)*b33))+aux9))))));
aux11=(-4.*((b12**2)*(1.+((-4.*(b34**2))+((2.*b44)+(b33*(2.+(4.*b44))))))))+aux10;
aux12=(-16.*((b13**2)*(b22*b44)))+((-8.*((b23**2)*b44))+((4.*(b33*b44))+((8.*(b22*(b33*b44)))+aux11)));
aux13=(-8.*(b22*(b34**2)))+((2.*b44)+((-8.*((b13**2)*b44))+((4.*(b22*b44))+aux12)));
aux14=(32.*(b13*(b14*(b22*b34))))+((16.*(b23*(b24*b34)))+((-4.*(b34**2))+aux13));
aux15=(-16.*((b14**2)*(b22*b33)))+((-8.*((b24**2)*b33))+((16.*(b13*(b14*b34)))+aux14));
aux16=(16.*((b13**2)*(b24**2)))+((2.*b33)+((-8.*((b14**2)*b33))+((4.*(b22*b33))+aux15)));
aux17=(16.*((b14**2)*(b23**2)))+((-32.*(b13*(b14*(b23*b24))))+((-4.*(b24**2))+aux16));
aux18=(2.*b22)+((-8.*((b13**2)*b22))+((-8.*((b14**2)*b22))+((-4.*(b23**2))+aux17)));
aux19=((0.+2.j)*((b23+((2.*(b11*b23))+(((0.+-4.j)*(b13*b23))+aux6)))*s2))/(1.+((-4.*(b13**2))+((-4.*(b14**2))+aux18)));
aux20=(b14*((-2.*(b23*b24))+(b34+(2.*(b22*b34)))))+(b12*(b23+((-2.*(b24*b34))+(2.*(b23*b44)))));
aux21=(-4.*((b23**2)*b44))+((2.*(b33*b44))+((4.*(b22*(b33*b44)))+((b12**2)*(2.+(4.*b44)))));
aux22=(-4.*((b24**2)*b33))+((8.*(b23*(b24*b34)))+((-2.*(b34**2))+((-4.*(b22*(b34**2)))+aux21)));
aux23=(4.*((b14**2)*b22))+((-2.*(b23**2))+((-8.*(b12*(b14*b24)))+(b33+((2.*(b22*b33))+aux22))));
aux24=((2.*(b14**2))+aux23)-(b11*(1.+((-4.*(b24**2))+((2.*b44)+(b22*(2.+(4.*b44)))))));
aux25=(b14*(b24+((2.*(b24*b33))+(-2.*(b23*b34)))))+(b13*(b23+((-2.*(b24*b34))+(2.*(b23*b44)))));
aux26=(-4.*((b23**2)*(1.+(2.*b44))))+(b22*(2.+((-8.*(b34**2))+((4.*b44)+(b33*(4.+(8.*b44)))))));
aux27=(16.*(b23*(b24*b34)))+((-4.*(b34**2))+((2.*b44)+((4.*(b33*b44))+aux26)));
aux28=-2.*(b11*(1.+((-4.*(b24**2))+((2.*b33)+((-8.*((b24**2)*b33))+aux27)))));
aux29=(4.*((b12**2)*(1.+((-4.*(b34**2))+((2.*b44)+(b33*(2.+(4.*b44))))))))+((-16.*(b12*aux25))+aux28);
aux30=(16.*((b13**2)*(b22*b44)))+((8.*((b23**2)*b44))+((-4.*(b33*b44))+((-8.*(b22*(b33*b44)))+aux29)));
aux31=(8.*(b22*(b34**2)))+((-2.*b44)+((8.*((b13**2)*b44))+((-4.*(b22*b44))+aux30)));
aux32=(-32.*(b13*(b14*(b22*b34))))+((-16.*(b23*(b24*b34)))+((4.*(b34**2))+aux31));
aux33=(16.*((b14**2)*(b22*b33)))+((8.*((b24**2)*b33))+((-16.*(b13*(b14*b34)))+aux32));
aux34=(-16.*((b13**2)*(b24**2)))+((-2.*b33)+((8.*((b14**2)*b33))+((-4.*(b22*b33))+aux33)));
aux35=(-16.*((b14**2)*(b23**2)))+((32.*(b13*(b14*(b23*b24))))+((4.*(b24**2))+aux34));
aux36=(-2.*b22)+((8.*((b13**2)*b22))+((8.*((b14**2)*b22))+((4.*(b23**2))+aux35)));
aux37=((0.+2.j)*aux24)+((0.+1.j)*(-1.+((4.*(b13**2))+((4.*(b14**2))+aux36))));
aux38=(-4.*(b13*(1.+((-4.*(b24**2))+((2.*b44)+(b22*(2.+(4.*b44))))))))+((8.*aux20)+aux37);
aux39=(b14*(b24+((2.*(b24*b33))+(-2.*(b23*b34)))))+(b13*(b23+((-2.*(b24*b34))+(2.*(b23*b44)))));
aux40=(-4.*((b23**2)*(1.+(2.*b44))))+(b22*(2.+((-8.*(b34**2))+((4.*b44)+(b33*(4.+(8.*b44)))))));
aux41=(16.*(b23*(b24*b34)))+((-4.*(b34**2))+((2.*b44)+((4.*(b33*b44))+aux40)));
aux42=(16.*(b12*aux39))+(2.*(b11*(1.+((-4.*(b24**2))+((2.*b33)+((-8.*((b24**2)*b33))+aux41))))));
aux43=(-4.*((b12**2)*(1.+((-4.*(b34**2))+((2.*b44)+(b33*(2.+(4.*b44))))))))+aux42;
aux44=(-16.*((b13**2)*(b22*b44)))+((-8.*((b23**2)*b44))+((4.*(b33*b44))+((8.*(b22*(b33*b44)))+aux43)));
aux45=(-8.*(b22*(b34**2)))+((2.*b44)+((-8.*((b13**2)*b44))+((4.*(b22*b44))+aux44)));
aux46=(32.*(b13*(b14*(b22*b34))))+((16.*(b23*(b24*b34)))+((-4.*(b34**2))+aux45));
aux47=(-16.*((b14**2)*(b22*b33)))+((-8.*((b24**2)*b33))+((16.*(b13*(b14*b34)))+aux46));
aux48=(16.*((b13**2)*(b24**2)))+((2.*b33)+((-8.*((b14**2)*b33))+((4.*(b22*b33))+aux47)));
aux49=(16.*((b14**2)*(b23**2)))+((-32.*(b13*(b14*(b23*b24))))+((-4.*(b24**2))+aux48));
aux50=(2.*b22)+((-8.*((b13**2)*b22))+((-8.*((b14**2)*b22))+((-4.*(b23**2))+aux49)));
aux51=(b13*((-2.*(b23*b24))+(b34+(2.*(b22*b34)))))+(b12*(b24+((2.*(b24*b33))+(-2.*(b23*b34)))));
aux52=-2.*(b13*((b14+((2.*(b14*b22))+((2.*(b24*b34))+(-2.*(b23*b44)))))-b23));
aux53=(4.*(b14*b23))+((4.*(b13*b24))+((-2.*b33)+((4.*(b34**2))+((-2.*b44)+(-4.*(b33*b44))))));
aux54=(2.*(b22*b34))+((4.*(b11*(b22*b34)))+((-4.*(b14*(b23*b34)))+(aux52+(b12*(-1.+aux53)))));
aux55=(4.*(b14*(b24*b33)))+(b34+((2.*(b11*b34))+((-4.*((b12**2)*b34))+aux54)));
aux56=(0.+2.j)*((2.*(b14*b24))+((-2.*(b23*b24))+((-4.*(b11*(b23*b24)))+aux55)));
aux57=(-4.*(b14*(1.+((-4.*(b23**2))+((2.*b33)+(b22*(2.+(4.*b33))))))))+((8.*aux51)+aux56);
aux58=(b14*(b24+((2.*(b24*b33))+(-2.*(b23*b34)))))+(b13*(b23+((-2.*(b24*b34))+(2.*(b23*b44)))));
aux59=(-4.*((b23**2)*(1.+(2.*b44))))+(b22*(2.+((-8.*(b34**2))+((4.*b44)+(b33*(4.+(8.*b44)))))));
aux60=(16.*(b23*(b24*b34)))+((-4.*(b34**2))+((2.*b44)+((4.*(b33*b44))+aux59)));
aux61=(16.*(b12*aux58))+(2.*(b11*(1.+((-4.*(b24**2))+((2.*b33)+((-8.*((b24**2)*b33))+aux60))))));
aux62=(-4.*((b12**2)*(1.+((-4.*(b34**2))+((2.*b44)+(b33*(2.+(4.*b44))))))))+aux61;
aux63=(-16.*((b13**2)*(b22*b44)))+((-8.*((b23**2)*b44))+((4.*(b33*b44))+((8.*(b22*(b33*b44)))+aux62)));
aux64=(-8.*(b22*(b34**2)))+((2.*b44)+((-8.*((b13**2)*b44))+((4.*(b22*b44))+aux63)));
aux65=(32.*(b13*(b14*(b22*b34))))+((16.*(b23*(b24*b34)))+((-4.*(b34**2))+aux64));
aux66=(-16.*((b14**2)*(b22*b33)))+((-8.*((b24**2)*b33))+((16.*(b13*(b14*b34)))+aux65));
aux67=(16.*((b13**2)*(b24**2)))+((2.*b33)+((-8.*((b14**2)*b33))+((4.*(b22*b33))+aux66)));
aux68=(16.*((b14**2)*(b23**2)))+((-32.*(b13*(b14*(b23*b24))))+((-4.*(b24**2))+aux67));
aux69=(2.*b22)+((-8.*((b13**2)*b22))+((-8.*((b14**2)*b22))+((-4.*(b23**2))+aux68)));