clear all;
r = [2.8:0.005:4.0];
for j = 1:length(r);
  y = iterate(r(j));
  ypts(:,j) = y;
end
plot(r,ypts,'k.'); 
axis([min(r) max(r) 0 1])