Re: Fractals with Mathematica
- To: mathgroup at smc.vnet.net
- Subject: [mg34306] Re: Fractals with Mathematica
- From: shollis at armstrong.edu (Selwyn Hollis)
- Date: Tue, 14 May 2002 04:11:17 -0400 (EDT)
- Organization: http://groups.google.com/
- References: <abdg2n$chk$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Shz Shz, Here is a "simple" example using Newton's method: newtstep[x_]:= x - (x^3-x)/(3x^2-1); closeQ[xold_, xnew_] := Abs[xold-xnew] < 0.0005; newton[x0_]:= FixedPoint[newtstep, x0, 30, SameTest->closeQ]; xmin = .4; xmax = .65; nx := 200; dx = (xmax - xmin)/nx; ymin = -.05; ymax = .05; ny := 300; dy = (ymax - ymin)/ny; lims = Flatten[Transpose[Abs[Table[newton[x + y*I], {x, xmin, xmax, dx}, {y, ymin, ymax, dy}]]]]; Show[Graphics[RasterArray[Partition[Map[Hue,0.67 + 0.33lims], nx+1]]]] --Selwyn