Raymond Toy pisze:
> Adam Majewski wrote:
>> Hi,
>>
>> I use bfallroots function :
>> https://apps.sourceforge.net/mediawiki/maxima/index.php?title=Image:Centers_9_new.png
>> Can it be done for higher periods ? ( for me (GCL) it fails for period
>> 10 and fpprec:150 or 256)
>>
> Assuming I did this right, I see that the polynomial for period 10 is of
> order 495. The largest coefficient has about 285 bits in it. So maybe
> you need to use fpprec of 300 or more?
I have changed to :
elseif p=10 then fpprec:300,
but there is error after computing numbers :
(%i11) for period from 9 thru period_Max do
(centers[period]:GiveCenters_bf(period),N_of_centers:length(centers[period])+N_of_centers)
(%o11) done
(%i12) load(draw)(%o12)
C:/PROGRA~1/MAXIMA~1.3/share/maxima/5.16.3/share/draw/draw.lisp
(%i13)
draw2d(file_name=centers_10,pic_width=1000,pic_height=1000,yrange=[-1.5,1.5],xrange=[-2.5,0.5],
title=concat(centers of,string(N_of_centers), hyperbolic components of
Mandelbrot set for periods 1-,string(period_Max)),
user_preamble=set size square;set key out;set key top
left,terminal='png,pic_width=1000,pic_height=1000,xlabel=re ,
ylabel=im,point_type=filled_circle,points_joined=false,key=9,color=gray,points(map(realpart,centers[9]),map(imagpart,centers[9])),
key=10,color=purple,points(map(realpart,centers[10]),map(imagpart,centers[10])))
Maxima encountered a Lisp error: Error in FORMAT [or a callee]: The
storage for STREAM is exhausted.Currently,
205 pages are allocated.Use ALLOCATE to expand the space.Automatically
continuing.To reenable the Lisp debugger set *debugger-hook* to nil.
I think that it is coused by bf number format :
9.4630276158882190662033636749[245 digits]42308247717006463653684639b-91*%i
-1.9959241519307532077928469860[243 digits]8039417293485824429601363747b0
I have prpbably to remove "[245 digits]" to get number
Adam