In[]:=
f[t_]=2
-
3
160
(-32+40
2
t
-20
3
t
+
5
t
)
0≤t<2
0
True
;​​distF=ProbabilityDistribution[f[t],{t,0,2}];​​distG=TransformedDistribution[x+y,{x,y}ProductDistribution[{distF,2}]];​​trunc=TruncatedDistribution0,
4
Sqrt[3]
,distG;​​PDF[trunc,t]​​​​​​
Out[]=
-
2187
7
(-4+t)
(80+536t+228
2
t
+28
3
t
+
4
t
)
1024(88347959-48611328
3
)
2≤t≤
4
3
4209975
36t
25
-
6
3
t
5
+
9
4
t
20
+
3
5
t
40
-
21
6
t
400
+
9
7
t
2240
+
3
8
t
4480
-
9
t
8960
+
11
t
1971200
88347959-48611328
3
0<t<2
0
True