function trirandom(n)
a=[0,0];
b=[1,0];
c=[0.5,sqrt(3)/2];
x=[0.5,0.5];
figure;
axis([-1,2,-1,2]);
s(1,1)=x(1);
s(1,2)=x(2);
for i=2:n
r=rand;
if (r<0.333333)
x=(a+x)/2;
elseif (r<0.666666)
x=(b+x)/2;
else
x=(c+x)/2;
end
s(i,1)=x(1);
s(i,2)=x(2);
end
plot(s(:,1),s(:,2),'.');

