domingo, 7 de febrero de 2016

Sierpinski triangle. Chaos game


###Python 3.4

import numpy as np
import matplotlib.pyplot as plt
import math as m
np.random.seed(5)

blueX=0
blueY=0
greenX=50
greenY=50
redX=100
redY=0
a=plt.plot(0,0,"o",color='b')
plt.plot(100,0,"o",color='r')
plt.plot(50,100*m.sin(30*m.pi/180),"o",color='g')

iniX=20
iniY=20
plt.plot(iniX,iniY,'.')
for i in range(50000):
    x=np.random.randint(1,4)
    if (x==1):
        iniX=(iniX+blueX)/2
        iniY=(iniY+blueY)/2
    elif x==2:
        iniX=(iniX+greenX)/2
        iniY=(iniY+greenY)/2
    else:
        iniX=(iniX+redX)/2
        iniY=(iniY+redY)/2
    plt.plot(iniX,iniY,'.',markersize=0.5)  
       

plt.axis([-1,101,-1,51])

F = plt.gcf()


F.savefig('Siperpinski.jpeg', dpi = (200),bbox_inches='tight')



No hay comentarios:

Publicar un comentario