大家看看这个程序运行的结果是否漂亮？

minjiecow

tree[set_, m_] :=
  Module[{\[Theta], mz, my, temp, k, i, p1, p2, n, a, b, c, d, e, f},
   \[Theta] = Pi/6;
   mz = {{Cos[\[Theta]], -Sin[\[Theta]]}, {Sin[\[Theta]],
      Cos[\[Theta]]}};
   my = {{Cos[-\[Theta]], -Sin[-\[Theta]]}, {Sin[-\[Theta]],
      Cos[-\[Theta]]}};;
   temp = set;
   p1 = {};(*存放树干*)
   p2 = {};(*存放树枝*)
   For[k = 1, k <= m, k++,
    n = Length[temp];
    p2 = {};
    For[i = 1, i <= n, i++,
     L = Part[temp, i];
     a = Part[L, 2];
     b = Part[L, 1];
     c = 1/3*a + 2/3*b;
     d = c + mz.(a - b)*(1/3);
     e = c + my.(a - b)*(1/3);
     AppendTo[p1, {Brown, Thickness[0.035/k], Line[{b, c}]}];
     AppendTo[p2, {c, a}];
     AppendTo[p2, {c, d}];
     AppendTo[p2, {c, e}];
     ];
    temp = p2;
    ];
   
   Show[Graphics[{
      p1,
      {RGBColor[0.109804, 0.423529, 0.172549], Line[p2]}
      }]
    ]
   ];
tree[{{{0., 0.}, {0., 1.}}}, 7]

HyperGroups

多了一个封号，
漂亮，楼主再接再励，弄个三维的漂亮的树

熔雨

怎么弄的啊，不明白。

chyanog

1. ctree = Compile[{{ls, _Real, 2}},
   
2.    Module[{a, b},
   
3.     {a, b} = ls;
   
4.     {{a, (a + b)/2}, {(a + b)/2, b},
   
5.      {(a + b)/2, {{1, 1}, {-1, 1}} . (b - a)/2.82 + (a + b)/2}, {(3*a + b)/4,
   
6.       ({{1.73, -1}, {1, 1.73}}.(a + 3 b) +   {{.26, 1}, {-1, .26}}.(3 a + b))/8}}],
   
7.    RuntimeAttributes -> {Listable}];
   
8. 
   
9. With[{p = Flatten[Nest[ctree, N@{{{0, 0}, {0, 1}}}, #], #] &[5]},
   
10. Graphics[Line[#,
    
11.      VertexColors -> (ColorData["AvocadoColors", #] & /@ #[[All,   2]])] & /@ Rescale[p]]
    
12. ]
    
13. 
    

复制代码