Shevchenko / 2sem / LR10_1

clc
clear
sx=[1, -1, -1, 1, 1]
sy=[1, 1, -1, -1, 1]

function[polar] = anglep(x1, y1, x2, y2);
dx = (x2-x1);
dy = (y2-y1);
r = sqrt((dx^2)+(dy^2));
sns = dy/r;
ang = asin(sns); %Находим значение относительного положения точек по координате игрик в радианах
polar = ang;
if dx<0 && dy>0
polar = pi()-ang; %Находим полярные углы из положения на окружности
end
if dy<0 && dx>0
polar = 2*pi()+ang;
end
end

for I=1:4; %Выполняем для каждой пары точек
x1=sx(I);
y1=sy(I);
x2=sx(I+1);
y2=sy(I+1);
p=anglep(x1, y1, x2, y2);
p
end
hp=compass(sx, sy); %Построение кругового графика