pbr=0.01;
a=1:1000;
fa=exp((a-1).*log(1-pbr));
plot(log(a),log(fa));