Дракон Хартера-Хайтвея

function [X,Y]=Dracon(Lmax)
% функция, возвращающая изображение дракона
%  -- порЯядок дракона
%Lmax= 'FX';
%Lmax= 'F';
Axiom='FX';
Newf='F';
Newx='X+YF+';
Newy='-FX-Y';
teta=pi/2;
alpha=0;
p=[0;0];
p=Coord(p,Lmax,Axiom,Newf,Newx,Newy,alpha,teta);
% обращение к функции, возвращающей координаты угловых точек дракона
M=size(p,2);
X=p(1:1,1:M);
% создание вектора, содержащего Х-е координаты угловых точек дракона
Y=p(2:2,1:M);
% создание вектора, содержащего Y-е координаты угловых точек дракона
plot(X,Y,Color,k);
set(gca,xtick, Ytick,[]);
set(gca,XColor,w,YColor,w);
function z=Coord(p,Lmax,Axiom,Newf,Newx,Newy, alpha,teta)
% функция, возвращающая координаты угловых точек дракона
Rule=DraconString(Lmax,Axiom,Newf,Newx,Newy, 1, '') ; % задание L-системы
M=length(Rule);
for i=l:M
Tmp=p(1:2, size(p,2) :size(p,2));
if Rule(i)==F % шаг вперед
R=[cos(alpha);sin(alpha)];
R=R/(2^Lmax);
p=cat(2,p,Tmp);
end;
if Rule(i) == '+' % увеличение угла, задающего направление движения
alpha=alpha+teta;
end;
if Rule(i) == '-' % уменьшение угла, задающего направление движения
alpha=alpha-teta;
end;
end;
z=p;
function z1=DraconString(Lmax,Axiom,Newf,Newx,Newy,n,tmp)
% функция, возвращающая L-систему
if n<=Lmax
if n==1
tmp=Axiom;
end;
M=length(tmp);
tmpl= '';
for i=l:M
if tmp(i)==X 
    tmpl=streat(tmpl,Newx); end;
 
if tmp(i)==Y
    tmpl=streat(tmpl,Newy); end;
if not(tmp(i)==F) &&not(tmp(i)==X) &&not(tmp(i)==Y) 
    tmpl=strcat(tmpl,tmp(i));
end;
tmp=tmpl;
n=n+l;
tmp=DraconString(Lmax,Axiom,Newf,Newx,Newy ,n, tmp);
end;
end;
z1=tmp;