Code covered by the BSD License  

Highlights from
Tutorial on solving DDEs with DDE23

  • exam1
  • exam1
  • exam2
  • exam2
  • exam3
  • exam3
  • exam4
  • exam4
  • exam5
  • exam5
  • exam6 This is a demonstration problem for CTMS/BD in
  • exam6 This is a demonstration problem for CTMS/BD in
  • exam7 An example from C. Marriott and C. DeLisle, Effects
  • exam7 An example from C. Marriott and C. DeLisle, Effects
  • exam8 This is the suitcase problem from Suherman, et al.,
  • exam8 This is the suitcase problem from Suherman, et al.,
  • exer1 Example 1 of K.W. Neves, Automatic integration
  • exer1 Example 1 of K.W. Neves, Automatic integration
  • exer2 Example of J.D. Farmer, Chaotic Attractors of an
  • exer2 Example of J.D. Farmer, Chaotic Attractors of an
  • exer3 Wheldon's model of chronic granuloctic leukemia
  • exer3 Wheldon's model of chronic granuloctic leukemia
  • exer5
  • exer5
  • exer6 Sample problem of ARCHI manual. The absolute error
  • exer6 Sample problem of ARCHI manual. The absolute error
  • exer7 Marchuk immunology model of E. Hairer, S.P. Norsett, and
  • exer7 Marchuk immunology model of E. Hairer, S.P. Norsett, and
  • prob1 This system of ODE's is taken from 'An Introduction to Nuermcial Methods
  • prob1 This system of ODE's is taken from 'An Introduction to Nuermcial Methods
  • prob2 This problem considers a cardiovascular model, which can be found in
  • prob2 This problem considers a cardiovascular model, which can be found in
  • prob2b This problem considers a cardiovascular model, which can be found in
  • prob2b This problem considers a cardiovascular model, which can be found in
  • prob3 This problem is epidemic model due to Cooke, more information can be
  • prob3 This problem is epidemic model due to Cooke, more information can be
  • prob4 This problem is an epidemic model due to Cooke et alia, more information
  • prob4 This problem is an epidemic model due to Cooke et alia, more information
  • prob5 This problem population growth model due to Cooke et alia, more information
  • prob5 This problem population growth model due to Cooke et alia, more information
  • View all files

Tutorial on solving DDEs with DDE23

by

 

22 Aug 2003 (Updated )

Solving delay differential equations with DDE23. Tutorial + Examples.

prob1
function [sol0,sol1] = prob1
%  This system of ODE's is taken from 'An Introduction to Nuermcial Methods
%  for Differential Equations', by J.M. Ortega and W.G. Poole (Reference 13 
%  of the tutorial).  There is additional information about predator-prey 
%  systems in 'Functional Differntial Equations' by J. Hale (Reference 7 of
%  the tutorial).

% Copyright 2002, The MathWorks, Inc.

history = [80; 30];
tspan = [0, 100];
opts = ddeset('RelTol',1e-5,'AbsTol',1e-8);
a =  0.25;
b = -0.01;
c = -1.00;
d =  0.01;

% Solve the ODEs that arise when there is no delay.
sol0 = dde23(@prob1f,[],history,tspan,opts,a,b,c,d);
   
% Solve the DDEs that arise when there is a delay of tau.
tau = 1;
sol1 = dde23(@prob1f,tau,history,tspan,opts,a,b,c,d);

figure
plot(sol0.y(1,:),sol0.y(2,:),sol1.y(1,:),sol1.y(2,:))
title('Problem 1. Solution with and without delay.')
xlabel('y_1(t)')
ylabel('y_2(t)')
legend('No delay',['Delay \tau = ',num2str(tau)],2)

%-----------------------------------------------------------------------

function v = prob1f(t,y,Z,a,b,c,d)
%PROB1F  The derivative function for Problem 1 of the DDE Tutorial.
v = zeros(2,1);
if isempty(Z)     % ODEs
   v(1) = a * y(1) + b * y(1) * y(2);
   v(2) = c * y(2) + d * y(1) * y(2);
else              % DDEs
   m = 200;
   ylag = Z(:,1);
   v(1) = a * y(1) * (1 - y(1) / m) + b * y(1)    * y(2);
   v(2) = c * y(2)                  + d * ylag(1) * ylag(2);
end

Contact us