Fiz esta rotina para o domínio 1,infinito.
Porém como eu faço para a curva da densidade começar no 1 assim como no
histograma?
d2=runif(100)^-1/alfa
d3=subset(d2,d2>=1)
hist(d3,probability=T,nclass=60,col='green')
dens=density(d3,n=200,na.rm=T)
points(dens,type='l',col=2)
abline(v=1,lty=4)