# recta: y= mx +b


# recta que pasa por dos puntos:
# (y-y1) = ( (y1-y2)/(x1-x2) ) * (x - x1)

clear all
clc
close all
hold
plot(0,1,'o')
plot(1,0,'o')
plot(0,0,'o')

xi = rand();
yi = rand();

p1x = cos(pi/3);
p1y = sin(pi/3);

p2x = 0;
p2y = 0;

p3x = 1;
p3y = 0;

n = 1;
N = 3000;

while n <= N
 punto = unidrnd(3);

 if punto == 1
  xi = (p1x + xi) / 2;
  yi = (p1y + yi) / 2;
 end

 if punto == 2
   xi = (p2x + xi) / 2;
   yi = (p2y + yi) / 2;
 end

 if punto == 3
   xi = (p3x + xi) / 2;
   yi = (p3y + yi) / 2;
 end

 plot(xi,yi,'.')
 n += 1;

end