% OCRTAN_KROG 	O"crtan krog
%	[S,r]=OCRTAN_KROG(T) vrne sredi"s"ce in radij trikotniku T
%	o"crtanega kroga. T je matrika 2x 3, 
%       kjer prva vrstica predstavlja
%	x koordinate. S=[x,y], r je "stevilo 

function [S,r]=ocrtan_krog(T);
  T=T'; %da uskladimo dimenzije
  %dolocimo koordinati razpolovisc stranic a in c
  a=T(2,:)-T(1,:);
  c=T(3,:)-T(1,:);
  Sa=T(1,:)+1/2*a;
  Sc=T(1,:)+1/2*c;
  %dolocimo pravokotnici na a in c skozi Sa in Sb v implicitni obliki
  %najprej pogledamo, ali je a vzporedna x osi
  if T(1,2)==T(2,2)
    Pa=[1,0,-Sa(1,1)];
  else
    %dolocimo smerni koeficient pravokotnice
    ka=-(T(2,1)-T(1,1))/(T(2,2)-T(1,2));
    Pa=[ka,-1,-ka*Sa(1,1)+Sa(1,2)];
  end
  
  %enako se za pravokotnico na c
  if T(1,2)==T(3,2)
    Pc=[1,0,-Sc(1,1)];
  else
    %dolocimo smerni koeficient pravokotnice
    kc=-(T(3,1)-T(1,1))/(T(3,2)-T(1,2));
    Pc=[kc,-1,-kc*Sc(1,1)+Sc(1,2)];
  end
  
  S=presek_premic(Pa,Pc); %izracunamo presecisce pravokotnic
  %dolocimo se r kot razdaljo od S do npr. T1
  vr=S-T(1,:);
  r=sqrt(vr*vr');
    
    
