分成六部分，同上題作法，做六次後再組合。

> restart;

> with(plots):

Warning, the name changecoords has been redefined

>

> m:=[[0,0],[1,0],[cos(Pi/3),sin(Pi/3)],[0,0]];

>

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a1:=plot(w,axes=none,scaling=constrained,color=brown):

> m:=[[0,0],[cos(2*Pi/3),sin(2*Pi/3)],[cos(Pi/3),sin(Pi/3)],[0,0]];

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a2:=plot(w,axes=none,scaling=constrained,color=brown):

> m:=[[0,0],[cos(2*Pi/3),sin(2*Pi/3)],[cos(Pi),sin(Pi)],[0,0]];

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a3:=plot(w,axes=none,scaling=constrained,color=brown):

> m:=[[0,0],[cos(4*Pi/3),sin(4*Pi/3)],[cos(Pi),sin(Pi)],[0,0]];

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a4:=plot(w,axes=none,scaling=constrained,color=brown):

> m:=[[0,0],[cos(4*Pi/3),sin(4*Pi/3)],[cos(5*Pi/3),sin(5*Pi/3)],[0,0]];

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a5:=plot(w,axes=none,scaling=constrained,color=brown):

> m:=[[0,0],[1,0],[cos(5*Pi/3),sin(5*Pi/3)],[0,0]];

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> a6:=plot(w,axes=none,scaling=constrained,color=brown):

> display(a1,a2,a3,a4,a5,a6);

>