-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathexp_well_def.hh
1296 lines (1296 loc) · 86.9 KB
/
exp_well_def.hh
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
gsl_matrix_complex_set(output, 0, 0, gsl_complex_rect(((2-2*pow(M_E,(a1+a2)*ra)+2*(a1+a2)*ra+pow(a1+a2,2)*pow(ra,2))*V)/(pow(a1+a2,3)*pow(M_E,(a1+a2)*ra)), 0));
gsl_matrix_complex_set(output, 0, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 0, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 1, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 1, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 1, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 2, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 2, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 2, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 3, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-(a1+a2)*ra*(720+(a1+a2)*ra*(360+(a1+a2)*ra*(120+(a1+a2)*ra*(30+(a1+a2)*ra*(6+(a1+a2)*ra))))))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 3, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 3, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 4, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 4, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 4, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 5, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 5, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 5, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 6, gsl_complex_rect(((2-2*pow(M_E,(a1+a2)*ra)+2*(a1+a2)*ra+pow(a1+a2,2)*pow(ra,2))*V)/(pow(a1+a2,3)*pow(M_E,(a1+a2)*ra)), 0));
gsl_matrix_complex_set(output, 6, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 6, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 7, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 7, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 7, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 8, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 8, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 8, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 9, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-(a1+a2)*ra*(720+(a1+a2)*ra*(360+(a1+a2)*ra*(120+(a1+a2)*ra*(30+(a1+a2)*ra*(6+(a1+a2)*ra))))))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 9, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 9, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 10, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 10, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 10, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 11, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 11, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 11, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 12, gsl_complex_rect(((2-2*pow(M_E,(a1+a2)*ra)+2*(a1+a2)*ra+pow(a1+a2,2)*pow(ra,2))*V)/(pow(a1+a2,3)*pow(M_E,(a1+a2)*ra)), 0));
gsl_matrix_complex_set(output, 12, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 12, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 13, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 13, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 13, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 14, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 14, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 14, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 15, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-(a1+a2)*ra*(720+(a1+a2)*ra*(360+(a1+a2)*ra*(120+(a1+a2)*ra*(30+(a1+a2)*ra*(6+(a1+a2)*ra))))))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 15, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 15, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 16, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 16, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 16, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 17, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 17, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 17, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 18, gsl_complex_rect(((2-2*pow(M_E,(a1+a2)*ra)+2*(a1+a2)*ra+pow(a1+a2,2)*pow(ra,2))*V)/(pow(a1+a2,3)*pow(M_E,(a1+a2)*ra)), 0));
gsl_matrix_complex_set(output, 18, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 18, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 19, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 19, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 19, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 20, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 20, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 20, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 21, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-(a1+a2)*ra*(720+(a1+a2)*ra*(360+(a1+a2)*ra*(120+(a1+a2)*ra*(30+(a1+a2)*ra*(6+(a1+a2)*ra))))))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 21, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 21, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 22, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 22, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 22, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 23, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 23, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 23, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 24, gsl_complex_rect(((2-2*pow(M_E,(a1+a2)*ra)+2*(a1+a2)*ra+pow(a1+a2,2)*pow(ra,2))*V)/(pow(a1+a2,3)*pow(M_E,(a1+a2)*ra)), 0));
gsl_matrix_complex_set(output, 24, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 24, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 25, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 25, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 25, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 26, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-720*(a1+a2)*ra-360*pow(a1+a2,2)*pow(ra,2)-120*pow(a1+a2,3)*pow(ra,3)-30*pow(a1+a2,4)*pow(ra,4)-6*pow(a1+a2,5)*pow(ra,5)-pow(a1+a2,6)*pow(ra,6))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));
gsl_matrix_complex_set(output, 26, 27, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 28, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 29, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 30, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 31, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 32, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 33, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 34, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 26, 35, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 0, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 1, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 2, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 3, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 4, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 5, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 6, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 7, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 8, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 9, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 10, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 11, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 12, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 13, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 14, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 15, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 16, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 17, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 18, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 19, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 20, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 21, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 22, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 23, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 24, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 25, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 26, gsl_complex_rect(0, 0));
gsl_matrix_complex_set(output, 27, 27, gsl_complex_rect(-(((-720+720*pow(M_E,(a1+a2)*ra)-(a1+a2)*ra*(720+(a1+a2)*ra*(360+(a1+a2)*ra*(120+(a1+a2)*ra*(30+(a1+a2)*ra*(6+(a1+a2)*ra))))))*V)/(pow(a1+a2,7)*pow(M_E,(a1+a2)*ra))), 0));