function dydt = dderhs(t,y,z)
dydt = 2*z/(1+z^5)-y;

end

% Main:
% res = dde23(@dderhs,2,0.5,[0,300])
% plot(sol.x,sol.y)