[comp.sys.apple] BASIC Mandelbrot program for regular hi-res

LMB7421@RITVAX.BITNET (03/26/88)

Standard values for the set are Rmin (Pmin) = -2.25
                                Rmax (Pmax) = 0.75
                                Imin (Qmin) = -1.5
                                Imax (Qmax) = 1.5
Choose other values from the boundary region of the set (try to read them
from the screen, get them from the Peitgen book, or send me e-mail)

--------------------------------------------------------------------
0  REM Freeware program  see line 10000 for more info
1  REM Mandelbrot version for the ][ series (][,][+,][e,//c,][GS)
10 INPUT "Input Pmin, Pmax, Qmin, Qmax: ";PN,PX,QN,QX
20 M = 100
22 MC = 100
25 HGR2
30 DP = (PX - PN) / 279 : DQ = (QX - QN) / 191
40 FOR Y = 0 TO 191: FOR X = 0 TO 279
50 P = PN + X * DP : Q = QN + Y * DQ: K = 0 : X0 = 0 : Y0 = 0
60 XK = X0 * X0 - Y0 * Y0 + P : YK = 2 * X0 * Y0 + Q
70 K = K + 1: X0 = XK : Y0 = YK : R = X0 * X0 + Y0 * Y0
80 IF R > M THEN CL = INT(K / 2) * 4 : GOTO 110
90 IF K = MC THEN CL = 0 : GOTO 110
100 GOTO 60
110 HCOLOR= CL : HPLOT X , Y : NEXT X : NEXT Y
10000 REM _______________________________________________
10010 REM |                                             |
10020 REM |         Mandelbrot Set grapher ][           |
10030 REM |                  by                         |
10040 REM |              Les Barstow                    |
10050 REM |          LMB7421@RITVAX.BITNET              |
10060 REM |    rutgers!rochester!ritcv!ultb!lmb7421     |
10070 REM |    292 Kimball Dr.  Rochester, NY 14623     |
10080 REM |                                             |
10090 REM +---------------------------------------------+
10100 REM ***Freeware.  Distribute, don't sell.