% Vhodni podatek: 2x3 matrika, katere stolpci predstavljajo koordinate
% oglisc trikotnika
O = [1 3 0;2 1 -1]; 
A = O(:,1);
B = O(:,2);
C = O(:,3);

% OCRTANA KROZNICA
[T1,s1] = simetrala_daljice(A,B);
[T2,s2] = simetrala_daljice(A,C);

p = presek_premic(T1, T2, s1,s2); %sredisce

R = norm(p-A); %polmer

% VCRTANA KROZNICA
[Q1,q1] = simetrala_kota(A,B,C);
[Q2,q2] = simetrala_kota(B,C,A);

q = presek_premic(Q1, Q2, q1,q2); %sredisce

a = norm(C-B); %dolzine stranic
b = norm(A-C);
c = norm(A-B);
s = (a+b+c)/2; %polovicni obseg
ploscina = sqrt(s*(s-a)*(s-b)*(s-c)); %ploscina
r = ploscina/s; %polmer: r = ploscina trikotnika / polovicni obseg

% PARAMETRIZACIJA KROZNIC ZA RISANJE
t = linspace(0,2*pi, 100);

x_o = p(1) + R*cos(t); %ocrtana
y_o = p(2) + R*sin(t);

x_v = q(1) + r*cos(t); %vcrtana
y_v = q(2) + r*sin(t);

% narisemo obe kroznici in njuni sredisci
hold on
plot(x_o,y_o,x_v,y_v)
plot(p(1),p(2), 'x')
plot(q(1),q(2), 'x')

% narisemo se trikotnik
xx = O(1,:);
yy = O(2,:);
xx(4) = xx(1);
yy(4) = yy(1);
line(xx,yy)
axis equal
hold off