Logo by bib - Contribute your own Logo!

END OF AN ERA, FRACTALFORUMS.COM IS CONTINUED ON FRACTALFORUMS.ORG

it was a great time but no longer maintainable by c.Kleinhuis contact him for any data retrieval,
thanks and see you perhaps in 10 years again

this forum will stay online for reference
News: Check out the originating "3d Mandelbulb" thread here
 
*
Welcome, Guest. Please login or register. April 24, 2024, 06:59:38 PM


Login with username, password and session length


The All New FractalForums is now in Public Beta Testing! Visit FractalForums.org and check it out!


Pages: [1]   Go Down
  Print  
Share this topic on DiggShare this topic on FacebookShare this topic on GoogleShare this topic on RedditShare this topic on StumbleUponShare this topic on Twitter
Author Topic: Appearances can be deceptive  (Read 2857 times)
0 Members and 1 Guest are viewing this topic.
element90
Strange Attractor
***
Posts: 298



WWW
« on: May 11, 2015, 12:29:17 PM »

The following formula produces something that appears to be a standard Mandelbrot as part of a larger fractal. It turns out not to be.

z = \frac{z^2+c}{\frac{1}{z^6}+\frac{1}{c}}+c+7\

The initial value of z is zero.






Clearly zero is a critical point but I don't see how it is derived:

f'(z) = \left( \frac{g(z)}{h(z)} \right)' = \frac{g'(z)h(z)-g(z)h'(z)}{[h(z)]^2}

g(z) = z^2+c
g'(z) = 2z
h(z) = \frac{1}{z^6}+\frac{1}{c}
h'(z) = -\frac{6}{z^7}

f'(z) = \frac{2z\left( \frac{1}{z^6}+\frac{1}{c} \right)-\left( z^2+c \right)\left( -\frac{6}{z^7} \right)}{[\frac{1}{z^6}+\frac{1}{c}]^2} = 0

Multiplying both sides by h(z) squared simplifies the formula to:

\frac{1}{z^7}\left( 8z^2+\frac{2z^8}{c}-6c \right) = 0

The polynomial can be solved for zero and will give location dependent solutions and hence critical points. How can zero be a solution, as it must be if it is a critical point?
Logged

Elelemt90 Fractals blog www.element90.wordpress.com
lkmitch
Fractal Lover
**
Posts: 238



« Reply #1 on: May 11, 2015, 05:59:06 PM »

Two things:  1) a critical point is a spot where the derivative = 0 or fails to exist.  In your last equation, the 1/z^7 factor is what makes 0 a critical point, since you'd be dividing by zero.  2) If you rewrite the first term of f(z) as a rational function (that is, one numerator polynomial and one denominator polynomial), you'll see in its derivative that every term has some extra z factors in the numerator, so the net effect is that z=0 would cause the derivative to be zero.
Logged
Adam Majewski
Fractal Lover
**
Posts: 221


WWW
« Reply #2 on: May 11, 2015, 06:02:41 PM »

IMHO one has to compute critical point for every point c

(%i1) f:(z^2+c)/(z^{-6} +c^{-1}) +c+7;
(%o1) (z^2+c)/(z^{−6}+c^{−1})+c+7
(%i2) d:diff(f,z,1);
(%o2) (2*z)/(z^{−6}+c^{−1})−({−6}*z^{−6}−1*(z^2+c))/(z^{−6}+c^{−1})^2
(%i3) solve(d=0);
solve: more unknowns than equations.
Unknowns given : 
[c,z]
Equations given: 
errexp1
 -- an error. To debug this try: debugmode(true);
(%i4) expand(d);
(%o4) −({−6}*z^({−6}+1))/(z^(2*{−6})+2*c^{−1}*z^{−6}+c^(2*{−1}))−({−6}*c*z^{−6}−1)/(z^(2*{−6})+2*c^{−1}*z^{−6}+c^(2*{−1}))+(2*z)/(z^{−6}+c^{−1})

Compare https://en.wikibooks.org/wiki/Fractals/Iterations_in_the_complex_plane/qpolynomials

HTH
Logged
element90
Strange Attractor
***
Posts: 298



WWW
« Reply #3 on: May 11, 2015, 08:42:24 PM »

Quote
IMHO one has to compute critical point for every point c

I know that's why I stated:

Quote
The polynomial can be solved for zero and will give location dependent solutions

i.e. the critical points have to determined for each location, I already do this in Neptune and in the new version of Saturn.

Quote
1) a critical point is a spot where the derivative = 0 or fails to exist.

I clearly know about the solutions to f'(z) = 0 but what does "fails to exist" mean?

Quote
2) If you rewrite the first term of f(z) as a rational function (that is, one numerator polynomial and one denominator polynomial), you'll see in its derivative that every term has some extra z factors in the numerator, so the net effect is that z=0 would cause the derivative to be zero.

I know how to do that.
Logged

Elelemt90 Fractals blog www.element90.wordpress.com
Pages: [1]   Go Down
  Print  
 
Jump to:  

Related Topics
Subject Started by Replies Views Last post
demoscene - fractal appearances Art Discussions Oppenheimer_ 7 4327 Last post July 06, 2011, 01:05:57 PM
by Oppenheimer_

Powered by MySQL Powered by PHP Powered by SMF 1.1.21 | SMF © 2015, Simple Machines

Valid XHTML 1.0! Valid CSS! Dilber MC Theme by HarzeM
Page created in 0.111 seconds with 24 queries. (Pretty URLs adds 0.005s, 2q)