a=126.2675455;
b=118.6876595;
tau=1.0;

history = [0.000; 0.00085]; %initial function
tspan = [0, 20]; 
opts = ddeset('RelTol',1e-5,'AbsTol',1e-7);

% Solve equation
sol = dde23('logdist', tau, history, tspan, opts, a, b);

plot(sol.x, sol.y(2,:), 'k');
%legend('x(t)=\int^t_0 y(s)ds', 'y(t)',2);
text(2,0.0025,'Initial condition is y(\theta)=0.00085, \theta\in [-1, 0]','FontSize',10)
title('Periodic solution of y\prime=\sigma y(1-y-\nu\int^1_0 y(t-s)ds)-\nuy;');
xlabel('time t');
ylim([0 0.003]);
ylabel('y', 'rotation', 0);