function [X Y] = init_circle(N)


for i=1:N
    deg = randi(360);
    rad = deg*pi/180;
    X(i) = cos(rad);
    Y(i) = sin(rad);
end

X = X';
Y = Y';