Thread Rating:
  • 1 Vote(s) - 5 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Iteration exercises: f(x)=x^2 - 0.5 ; Fixpoint-irritation...
#23
(10/20/2017, 07:55 PM)Gottf ried Wrote:
(10/20/2017, 06:00 PM)sheldonison Wrote: For example, here is the pari-gp program for the formal inverse schroeder function.  I don't know how to turn this into a matrix function, but not many programming languages support the powerful polyonomial functions that pari-gp has.

Code:
formalischroder(fx,n) = {
  local(lambda,i,j,z,f1t,f2t,ns,f1s);
  lambda = polcoeff(fx,1);
  f1t=x;
  i=2;
  while (i<=n,
    f1s=f1t;
    f1t=f1t+acoeff*x^i+O(x^(i+1));
    f2t=subst(f1t,x,lambda*x)-subst(fx+O(x^(i+1)),x,f1t);
    z = polcoeff(f2t, i);
    z = subst(z,acoeff,x);
    ns=-polcoeff(z,0)/polcoeff(z,1);
    f1t=f1s+ns*x^i;
    i++;
  );
  return(Pol(f1t));
}
fz1=x^2+(1-sqrt(3))*x;
[size=small][font=Monaco, Consolas, Courier, monospace]lambda1=polcoeff(fz1,1);[/font][/size]
[size=small][font=Monaco, Consolas, Courier, monospace]fs1=formalischroder(fz1,20);[/font][/size]
superfunction1(z)=subst(fs2,x,lambda2^z);

Sheldon - I find some unexplained terms: what is "acoeff" / how is this defined before it is queried for "f1t"?

the same with "fz2" in "lambda1=pol..." and "fs2" in "superfuncion...subst(fs2..."
acoeff is just like an unknown; like "x" in the equations "x^3+2*x^2+3*x+1".  In the routine, Pari-gp will treat acoeff as an unknown unassigned variable just like "x", that happens to be the multiplier for the "i'th" coefficient we are trying to determine.  

The others are typos:  I was originally going to post both superfunctions, from both fixed points....  The other fixed point equation would be fz2=x^2+(1+sqrt(3))*x.  Depending on whether |lambda|>1 or <1, determines where the superfunction equation converges, either negative or positive values of .

Actually, the other fixed point also converges better, so that's what I meant to post.  Here superfunction1 converges pretty good for
Code:
fz1=x^2+(1+sqrt(3))*x;
lambda1=polcoeff(fz1,1);
fs1=formalischroder(fz1,20);
superfunction1(z)=subst(fs1,x,lambda1^z);
- Sheldon
Reply


Messages In This Thread
RE: Iteration exercises: f(x)=x^2 - 0.5 ; Fixpoint-irritation... - by sheldonison - 10/20/2017, 08:19 PM

Possibly Related Threads...
Thread Author Replies Views Last Post
  (Again) fixpoint outside Period tommy1729 2 2,630 02/05/2017, 09:42 AM
Last Post: tommy1729
  Polygon cyclic fixpoint conjecture tommy1729 1 2,349 05/18/2016, 12:26 PM
Last Post: tommy1729
  The " outside " fixpoint ? tommy1729 0 1,577 03/18/2016, 01:16 PM
Last Post: tommy1729
  2 fixpoint pairs [2015] tommy1729 0 1,874 02/18/2015, 11:29 PM
Last Post: tommy1729
  [2014] The secondary fixpoint issue. tommy1729 2 3,806 06/15/2014, 08:17 PM
Last Post: tommy1729
  Simple method for half iterate NOT based on a fixpoint. tommy1729 2 3,533 04/30/2013, 09:33 PM
Last Post: tommy1729
  Iteration exercises: Lucas-Lehmer-test and Schröder-function Gottfried 0 2,768 04/04/2012, 06:17 AM
Last Post: Gottfried
  Iteration series: Different fixpoints and iteration series (of an example polynomial) Gottfried 0 2,979 09/04/2011, 05:59 AM
Last Post: Gottfried
  Fractional iteration of x^2+1 at infinity and fractional iteration of exp bo198214 10 17,287 06/09/2011, 05:56 AM
Last Post: bo198214
  2 fixpoint failure tommy1729 1 2,969 11/13/2010, 12:25 AM
Last Post: tommy1729



Users browsing this thread: 1 Guest(s)