我正在尝试绘制f(x)= r x (1-x)的蜘蛛网图,其中r = 3.我目前正在尝试使用以下代码:
function cobweb(f,a,b,x0,x1,N)
x(1)=0.2; % plot orbit starting at x0
for i=1:100
x(i+1)= 3*x(i)*(x(i)+1);
plot([x(i),x(i)],[x(i),x(i+1)]);
hold on
plot([x(i),x(i+1)],[x(i+1),x(i+1)]);
hold on
end
hold off
但我得到的情节完全错了!