Let x be a random variable with probability density function (pdf) x ~ Gamma(1, 1). How to plot the pdf of z where z=x^2/(1+sin(x))? This is what I have tried, I am not sure how to get z from x.
x = np.linspace (0, 100, 200)
y1 = stats.gamma.pdf(x, a=1, loc=1)
plt.plot(x, y1, "y-", label=(r'$\alpha=29, \beta=3$'))