Search code examples
wolfram-mathematicafixed-point-iteration

Finding the Fixed Points of an Iterative Map


I need to find fixed points of iterative map x[n] == 1/2 x[n-1]^2 - Mu.
My approach:

Subscript[g, n_ ][Mu_, x_] :=  Nest[0.5 * x^2 - Mu, x, n]

fixedPoints[n_] := Solve[Subscript[g, n][Mu, x] == x, x]

Plot[
  Evaluate[{x, 
   Table[Subscript[g, 1][Mu, x], {Mu, 0.5, 4, 0.5}]}
  ], {x, 0, 0.5}, Frame -> True]

Solution

  • I'll change notation slightly (mostly so I myself can understand it). You might want something like this.

    y[n_, mu_, x_] := Nest[#^2/2 - mu &, x, n]
    fixedPoints[n_] := Solve[y[n, mu, x] == x, x]
    

    The salient feature is that the "function" being nested now really is a function, in correct format.

    Example:

    fixedPoints[2]
    
    Out[18]= {{x -> -1 - Sqrt[-3 + 2*mu]}, 
              {x -> -1 + Sqrt[-3 + 2*mu]}, 
              {x ->  1 - Sqrt[ 1 + 2*mu]}, 
              {x ->  1 + Sqrt[ 1 + 2*mu]}}
    

    Daniel Lichtblau