with(plots):
r:=t->exp(cos(t))-2*cos(4*t)+sin(t/12)^5:
K:=300:for i from 1 to K do ti:=i*50*Pi/K:
curve[i]:=polarplot(r(t),t=0..ti,numpoints=3000*i/K) od:
curve:=display(seq(curve[i],i=1..K),insequence=true):
text:=textplot([0,3,"Butterfly"],color=grey):
display(curve,text,scaling=constrained,tickmarks=[6,6]);
评论