I want to draw a bifurcation diagram for system of four delay differential equations

clc

lags = 0.6;

sol = dde23(@ddex1de,lags,@ddex1hist,[0,8]);

t = sol.x;

y = sol.y;

plot(t,y(1,:));

function dydt = ddex1de(t,y,Z)

ylag1 = Z(:,1);

k1 = 2;

k_1 = 1;

k2 = 1;

k3 = 1;

k_3 = 1;

k4 = 2;

k5 = 1;

E1 = 1;

E2 = 1;

dydt = [-k1*y(1)*E1+k_1*y(2)+k4*y(3)-k5*ylag1(1);

k1*y(1)*E1-(k_1+k2)*y(2);

k3*y(4)*E2-(k_3+k4)*y(3);

k2*y(2)+k_3*y(3)-k3*y(4)*E2+k5*ylag1(1)];

end

function S = ddex1hist(t)

S = [5;0;0;0];

end