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