I'm trying to find the result of the sum of two PDF that is the square of two continuous uniform distribution (rectangular distribution) that have $\mu =0$ and $\sigma =1$. The square of this continuous uniform distribution is defined from 0 to 3 and is $$p(x)= \frac{1}{2 \sqrt{3x}}$$ Now i need to find the PDF of the sum of two random variables with p(x) as distribution. The convolution integral is $$p(u)=\int_0^3 p(u-t)p(t)dt=\int_0^3 \frac{dt}{12 \sqrt{t(u-t)}}= \frac{arcsin \left( \sqrt{ \frac{3}{u}} \right)}{6}$$ This funcion is not defined on R from 0 to 6. If i extract the real part of $(p(u))$ i obtain the correct PDF from 0 to 3 but the part from 3 to 6 is wrong and p(u) is not normalized from 0 to 6. I simulate this PDF with python, here is an image
Simulation of PDF: https://imgur.com/ImqT7ya
Plot of the PDF: https://imgur.com/khg8rGM
what did I do wrong?