1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
| def transGlob(L):
import numpy as np
import pylab as py
theta=np.pi/3
w=np.cos(theta)+1j*np.sin(theta)
n=len(L)
Ls=[]
for i in range(0,n-1):
za,z2,z3,z4=transEle(L[i],L[i+1],w)
for z in (za,z2,z3,z4):
Ls.append(z)
Ls.append(L[len(L)-1]) #Le dernier point
return Ls
def transEle(za,zb,w):
h=(zb-za)/3
z2=za+h
z4=z2+h
z3=z2+h*w
return za,z2,z3,z4
import numpy as np
import pylab as pl
L=[0,1]
k=0
for x in range(1,7):
k=k+1
print(k)
L=transGlob(L)
pl.plot(np.real(L),np.imag(L),lw=3,color='black')
pl.show()
pl.close() |
Partager