Mort de Benoît Mandelbrot

' Adaptation en Panoramic de :
' http://snipplr.com/view/15632/gwbasic-mandelbrot-set-fractal-explorer/
DIM NbIteration,OX,OY,OU,D,X,Y,X1,Y1,X2,Y2,IT,YY,XX,A
dim p(9,3),i,j

LABEL Calc
LABEL Palette: gosub Palette

LABEL bClick

DIM W,H: W=400 : H=200
WIDTH 0,W+10 : HEIGHT 0,H+80

DIM Image:Image=number_objects+1:PICTURE Image
WIDTH Image,w:Height Image,H
2D_TARGET_IS Image


DIM Gauche:Gauche=number_objects+1:BUTTON Gauche: HINT Gauche,"Gauche"
CAPTION Gauche,"G": ON_CLICK Gauche,bClick : TOP Gauche,H+10
LEFT Gauche, 10:WIDTH Gauche,20
DIM Droite:Droite=number_objects+1:BUTTON Droite : HINT Droite,"Droite"
CAPTION Droite,"D": ON_CLICK Droite,bClick : TOP Droite,TOP(Gauche)
LEFT Droite, LEFT(Gauche)+WIDTH(Gauche)+5:WIDTH Droite,WIDTH(Gauche)
DIM Haut:Haut=number_objects+1:BUTTON Haut: HINT Haut,"Haut"
CAPTION Haut,"H": ON_CLICK Haut,bClick : TOP Haut,TOP(Gauche)
LEFT Haut, LEFT(Droite)+WIDTH(Droite)+5:WIDTH Haut,WIDTH(Gauche)
DIM Bas:Bas=number_objects+1:BUTTON Bas: HINT Bas,"Bas"
CAPTION Bas,"B": ON_CLICK Bas,bClick : TOP Bas,TOP(Gauche)
LEFT Bas, LEFT(Haut)+WIDTH(Haut)+5:WIDTH Bas,WIDTH(Gauche)
DIM ZoomIn:ZoomIn=number_objects+1:BUTTON ZoomIn: HINT ZoomIn,"Zoom In"
CAPTION ZoomIn,"ZI": ON_CLICK ZoomIn,bClick : TOP ZoomIn,TOP(Gauche)
LEFT ZoomIn, LEFT(Bas)+WIDTH(Bas)+5:WIDTH ZoomIn,WIDTH(Gauche)
DIM ZoomOut:ZoomOut=number_objects+1:BUTTON ZoomOut: HINT ZoomOut,"Zoom Out"
CAPTION ZoomOut,"ZO": ON_CLICK ZoomOut,bClick : TOP ZoomOut,TOP(Gauche)
LEFT ZoomOut, LEFT(ZoomIn)+WIDTH(ZoomIn)+5:WIDTH ZoomOut,WIDTH(Gauche)


DIM TimerCalc:TimerCalc=number_objects+1:TIMER TimerCalc
TIMER_INTERVAL TimerCalc,20

REM Subroutine -- Set up Palette
DATA 0,0,0      : ' black
DATA 227,91,216  : ' magenta
DATA 0,0,255    : ' blue
DATA 0,255,0    : ' green
DATA 224,255,255 : ' cyan
DATA 255,0,0    : ' red
DATA 255,140,0  : ' orange
DATA 255,255,0  : ' yellow
DATA 255,255,255 : ' white

TO_FOREGROUND 0
NbIteration=40
OX=2:OY=1:OU=100
OU=OU*2:OX=OX-160/OU:OY=OY-100/OU
GOSUB Palette : ' init de la palette de couleur
ON_TIMER TimerCalc,Calc


END

bClick:
  select (NUMBER_Click-Gauche)
  case 0
    OX=OX+80/OU : ' gauche
    X=0:Y=0
  case 1
    OX=OX-80/OU : ' Droite
    X=0:Y=0
  case 2
    OY=OY+50/OU : ' Haut
    X=0:Y=0
  case 3
    OY=OY-50/OU : ' Bas
    X=0:Y=0
  case 4
    OU=OU*2:OX=OX-160/OU:OY=OY-100/OU : ' ZoomIn
    X=0:Y=0
  case 5
    OX=OX+160/OU:OY=OY+100/OU:OU=OU*.5 : ' ZoomOut
    X=0:Y=0
  end_select
return


Calc:
  IF X=(W-1) THEN RETURN
  IF Y=(H-1)
    X=X+1:Y=0
  ELSE
  IF Y<(H-1)
    Y=Y+1
  END_IF
  END_IF
  X1 = X/OU-OX : XX=X1
  Y1=Y/OU-OY:YY=Y1
  X2 = XX*XX
  Y2 = YY*YY
  FOR IT=0 TO NbIteration
    IF (X2+Y2) > 4 THEN EXIT_FOR
    YY = 2*XX*YY + Y1
    XX = X2 - Y2 + X1
    X2 = XX*XX
    Y2 = YY*YY
  NEXT IT
  A=IT-INT(IT/5)*5
  2d_pen_color p(a,0),p(a,1),p(a,2)
  2d_point x,y
  CAPTION 0,"Mandelbrot Explorer"+" ("+str$(160/OU-OX)+","+str$(100/OU-OY)+") Zoom:"+str$(OU/100)+" X:"+STR$(X)+" Y:"+STR$(Y)+" NbIter:"+str$(NbIteration)
return


Palette:
  RESTORE
  for i=0 to 8
    for j=0 to 2
      read p(i,j)
    next j
  next i
RETURN