带时滞微分方程分岔图

请问matlab怎么画含有时滞的二阶微分方程的分岔图？
照猫画虎也没弄对

f=0:0.01:5;
tt=0.05;tspan=[-5 5];
options=odeset('RelTol',1e-7);
hold on
for j=1:length(f)
    sol=dde23(@fc,lags,ones(2,1),tspan,options);
    i=500:100:1500
plot(f(j),x(i,1),'k.','LineWidth',10)
end

function dydt=fc(t,x,z)    %函数主体
 global f
dydt=[x(2);-0.2*x(2)+x(1)-0.5*(x(1)^3)+f*cos(t)+0.05*z(1,1)];
 end