a = 0.8;
b = 0.7;
r = 0.08;

% Compute steady state solution [v0; w0] for m = 0:
v0 = 2.417960226013935;
w0 = (v0 + a) / b;
history = [0.4*v0 ; 1.8*w0];

for m = [+10 -10]
   sol = dde23('prob3f',20,history,[0, 60],[],v0,a,b,r,m);
   figure
   plot(sol.x,sol.y)
   title(['Problem 3. Plant''s Neuron Interaction Model' ...
          ' with m = ',num2str(m),'.'])
   xlabel('time t')
end