The Museum of HP Calculators The Museum of HP Calculators
Copyright (C) 2002 Glen Kilpatrick
Distributed under GNU General Public License
This program is supplied without representation or warranty of any kind. The author and The Museum of HP Calculators therefore assume no responsibility and shall have no liability, consequential or otherwise, of any kind arising from the use of this program material or any part thereof.
I hardly need introduce Benoit B. Mandelbrot's fascinating discovery. A Google search on "Mandelbrot Set" turned up "about 35,900" hits as I was writing this; many of them will probably give a better introduction than I can, so I'm going to completely pass over what it is, and just give code to compute it.
I sold my HP-85B to Museum of HP Calculators years ago, so may not be very useful for support, but you're always welcome to ask....
None given. This is a "hardcoded" example to be executed directly. Exploration of the Mandelbrot set is done by varying the DATA statement at line 210, then re-executing the program.
10 ! HP-85B (Benoit B.) Mandelbrot Set (ported from HP-28S 6→19SEP92)
20 ! Copyright 1992 / Glen A. Kilpatrick / All Rights Reserved
30 ! 1:1 aspect for 132*73 HP-28S is (-6.1015625,-1.25),(4.6015625,1.25)
40 CLEAR @ INTEGER B,N.P,R ! ?0/1/8/9 ≤> Real/Imaginary/Lower/Upper
50 T8=TIME @ D8=DATE @ READ X8,X9,Y8,Y9,B,R,P ! 1:1 aspect for 256*192
60 PEN P @ A=(X9-X8)/(Y9-Y8) @ A=A*192/256 @ IF A≥1 THEN 70 ELSE 80
70 M=Y8+(Y9-Y8)/2 @ D=M-Y8 @ D=D*A @ Y8=M-D @ Y9=M+D @ GOTO 90 ! *H
80 M=X8+(X9-X8)/2 @ D=M-X8 @ D=D/A @ X8=M-D @ X9=M+D ! INV *W
90 DISP INT(.5+SQR(100/9+6.25)/SQR((X9-X8)^2+(Y9-Y8)^2));"magnification"
100 DISP "PMIN=(";X8;",";Y8;")" @ DISP "PMAX=(";X9;",";Y9;")"
110 DISP "RES=";R;"BntB=";B @ WAIT 99999 @ SCALE X8,X9,Y8,Y9 @ GCLEAR
120 FOR C0=X8 TO X9 STEP R*(X9-X8)/256 ! outer loop runs left → right
130 FOR C1=Y8 TO Y9 STEP R*(Y9-Y8)/192 ! inner loop runs bottom → top
140 Z0,Z1=0 @ FOR N=1 TO B @ T0=Z0 @ Z0=Z0*Z0-Z1*Z1+C0 @ Z1=2*T0*Z1+C1
150 IF Z0*Z0+Z1*Z1>4 THEN 170 ! set exclusion loop (|R|>2 => not in set)
160 NEXT N @ PENUP @ PLOT C0,C1 @ PENUP ! set inclusion, plot this point
170 NEXT C1 @ NEXT C0 @ D9=DATE @ T9=TIME ! done; show hours, toggle axes
180 DISP (D9-D8)*24+INT(.5+10*((24*(T9-T8)/3600)MOD24))/10;"hours"
190 XAXIS Y8,(X9-X8)/10 @ YAXIS X9,(Y9-Y8)/10 @ P=-P @ PEN P
200 YAXIS X8,(Y9-Y8)/10 @ XAXIS Y9,(X9-X8)/10 @ WAIT 99999 @ GOTO 190
210 DATA -2.41666666667,.91666666667,-1.25,1.25,100,1,-1
220 END
Go back to the software library