function [xx,yy]=tree_leaf(N)
%IFS生成“带梗的植物叶子”(建议运行10000次以上)
x=0;y=0;bbs=rand(1,N);
AA=[0,0,0.16,0,0,0;0.85,-2.5/180*pi,0.85,-2.5/180*pi,0,1.6;...
   0.3,49/180*pi,0.34,49/180*pi,0,1.6;0.3,120/180*pi,0.37,-50/180*pi,0,0.44];
figure;
xx=zeros(N,1); yy=zeros(N,1);
for ss=1:N;
   if bbs(1,ss)<=0.005;
       [x,y]=IFS(x,y,AA(1,1),AA(1,2),AA(1,3),AA(1,4),AA(1,5),AA(1,6));
   elseif bbs(1,ss)<=0.805;
       [x,y]=IFS(x,y,AA(2,1),AA(2,2),AA(2,3),AA(2,4),AA(2,5),AA(2,6));
   elseif bbs(1,ss)<=0.9025;
       [x,y]=IFS(x,y,AA(3,1),AA(3,2),AA(3,3),AA(3,4),AA(3,5),AA(3,6));
   else
       [x,y]=IFS(x,y,AA(4,1),AA(4,2),AA(4,3),AA(4,4),AA(4,5),AA(4,6));
   end
   xx(ss)=x; yy(ss)=y;
end
plot(xx,yy,'.g','markersize',2);
set(gcf,'color','w')
axis square off;
title('带梗的植物叶子')
function [xp,yp]=IFS(x,y,r,thita,s,phi,h,k)
xp=r*x*cos(thita)-s*y*sin(phi)+h;
yp=r*x*sin(thita)+s*y*cos(phi)+k;
return