function sol =dglv1



sol = dde23(@dglv, [50],ones(3,1),[0,5]);
figure
plot(sol.x,sol.y);



 
end