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