A= 0.1529; B=0.1933; prop(1) = 0.01; bdist(1) = betapdf(prop(1),A,B); for i = 2:100 prop(i) = prop(i-1)+.01; bdist(i) = betapdf(prop(i),A,B); end plot(prop,bdist);