Re: Mathematica 6.0 Desperation
- To: mathgroup at smc.vnet.net
- Subject: [mg83686] Re: Mathematica 6.0 Desperation
- From: congruentialuminaire at yahoo.com
- Date: Wed, 28 Nov 2007 05:35:24 -0500 (EST)
- References: <figtq9$fa7$1@smc.vnet.net>
Hello Jesus:
I have seen myriad examples "hang" V6. Some I have debugged, many
others have not had my attention yet.
I changed the code above to this:
Clear[JuliaTest];
JuliaTest[x_, y_, n_, c_] := Module[{z, num = 0},
z = x + y I;
While[Abs[z] < 2.0 && num < n, z = z^2 + c; num++];
num
]
Then this:
DensityPlot[
JuliaTest[x, y, 50, -0.75 + 0.5 I], {x, 1.24, 1.27}, {y, -0.13,
-0.1}, PlotPoints -> 50, Mesh -> False, ColorFunction -> Hue]
Then the "Julia Set" is shown.
I did not check much to see why the old code does not work. Note that
I removed Complile[].
HTH
Regards..Roger W.
On Nov 27, 3:07 am, "Jesus M. Goiri" <jesgoir... at mac.com> wrote:
> Dear Sirs
>
> I have tried some items that seem to hang Mathematica 6.0. However
> they work perfectly fine in the 5.2 version. I have tried to
> discover the flaw and I have not been able to do so. The problem
> pertains at some code I extracted from a book about Mathematica
> applications by Boccara, particularly the examples that relate to
> fractals, Julia and Mandelbrot Sets.
>
> I have below the lines of code. I you try them, I'd wish to see if
> you are able to make them work.
>
> Yours sincerely,
>
> Jesus M, Goiri
>
> Clear[JuliaTest]
> JuliaTest = Compile[{x, y, {n, _Integer}, {c, _Complex}},
> Module[{z, num = 0}, z = x + y I;
> While[Abs[z] < 2.0 && num < n, z = z^2 + c; num++]; num]]
> DensityPlot[JuliaTest[x, y, 50, -0.75 + 0.5 I],
> {x, 1.24, 1.27}, {y, -0.13, -0.1}, PlotPoints -> 500,
> Mesh -> False, ColorFunction -> Hue]