function sol=simpledde
tau=pi/2+0.1;
options=ddeset('RelTol', 1e-5, 'AbsTol', 1e-5);
sol = dde23(@ddede, [tau], @ddehist, [0,200],options);   
plot(sol.x, sol.y);

function y = ddehist(t)
   y=[0.1];
    
function dydt = ddede(t,y, yl)

  dydt = -yl-y^3;

