fractal package



On 11/26/2012 11:09 PM, Doug Stewart wrote:
> Now for bfloats
>
> plot3d (logistic_set, [x, 3.850000000000000000001b0, 
> 3.850000000000000000009b0], [y, -.000005, .000005],
>           [gnuplot_preamble, "set view map"],
>                  [gnuplot_pm3d, true],
>                  [grid, 150, 150])$
>
> plotting: range must be of the form [variable, min, max]; found: [x, 
> 3.85b0, 3.85b0]
>  -- an error. To debug this try: debugmode(true);
>
> How do I make this work with bfloats????
>
>
plot3d does not accept bfloats. For you purpose I think you want to do 
your calculations with a very high precision, and then create for 
instance a matrix that will represent the pixels in the plot you want to 
make. The bfloat numbers that you compute with your logisitic_set 
function should then be rounded to integers that will label the rows and 
columns in the matrix, and the color for the pixel. I think that using 
rational numbers would be better that using bfloats, but you can used 
either.

You can use as example  the function mandelbrot, that I wrote before the 
fractals module was contributed and whose source code you will find in 
share/dynamics/complex_dynamics

(? mandelbrot will show you some documentation and you can see an 
example in ? Functions and Variables for dynamics)

Regards,
Jaime