function [ X ] = rpareto( n,a )

U=rand(1,n);
X=(1-U)^(-1/a);

endfunction

