
###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