绘制递归树

MeteoInfo

脚本：

1. def f(x, y, d, R, n, N):
   
2.     if n < N:
   
3.         b = zeros(3)
   
4.         if n > 3:
   
5.             b = [.5,.2,.6]
   
6. 
   
7.         X = R * sin(np.deg2rad(d)) + x
   
8.         Y = R * cos(np.deg2rad(d)) + y
   
9.         plot(X, Y, color=b, linewidth=np.linalg.norm(R)/15)
   
10.         A = array([0,-48,-25,43,31])
    
11.         for i in range(1, 5):
    
12.             r = (i + 1) % 2
    
13.             L = .58 + r * .07
    
14.             f(X, Y, d+A, L*R, n+1, 7-r)
    
15. 
    
16. ax = axes(outerposition=[0,0,1,1], axis=False)
    
17. f(0, 0, array([0,0,-5.76,0,6.58,0]), array([0,55,99,55,105,55]), 0, 1)

def f(x, y, d, R, n, N):<br />     if n < N:<br />         b = zeros(3)<br />         if n > 3:<br />             b = [.5,.2,.6]<br /> <br />         X = R * sin(np.deg2rad(d)) + x<br />         Y = R * cos(np.deg2rad(d)) + y<br />         plot(X, Y, color=b, linewidth=np.linalg.norm(R)/15)<br />         A = array([0,-48,-25,43,31])<br />         for i in range(1, 5):<br />             r = (i + 1) % 2<br />             L = .58 + r * .07<br />             f(X<i>, Y<i>, d+A<i>, L*R, n+1, 7-r)<br /> <br /> ax = axes(outerposition=[0,0,1,1], axis=False)<br /> f(0, 0, array([0,0,-5.76,0,6.58,0]), array([0,55,99,55,105,55]), 0, 1)

VV09

厉害厉害~~膜拜大神！