• 0 Vote(s) - 0 Average
• 1
• 2
• 3
• 4
• 5
 Find all fixed points of exp[b] MorgothV8 Junior Fellow Posts: 18 Threads: 6 Joined: Dec 2012 08/30/2014, 01:41 PM (This post was last modified: 08/30/2014, 02:48 PM by MorgothV8.) Hi, maybe too easy but I have problems. I waant to find set of all fixed points for complex z^x. I know that I can do: L = fixed_point(b^z) for all b in Complex plane. I'm doing: fix_point(b) := Lambert_W(a), where a = -ln(b) Labmert_W(a) := Tania(t), where t = ln(a)-1 Tania(t): 4 different implementations (after citizendium), depending on t values (imag/real). I can get fixed points quite OK, for: b=e, b=sqrt(2), b=2. I'm computing: L=fix_point(b) fL = f(L) = pow(b, L) Then I check |L-fL| and get values < 1e-19, which is OK. My only problem is: I want ALL fixed points for ALL complex bases b f(b) = set(L), set(L) = set of fixed points not just one. Problem is that for b=sqrt(2) i get L=4, which is OK, but I want to get L=2 --> how to do that? Anybody will help? My goal is to have complex map of f(z) = set { fixed points of pow(z) } I mean I want multivalued mapping. EDIT: It behaved really bad for bases near e^(1/e), specially equal and little below. This was due to log of values with negative real part and almost zero imag part - which returned really bad results (because log of negative real is cut line), but more itarations helped. Code updated to "adjust" iterations untill they behave OK. So there is limit 100 iterations, but after each one I check added value (in Tania function), if value drops |s| < EPS (1e-16) then abort iterations and return Tania result. This helped for e^(1/e) and its vicinity but.... there is NO convergence of Tania function (this is the one from citizendium) when b = PI/3 (I wanted some special PI related value, and less than e^(1/e) and I tested PI/3 and Tania is not converging) Checked such values: //z = C(-1.0, 0.); //z = C(0., 0.); //z = C(1., 0.); //z = C(exp(1.), 0.); //z = C(exp(2.), 0.); //z = C(exp(-exp(1.)), 0.0); //z = C(exp(1./exp(1.)), 0.0); //z = C(exp(1./exp(1.)-0.00000001), 0.0); //z = C(exp(1./exp(1.)+0.00000001), 0.0); //z = C(sqrt(2.), 0.0); //z = C(M_PI/3., 0.0); Updated code Attached Files   math4.zip (Size: 1.54 KB / Downloads: 632) Fuji GSW690III Nikon D3, Nikkors 14-24/2.8, 24/1.4, 35/2, 50/1.4, 85/1.4, 135/2, 80-200/2.8 « Next Oldest | Next Newest »

 Messages In This Thread Find all fixed points of exp[b] - by MorgothV8 - 08/30/2014, 01:41 PM RE: Find all fixed points of exp[b] - by MorgothV8 - 09/17/2014, 04:17 PM RE: Find all fixed points of exp[b] - by jaydfox - 09/17/2014, 11:10 PM RE: Find all fixed points of exp[b] - by MorgothV8 - 09/18/2014, 04:31 AM RE: Find all fixed points of exp[b] - by MorgothV8 - 09/24/2014, 03:53 PM RE: Find all fixed points of exp[b] - by jaydfox - 09/26/2014, 04:47 PM RE: Find all fixed points of exp[b] - by mike3 - 09/29/2014, 01:49 AM RE: Find all fixed points of exp[b] - by mike3 - 09/29/2014, 02:41 PM RE: Find all fixed points of exp[b] - by mike3 - 10/07/2014, 01:32 AM RE: Find all fixed points of exp[b] - by mike3 - 10/07/2014, 01:50 AM RE: Find all fixed points of exp[b] - by Gottfried - 10/07/2014, 11:00 AM

 Possibly Related Threads... Thread Author Replies Views Last Post Kneser-iteration on n-periodic-points (base say \sqrt(2)) Gottfried 11 4,290 05/05/2021, 04:53 AM Last Post: Gottfried On n-periodic points of the exp() - A discussion with pictures and methods Gottfried 1 2,822 06/10/2020, 09:34 AM Last Post: Gottfried fixed point formula sheldonison 6 18,755 05/23/2015, 04:32 AM Last Post: mike3 Attempt to find a limit point but each step needs doubling the precision... Gottfried 15 33,686 11/09/2014, 10:25 PM Last Post: tommy1729 sexp(strip) is winding around the fixed points Kouznetsov 8 20,700 06/29/2009, 10:05 AM Last Post: bo198214 An error estimate for fixed point computation of b^x bo198214 0 4,092 05/31/2008, 04:11 PM Last Post: bo198214

Users browsing this thread: 1 Guest(s)