pbr=0.01;
a=1:1000;
fa=(1-pbr).^(1-a);
plot(a,fa);