Services & ResourcesWolfram ForumsMathGroup Archive

[Date Index] [Thread Index] [Author Index]

question about NDelayDSolve

• To: mathgroup at smc.vnet.net
• Subject: [mg69703] question about NDelayDSolve
• From: aitor69gonzalez at gmail.com
• Date: Thu, 21 Sep 2006 07:29:36 -0400 (EDT)

Hello,

I am using NDelayDSolve to solve a set of delay differential equations
that oscillate (see below). I have noticed that the range of
integration (either limit1=400 or limit2=1000) of this function exerts
a big influence on the result. This leads me to ask you, which is the
best result if any? Could you give me some tips on how to use this
function rigorously?
Thank you in advance.

This is the code:

SetDirectory["/your/mathematica/package/directory"];
<< NDelayDSolve.m
a = 4.5; b = 0.23; c = 0.23; k = 33; p0 = 40; Tm = 15; Tp = 2; limit1 =
400; \
limit2 = 1000;
range1 = NDelayDSolve[{p'[t] == a*m[t - Tp] - b*p[t],
        m'[t] == k/(1 + (p[t - Tm]^2/p0^2)) - c*m[t]}, {m -> (0 &),
        p -> (0 &)}, {t, 0, limit1}, AccuracyGoal -> 0, PrecisionGoal
-> 0];
range2 = NDelayDSolve[{p'[t] == a*m[t - Tp] - b*p[t],
        m'[t] == k/(1 + (p[t - Tm]^2/p0^2)) - c*m[t]}, {m -> (0 &),
        p -> (0 &)}, {t, 0, limit2}, AccuracyGoal -> 0, PrecisionGoal
-> 0];
Plot[{Evaluate[p[t] /. range1], Evaluate[p[t] /. range2]}, {t, 0,
limit1}, 
  PlotStyle -> {RGBColor[0, 1, 0], RGBColor[1, 0, 0]}]

• Follow-Ups:
  • Re: question about NDelayDSolve
    • From: "Chris Chiasson" <chris@chiasson.name>