Sierpinski probabiliste

[P.Fradin]

L'idée m'avait été montrée par Alphonse dans un ancien forum, je viens de retomber dessus aujourd'hui:

C'est un bitmap avec 20 000 points:

Code: Tout sélectionner

[NewBitmap(),
 a:=-4-4*i, b:=4-4*i, c:=4*i, m:=Rand()+i*Rand(),
 Pixel(
     for k from 1 to 20000 do
           x:=3*Rand(),
           if x<1 then m:=(a+m)/2, Scr2Pixel(m), green
           elif x<2 then m:=(b+m)/2, Scr2Pixel(m), blue
           else m:=(c+m)/2, Scr2Pixel(m), red
           fi
     od)
]


L'export au format bmp puis png (bouton capture):

SierpinskiProba.png (6.55 Kio) Consulté 8112 fois

[P.Fradin]

Après le triangle, le tapis:

Code: Tout sélectionner

[
 NewBitmap(),
 a:=-4-4*i, b:=4-4*i, c:=4+4*i,d:=-4+4*i,
 c1:=(a+b)/2, c2:=(b+c)/2, c3:=(c+d)/2, c4:=(d+a)/2,
  m:=Rand()+i*Rand(),
 Pixel(
     for k from 1 to 100000 do
           x:=8*Rand(),
           if x<1 then m:=(2*a+m)/3, Scr2Pixel(m), red
           elif x<2 then m:=(2*c1+m)/3, Scr2Pixel(m), blue
           elif x<3 then m:=(2*b+m)/3, Scr2Pixel(m), red
           elif x<4 then m:=(2*c2+m)/3, Scr2Pixel(m), blue
           elif x<5 then m:=(2*c+m)/3, Scr2Pixel(m), red
           elif x<6 then m:=(2*c3+m)/3, Scr2Pixel(m), blue
           elif x<7 then m:=(2*d+m)/3, Scr2Pixel(m), red
           elif x<8 then m:=(2*c4+m)/3, Scr2Pixel(m), blue
           fi
     od)
]


SierpinskiProba.png (13.73 Kio) Consulté 8104 fois