function triangoli(P,Q,R,n)

% Disegna i triangoli T, M(T), M(M(T)), ..., M^n(T) dove T e' il
% triangolo di vertici PQR e, dato un triangolo R, M(R) e' il
% triangolo che ha vertici nei punti medi ti R.

% Esempio:
%   >> triangoli([-1,1,1], [1,3,0], [2,2,3], 5);

if n<=0
    disp('n deve essere un intero positivo');
    return
end

hold on;
for i=1:n
    % Visualizza il triangolo.
    M = [P; Q; R; P].';
    plot(M(1,:), M(2,:));
    % Calcola i vertici del nuovo triangolo.
    P1=(P+Q)/2;
    Q1=(Q+R)/2;
    R1=(R+P)/2;
    P=P1;
    Q=Q1;
    R=R1;
end
hold off;