%program to generate Pareto random numbers
%according to the inversion method.
points=1000; %generate 1000 points
a = 2;
b = 2.5;

x=rand(points,1);
for i=1:points,
   y1=1-x(i);
   y2=-(1/b)*(log(y1));
   y3(i)=a*exp(y2);
end
hold on 
subplot(1,2,1)
plot(y3,'k');
box off
axis square 
xlabel('Sample index')
ylabel('Random number value')
%generate a histogram of the data to view the pdf
x = 0:0.1:1;
subplot(1,2,2)
hist(y3)
box off
axis square 
xlabel('Bins')
ylabel('Bin count')

