sierpinski

841 days ago by kirubakaran

space = plot(point((0,0))) def etri(length=1,x=0,y=0): p1 = (x , y) p2 = (x+length , y) p3 = (x+length/2 , y+length) l1 = line((p1,p2)) l2 = line((p2,p3)) l3 = line((p3,p1)) return (l1+l2+l3) def sier(l,x,y,i): global space if i==0: return else: space = space + etri(l,x,y) sier(l/2 , x , y , i-1) sier(l/2 , x+l/4 , y+l/2 , i-1) sier(l/2 , x+l/2 , y , i-1) sier(4,0,0,7) show(space) 
       
space = polygon([(0,0,0),(0,0,0),(0,0,0)]) def etri3(l=1,x=0,y=0,z=0): return polygon([(x,y,z),(x+l/2,y+l,z),(x+l,y,z)],opacity=0.75) def sier(l,x,y,i): global space if i==0: return else: space = space + etri3(l,x,y,i) sier(l/2 , x , y , i-1) sier(l/2 , x+l/4 , y+l/2 , i-1) sier(l/2 , x+l/2 , y , i-1) sier(4,0,0,6) show(space)