Group Abstract Group Abstract

Message Boards Message Boards

0
|
4.1K Views
|
4 Replies
|
4 Total Likes
View groups...
Share
Share this post:
GROUPS:

FixedPoint

Posted 10 years ago

How to find Fixedpoint of the famous Mackey Glass Equation using Mathematica code

x'[t] == (1/4) x[t - 15]/ (1 + x[t - 15]^10) - x[t]/10
POSTED BY: Dia Ghosh
4 Replies
Posted 10 years ago

Thank you once again ? I have one more question, you have used the command Findroot, can we use the command Fixedpoint? Is there any deference between root and fixed point?

POSTED BY: Dia Ghosh

Hi,

same idea. Again the right hand side has to be zero. The Bessel function does not make this easier, but it is possible to get this numerically.

It is a good idea to first plot the right hand side - setting x[t-15] and x[t] simply to x.

Plot[(1/4) BesselJ[1, x] - x/10, {x, -5, 5}]

enter image description here

So, it looks as if there are three zeros: one at 0 one close to 1 and one close to -1.

We can get that in one line:

FindRoot[(1/4) BesselJ[1, x] - x/10 == 0, {x, {-1, 0, 1}}]

where the -1,0,1 are the starting values for our search. We obtain:

({x -> {-1.31102, 0., 1.31102}})

so three fixed points. If you are unhappy about guessing the first points, you can start at lots of random starting points:

ToExpression /@ DeleteDuplicates[ToString /@ Chop[FindRoot[(1/4) BesselJ[1, x] - x/10 == 0, {x, #}]] & /@ RandomReal[{-10, 10}, 100]]
(*{{x -> -1.31102}, {x -> 0}, {x -> 1.31102}}*)

The ToString and ToExpression step is not nice, but otherwise DeleteDuplicates does not work.

Cheers,

Marco

POSTED BY: Marco Thiel
Posted 10 years ago

Thank you for your reply. Now instead of x[t - 15] if anybody has BesselJ[1, x[t - 15]] then ? I mean x'[t] == (1/4) BesselJ[1, x[t - 15]] - x[t]/10

POSTED BY: Dia Ghosh
POSTED BY: Marco Thiel
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard