Clear["Global ‘*"]​​SetDirectory[NotebookDirectory[]];​​f[x_]=Module{coeffi,clist,inner,inlist,Integral0,Integral1,Integral2,Integral3,Integral4,Integral5,Integral6,Ke,KeN},​​coeffi=CoefficientExpand
3
x
1
3
-
x
1
2
+
x
2
3
x
1
2
+
x
2
3
-
2
x
2
3
+
x
3
3
-
x
1
2
+
x
2
3
+
x
3
3
x
1
2
+
x
2
3
+
x
3
3
-
3
x
3
4
+
x
4
3
-
2
x
2
3
+
x
3
4
+
x
4
3
-
x
1
2
+
x
2
3
+
x
3
4
+
x
4
3
x
1
2
+
x
2
3
+
x
3
4
+
x
4
,
x
1
,3;​​clist=coeffi/.PlusList;​​inner=Inner[Times,clist/.{
x
2
->1,
x
3
->1},clist,Plus];​​inlist=inner/.Plus->List/.
x
2
->
2
r
2
,
x
3
->
2
r
3
;​​(*fost=Times[inlist,
r
2
r
3
r
4
]*)​​Integral0=Integrateinlist*
r
3
*Exp-
3
2
r
3
8
2
l
0
,{
r
3
,0,Infinity},Assumptions(
l
0
>0);​​Integral1=IntegrateIntegral0*
r
2
*Exp-
2
r
2
3
2
l
0
,{
r
2
,0,Infinity},Assumptions(
l
0
>0);​​Integral2=IntegrateIntegral1*
r
4
*Exp-
2
2
r
4
5
2
l
0
,{
r
4
,0,Infinity},Assumptions(
l
0
>0);​​Integral3=IntegrateIntegral2*
r
5
*Exp-
5
2
r
5
2
2
l
0
,{
r
5
,0,Infinity},Assumptions(
l
0
>0);​​Integral4=Plus@@Integral3;​​Integral5=Integrate
6
r
1
*Exp-
2
r
1
4
2
l
0
,{
r
1
,0,Infinity},Assumptions(
l
0
>0);​​Integral6=Integrate
7
r
1
*Exp-
2
r
1
4
2
l
0
,{
r
1
,0,Infinity},Assumptions(
l
0
>0);​​Ke=Times[Integral4,Integral5];​​KeN=Times[Integral4,Integral6];​​Export["sortie0.nb",Ke,"List"];​​Export["sortie0N.nb",KeN,"List"];​​Print[Ke];​​Print[KeN];
120
π
7
l
0
1563182913960349320038400000
62
l
0
+1578062819114272486195200000
60
l
0
2
x
4
+442318249100057121792000000
58
l
0
4
x
4
+70155086383161645465600000
56
l
0
6
x
4
+6595763863461128294400000
54
l
0
8
x
4
+751955283352842928128000
52
l
0
10
x
4
+95247002909672570880000
50
l
0
12
x
4
+5387410790582059008000
48
l
0
14
x
4
+384733907878957056000
46
l
0
16
x
4
+15775328416038912000
44
l
0
18
x
4
+1302883989700608000
42
l
0
20
x
4
+1688579923968000
38
l
0
24
x
4

768
8
l
0
1563182913960349320038400000
62
l
0
+1578062819114272486195200000
60
l
0
2
x
4
+442318249100057121792000000
58
l
0
4
x
4
+70155086383161645465600000
56
l
0
6
x
4
+6595763863461128294400000
54
l
0
8
x
4
+751955283352842928128000
52
l
0
10
x
4
+95247002909672570880000
50
l
0
12
x
4
+5387410790582059008000
48
l
0
14
x
4
+384733907878957056000
46
l
0
16
x
4
+15775328416038912000
44
l
0
18
x
4
+1302883989700608000
42
l
0
20
x
4
+1688579923968000
38
l
0
24
x
4

1013147954380800000
π
43
l
0
1013147954380800000
π
43
l
0