混沌动力学分岔图，代码有误，求教~

antonio1111

function s=der2(t,y,Z,flag,f);
ylag1=Z(:,1);
s=[493.8*y(2)*(y(1)+1)-1.2*y(1)*(y(1)+1);
1.458*4*f/3*(1+ylag1(1))-1.458*y(2)-0.486*y(1)-y(2)*(y(1)+1)+0.00243*y(1)*(y(1)+1)];
f=0:0.001:0.2;
hold on
for i=1:length(f)
sol=dde23(@der2,[2.7],[-0.1,-0.1],[0:0.001:10],[],f(i));
y1=sol.y(1,:);
plot(f(i),y1(1:10:end),'k.')
end
我想画的是关于反馈系数的分岔图，模拟了一个例子，但是报错，有没有人知道怎么改？