function [ X ] = rpareto( n,a )
%UNTITLED5 Summary of this function goes here
%   Detailed explanation goes here
U=rand(n,1);
X=(1-U).^(-1/a);

end

