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
from Tutorial on solving DDEs with DDE23 by Jacek Kierzenka
Solving delay differential equations with DDE23. Tutorial + Examples.

exam6 
function sol = exam6
% This is a demonstration problem for CTMS/BD in 
% L. Tavernini, Continuous--Time Modeling and 
% Simulation, Gordon and Breach, Amsterdam, 1996.  
% It illustrates how to deal with the unusual
% situation of a solution discontinuous at the 
% initial point. Tavernini uses an initial y of 
% 19.00001, but we use 19.001 so as to see the 
% cyclic behavior sooner.

% Copyright 2004, The MathWorks, Inc.

  % Known parameters, visible in nested functions.
  r = 3.5;
  m = 19;

  options = ddeset('RelTol',1e-4,'AbsTol',1e-7,...
                   'InitialY',19.001);
  sol = dde23(@exam6f,0.74,19,[0, 40],options);
  figure
  plot(sol.x,sol.y);
  title('Population of Lemmings--Time Series')
  xlabel('time t');
  ylabel('y(t)');
  figure
  plot(sol.y,sol.yp)
  title('Population of Lemmings--Phase Plane')
  xlabel('y(t)');
  ylabel('y''(t)');
  
  %-----------------------------------------------------------------------
  % Nested function
  %
  
  function yp = exam6f(t,y,Z)
  %EXAM6F  The derivative function for the Example 6 of the DDE Tutorial.
    yp = r*y*(1 - Z/m);
  end % exam6f

  %-----------------------------------------------------------------------
    
end  % exam6

Contact us