random number with p(x)= x^(-a) distribution

How can I generate random number within [1,inf] with below distribution in MATLAB:

p(x)= x^(-a)