Computes the pdf, cdf, value at risk and expected shortfall for the power function I distribution given by
$$\begin{array}{ll}
&\displaystyle
f (x) = a x^{a - 1},
\\
&\displaystyle
F (x) = x^a,
\\
&\displaystyle
{\rm VaR}_p (X) = p^{1 / a},
\\
&\displaystyle
{\rm ES}_p (X) = \frac {p^{1 / a}}{1 / a + 1}
\end{array}$$
for \(0 < x < 1\), \(0 < p < 1\), and \(a > 0\), the shape parameter.