%xpowery.m/created by PJNahin for "Duelling Idiots"(1/6/99)
%This m-file plots the probability density function of X^Y,
%where X and Y are independent, and both uniform from 0 to 1.
%
%
z=.0001:.01:1;
Z=expint(-log(z))./z;
plot(z,Z)
grid
title('Fig.13.2 - PDF of Z = X^Y, with X and Y uniform from 0 to 1')
ylabel('Probability Density')
xlabel('z')
figure(1)