function critical
clf;
v=0.5;
r=linspace(3.55,4,100);
hold on
for i=1:10
    v=myfn(r,v);
    plot(r,v);
end

hold off

function myval=myfn(r,x)
myval=r.*x.*(1-x);