SET MODE "color"
SET WINDOW 0,1026,0,750
SET BACKGROUND COLOR "white"
LET x=rnd/3
LET y=rnd/3
LET a=0
LET b =0
LET s1=310
LET s2 =s1*750/1026
RANDOMIZE
PRINT"  double diaxial Cantor set  I.F.S."
PRINT" BY R.L.BAGULA 20 Sept 1998 ©"
FOR n= 1  TO  2560000
    LET a =RND
    IF  a<= 1/4   THEN
       LET x1=x*y/9
       LET y1=(x/3+2/3)*(y/3+2/3)
    END IF
    IF  a<= 2/4 AND a>1/4 THEN
       LET x1=(x/3+2/3)
       LET y1= y/3
    END IF
    IF  a<= 3/4 AND a>2/4 THEN
       LET x1= x/3
       LET y1= y/3+2/3
    END IF
    IF  a<= 1 AND a>3/4 THEN
       LET x1=(x/3+2/3)*(y/3+2/3)
       LET y1= x*y/9
    END IF
    SET COLOR 256-32*mod(4*a,4)
    LET x=x1*sqr(2.5)
    LET y=y1*sqr(2.5)
    IF n>10 THEN PLOT 1026/2+s1*x-250,750/2+s2*y-200
NEXT n
END




PRINT "input Sierpinski number 3 to 12: try 4"
INPUT m
SET MODE "color"
SET WINDOW 0,1026,0,750
SET COLOR MIX (1) 0,0, 0
SET BACKGROUND COLOR "white"
LET x=1
LET y=1
LET c=0
LET s1=175
LET s2 =(s1)*750/1026
DIM a(25),b(25)
FOR i=1 to m
    LET a(i)=cos(2*pi*i/m)
    LET b(i)=sin(2*pi*i/m)
NEXT i
PRINT "  product -switch SIERPINSKI I.F.S. of first kind "
PRINT "  BY R.L.BAGULA 23 Sept 1998© copy rights reserved"
RANDOMIZE
PRINT "  M=";m
FOR n= 1  TO  2000000
    LET c =RND
    LET l=1+int(c*m)
    REM 50% switch between product and ordinary Sierpinski of first kind
    IF rnd>.5 then
       LET x1=(x/2+a(l))*(y/2+a(l))
       LET y1=(x/2+b(l))*(y/2+b(l))
    ELSE
       LET x1=(x/2+a(l))
       LET y1= (y/2+b(l))
    END IF
    LET x=x1
    LET y=y1
    SET COLOR 255
    IF n>10 THEN PLOT 1026/2+s1*x-150  ,750/2+s2*y-100
NEXT n
END