SET MODE "color"
SET WINDOW 0,1920,0,1024
SET BACKGROUND COLOR "white"
LET x=.3
LET y=.2
LET a=0
LET b =0
LET s1=550
LET s2 =s1*1024/1920
RANDOMIZE
PRINT" KISSING CIRCLES OF A RATIONAL JULIA: ROGER'S RINGS"
PRINT" "
PRINT" BY R.L.BAGULA 21 Feb 2002 ©"
FOR n= 1  TO  2560000
    LET a =RND
    LET b=RND
    LET R=SQR(X^2+Y^2)
    LET W=ANGLE(X,Y)
    LET x2=R/(R+0.5)-X
    LET y2=0.5/(R+0.5)-Y
    LET r2=sqr(x2*x2+y2*y2)
    LET w2=angle(x2,y2)
    IF a <= 1/3 THEN
       LET x1=-r2^.5*cos(.5*w2)
       LET y1=-r2^.5*sin(.5*w2)
    END IF
    IF  a<= 2/3 AND a>1/3 THEN
       LET x1=r2^.5*cos(.5*w2)
       LET Y1=r2^.5*sin(.5*w2)
    END IF
    IF  a<= 1 AND a>2/3 THEN
       LET W=ANGLE(X,Y)
       LET x1= SIN(2*W)/2
       LET y1= COS(2*W)/2
    END IF

    SET COLOR 255-64*int(mod(3*a,3))
    LET x=x1
    LET y=y1
    IF n>10 THEN PLOT 1920/2+s1*x,1024/2+s2*y
NEXT n
END