REM Fractal.BAS
REM Copyright (C) 1996 Erik Melin, Melin Software
REM If you use this code, please send an E-mail to melin@mbox302.swipnet.se


SCREEN 13
nx% = 319
ny% = 199
niter% = 100

'last = 1000

REM To zoom in, modify these!
x0 = -1
y0 = -1
x1 = 1
y1 = 1


dx = (x1 - x0) / nx%
dy = (y1 - y0) / ny%

FOR px% = 0 TO nx%
 FOR py% = 0 TO ny%
  cx = x0 + px% * dx
  cy = y0 + py% * dy
  x = 0
  y = 0
  xsqr = 0
  ysqr = 0

  FOR i% = 0 TO niter% - 1
   IF xsqr + ysqr > 4 THEN EXIT FOR
   xsqr = x * x
   ysqr = y * y
   xx = xsqr - ysqr + cx
   yy = (x * y) * 2 + cy
   x = xx
   y = yy
  NEXT i%
  IF i% = niter% THEN i% = 0 ELSE i% = i% MOD 256

  PSET (px%, py%), i%
 NEXT py%
NEXT px%

REM End Of File