% 1.731 Monte Carlo approximation to derived distribution
function derived_dist1
% Program for deriving distribution of y=x^2 when y is uniformly distributed
nrep=1000
% begin replicate loop
for i=1:nrep
    x(i)=unifrnd(-1,1);
    y(i)=x(i)^2;
end
close all
xmean=mean(x);
ymean=mean(y);
stdx=std(x)
stdy=std(y)
figure
hist(x,20)
figure
hist(y,20)
figure
cdfplot(x)
figure
cdfplot(y)
return